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EXECUTIVE  SUMMARY 


This  is  the  final  report  of  the  O.N.R.  sponsored  research  program 
entitled  "Study  of  Residual  Stresses  and  Distortion  in  Structural 
Weldments  in  High  Strength  Steels".  All  five  tasks  of  this  program, 
which  was  initiated  on  December  1,  1977,  have  been  completed  as  proposed. 

Task  1  deals  with  butt  welds  in  thick  HY-130  steel  plates.  Experi¬ 
mental  measurements  of  temperatures,  thermal  strains,  residual  stresses 
and  distortions  were  performed  for  a  variety  of  welded  specimens  during 
the  first  three  years  of  the  program.  Fracture  experiments  were  also 
performed  on  simply  restrained  welded  specimens.  Analytical  studies 
were  conducted,  parallel  to  the  experiments,  and  were  completed  during 
the  fourth  year  of  the  program.  Specif ical1y,  heat  transfer  analyses 
were  performed  using  both  a  newly  developed  analytical  solution  based  on 
a  moving  finite  heat  source  model,  and  the  multipurpose  nonlinear  heat 
transfer  finite  element  computer  program  ADINAT,  which  was  modified  to 
take  into  account  latent  heat  effects. 

Calculations  of  the  thermal  strain,  transient  stress,  and  residual 
stress  distributions  developed  during  the  multipass  GMA  welding  of  thick 
plates,  were  also  performed  using  the  thermal-elastic-plastic  constitutive 
model  of  the  multipurpose  finite  element  program  ADINA.  The  program  was 
enhanced  by  incorporating  the  phase  transformation  strains  developed 
during  the  steel  allotropic  phase  changes. 

A  hybrid  crack  element  was  also  developed,  aimed  at  facilitating  an 
inexpensive  evaluation  of  the  fracture  characteristics  of  weldments  in 
the  presence  of  a  residual  stress  field.  Using  this  element,  linear 
elastic  finite  element  analyses  of  the  fracture  experiments  were  performed. 
Both  the  crack  opening  displacement  and  the  strain  energy  released  during 
the  propagation  of  a  crack  in  a  residual  stress  field  were  calculated. 

Finally,  results  from  all  of  the  afc rementioned  analyses  were 
compared  with  data  obtained  from  experimental  measurements.  Good  correla¬ 
tion  was  generally  observed. 


Task  2  deals  with  girth  welds  in  HY-130  unstiffened  cylindrical 
shells.  Temperatures  and  strains  during  seam  and  girth  welding  of  the 
cylinders  were  first  measured.  The  axial  and  circumferential  residual 
stress  distributions  were  then  determined  and  compared  with  data  from 
other  investigators  and  analytical  predictions.  Task  2  was  completed 
during  the  third  year  of  this  program. 

Task  3,  aimed  at  developing  information  on  weldments  in  bronze,  was 
carried  in  its  entirety  during  the  fourth  year.  A  literature  survey  was 
first  conducted  to  obtain  material  properties  on  nickel-aluminum  bronze. 
Experimental  studies,  consisting  of  bead-on-plate  and  edge  welds  on  one 
quarter  inch  thick  plates,  were  then  performed  to  determine  the  resulting 
temperature,  strain,  and  displacement  changes.  The  experimental  results 
were  finally  compared  with  baseline  data  to  study  the  effects  of  restraint 
block  welding,  preheating,  postheating,  and  grinding  on  the  final  distor¬ 
tion. 

In  Task  4,  the  mechanisms  and  the  effectiveness  of  thermal  stress 
velieving  treatments  for  weldments  are  studied.  A  literature  survey 
covering  both  thermal  and  mechanical  stress  relieving  treatments  was 
initially  performed.  Then  the  M.I.T.  one-dimensional  program  was 
properly  modified  to  calculate  strains  and  stresses  during  and  after 
uniform-temperature  heat  treatments  (taking  into  account  stress  relaxation 
due  to  creep).  The  analytical  predictions  were  verified  by  experiments 
performed  on  304  stainless  steel  plates.  Further  analysis  was  performed 
to  evaluate  the  effectiveness  of  line  heating  with  a  flame.  Finally, 
experiments  with  HY-130  specimens  were  conducted  during  the  six-month 
extension  of  the  program. 

Under  Task  5,  manuals  have  been  prepared  for  the  computer  programs 
related  to  welding  that  have  been  developed  or  improved  during  the  course 
of  the  present  research  program.  This  task  was  also  completed  during  the 
final  six-month  extension. 

All  in  all,  therefore,  the  proposed  program  of  research  was  carried 
out  as  originally  proposed  but  with  a  delay  of  approximately  six  months. 
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This  final  report,  issued  at  the  end  of  the  program  extension,  covers  all 
five  tasks.  However,  since  only  a  brief  summary  of  the  work  on  the  first 
three  tasks  is  presented  here,  this  report  should  be  considered  as  a 
companion  to  the  three  technical  progress  reports  previously  published. 
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1 .  INTRODUCTION 

The  objective  of  this  research  program,  which  started  on  December  1, 
1977,  is  to  experimentally  and  analytically  study  temperatures,  thermal 
strains,  residual  stresses  and  distortion  in  structural  weldments  in 
high-strength  steel. 

This  program  includes  (1)  generation  of  experimental  data,  and  (2) 
development  of  analytical  systems.  HY-130  is  the  primary  material  to  be 
investigated;  a  limited  number  of  experiments  are,  however,  to  be  conducted 
using  low-carbon  steel  specimens.  Experiments  are  to  be  made  using  the 
multipass  gas  metal  arc  (GMA)  and  electron  beam  (EB)  one-pass  welding 
processes . 

More  specifically, 

Materials:  HY-130  steel  (with  some  control  specimens  in  low-carbon 
steel) 

Thickness:  1  inch  (with  some  specimens  1/2"  thick) 

Welding  processes:  Electron  beam  (EB)  and  multipass  gas  metal  arc 

(GMA)  processes 

Specimen  geometries:  The  following  types  of  specimens  are  considered 

a.  Butt  welds  in  thick  plates 

(1)  Unrestrained  butt  welds 

(2)  Simple  restrained  butt  welds 

(3)  Restrained  welds  simulating  practical  weldments 

b.  Girth  welds  of  cylindrical  shells 

(4)  Girth  welding  along  the  groove  of  a  cylindrical  shell 

(5)  Girth  welding  between  unstiffened  cylindrical  shells 

The  work  to  be  performed  has  been  divided  into  two  main  tasks: 

Task  1:  Research  on  butt  welds  in  thick  plates 

Task  2:  Research  on  girth  welds  in  cylindrical  shells 
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The  program  was  initially  proposed  to  cover  a  three-year  period. 

Table  1.1  illustrates  the  tasks  and  phases  of  this  study.  The  planned 
progress  for  each  task  appears  in  the  original  proposal,  dated  July  1977, 
where  also  a  more  detailed  description  of  each  task  can  be  found. 

Towards  the  end  of  the  initial  three-year  period  a  proposal,  dated 
November  1980,  was  sent  to  O.N.R.  requesting  a  one-year  extension  of  the 
research  program  to  November  30,  1981.  Under  this  proposal,  which  was 
accepted,  work  on  three  additional  tasks  would  be  done  as  follows  (Table 
1.2): 

Task  3:  Development  of  information  on  weldments  in  bronze 

Task  4:  Research  on  thermal  stress  relieving  of  weldments 

Task  5:  Development  of  improved  computer  programs 

For  the  most  part  the  program  has  been  carried  out  as  originally 
proposed.  For  the  completion  of  the  final  two  tasks,  however,  a  further 
six-month  extension  was  considered  necessary  and  was  requested  in  a 
proposal  dated  October  1981. 

Work  on  all  five  tasks  has  been  completed  by  May  31,  1982  and  is 
reported  here. 

Three  progress  reports,  dated  October  10,  1978,  August  31,  1979, 

July  15,  1980  and  three  technical  progress  reports,  dated  November  30, 

1979,  November  30,  1980  and  November  30,  1981  (covering  the  work  performed 
during  the  first  two,  third  and  fourth  years  respectively)  have  already 
been  issued. 

This  final  report,  issued  at  the  end  of  the  program  extension,  covers 
all  five  tasks.  However,  since  details  on  the  first  three  tasks  have 
already  been  given  in  the  three  technical  progress  reports,  only  a 
summary  will  be  presented  here.  Two  recent  papers  reviewing  the  work  on 
Task  1  are  also  included  in  Appendices  A  and  B.  Tasks  4  and  5,  which  were 
completed  during  the  six-month  extension  of  the  program  will  be  covered 
in  detail  in  this  report.  Therefore. this  final  report  should  be  considered 
as  a  companion  to  the  aforementioned  three  technical  progress  reports. 
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TABLE  1.1 


TASKS  AND  PHASES  OF  THE  PROGRAM  INITIALLY  PROPOSED 


First 

Year 


'Second 


Third 

Year 


Task  1:  Thick  Plate 


Develop  details  of  research  plan 


1.2  Experiments  on  unrestrained 
butt  welds 


1.3  Analysis  of  data  obtained  in 

1.2 


1.4  Develop  details  of  research  plan 


1.5  Measurement  of  residual  stresses 
in  specimens  made  in  1.2 


.6  Experiment  on  simple  restrained 
butt  welds 


1.7  Analysis  of  data  obtained  in 

1.5  and  1.6 


1.8  Develop  details  of  research  plan 


1.9  Experiment  on  restrained  cracking 
test  specimens 


1.10  Measurement  of  strain  energy 
release  on  specimens  made  in 

1 . 6  and  1 . 9 


1.11  Analysis  of  data  obtained  in 
1.9  and  1.10 


Preparation  of  the 


Information  on  residual  stresses 
and  strain  energy  release  in 
specimens  used  for  weld  cracking 
and  stress-corrosion  cracking 


Task  2:  Cylindrical  Shell 


2.1  Develop  details  of  research 
plan 


2.2  Experiment  on  girth  welding 
along  groove  of  a  cylindrical 
shell 


2.3  Analysis  of  data  obtained  in 

2.2 


2.4  Develop  details  of  research 
plan 


2.5  Measurement  of  residual 

stresses  in  specimens  made  in 

2.2 


2.6  Experiment  on  butt  welds 
between  unstiffened 
cylindrical  shells 


2.7  Analysis  of  data  obtained  in 
in  2.5  and  2 . 6 


2.8  Develop  details  of  research 
plan 


2.10  Measurement  of  residual 

stresses  in  specimens  made 
in  2.6 


2.11  Analysis  of  data  obtained 
in  2.10 


Final  Report 


Information  on  residual 
stresses  and  distortion  of 
welded  cylindrical  shells 


TABLE  1.2 


TASKS  TO  BE  COMPLETED  IN  THE  18-MONTH 
EXTENSION 


Task  3:  Bronze  Weldments 


3.1  Literature  Review 

3.2  Generation  of  Experimental  Information 


Task  4:  Thermal  Stress  Relieving 

4.0  Literature  Review 

4.1  Initial  Analysis  and  Verification  Tests 
(304  Stainless  Steel) 

4.2  Further  Analysis  and  Experiments 
(HY-130  Steel) 

Task  5:  Improved  Computer  Programs 


5.1  Development  of  Improved  Manuals  of  Independent 
Welding  Programs 

5.2  Development  of  Manuals  of  Programs  which  are 
Compatible  with  the  ADINA  Program 
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2.  WORK  ACCOMPLISHED  ON  THE  FIRST  THREE  TASKS 

2 . 1  Work  Accomplished  on  Task  1  -  Thick  Plates 

2.1.1  General  Status 

The  objective  of  Task  1  of  this  research  program  is  to  experimental¬ 
ly  and  analytically  study  temperatures,  thermal  strains,  residual  stresses 
and  distortion  in  thick  butt  welded  high  strength  steel  (specirically 
HY-130)  plates.  All  eleven  steps  (see  Table  1.1)  of  this  task  have  been 
completed  successfully  as  originally  planned. 

Most  of  the  experimental  work  has  been  completed  during  the  first 
three  years  of  the  research  program.  Details  of  the  experimental  proce¬ 
dures  and  the  obtained  results  can  be  found  in  the  first  two  technical 
progress  reports  to  O.N.R.  However,  the  final  two  steps  of  Task  1  were 
completed  during  the  fourth  year,  and  were  covered  in  the  third  technical 
progress  repor: 

The  anal  vs  is  of  the  experimental  data  and  the  development  of  computer 
programs  was  also  completed  during  the  fourth  year.  Again,  details  are 
provided  in  the  third  technical  progress  report.  In  what  follows,  the 
work  performed  on  Task  1  is  briefly  summarized.  A  comprehensive  overview 
is  also  provided  in  the  paper: 

Papazogiou,  V..J.,  Masubuchi,  K.  ,  Gonsalves,  E.  ,  and  Imakita,  A., 
"Residual  Stresses  Due  to  Welding:  Computer-Aided  Analysis  of  their 
Formation  and  Consequences",  presented  at  the  SNAME  1982  Annual 
Meeting,  New  York,  N.Y.,  November  18-20,  1982, 

which  is  included  in  Appendix  B  of  this  report. 

2.1.2  F.xpe  r  imental  Program 

The  following  number  of  experiments,  classified  according  to  types 
of  measurements  taken,  specimen  configuration,  and  welding  process  used, 
were  performed  during  the  first  four  years  of  this  research  project 
(December  1,  1977  to  November  30,  1981). 
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A.  Measurement  of  temperature  and  thermal  strain  distributions  during 
welding. 

A.l  Unrestrained  specimens 

A. 1.1  Gas  Metal  Arc  Welding 

(a)  One  24"  x  24"  x  1"  SAE  1020  specimen 

(b)  Four  24"  x  24"  x  1"  HY-130  specimens 
A.  1.2  Electron  Beam  Welding 

(a)  Two  24"  x  24"  x  1"  SAE  1020  specimens 

(b)  Two  24"  x  24"  x  1"  HY-130  specimens 

(c)  One  24  5/16"  x  26"  x  1"  HY-130  specimen 
A. 2  Clamped  specimens 

A. 2.1  Laser  Beam  Welding 

(a)  Two  11"  x  14  3/4"  x  1"  HY-130  specimens 
A. 3  Simply  restrained  specimens^ 

A. 3.1  Gas  Metal  Arc  Welding 


(a) 

Two 

30" 

X 

30” 

X 

1/2 

"  SAE  1020  specimens 

(b) 

Two 

30" 

X 

30" 

X 

1" 

SAE  1020  specimens 

(c) 

Two 

30" 

X 

30" 

X 

2" 

SAE  1020  specimens 

(d) 

Two 

30" 

X 

30" 

X 

7/8 

"  HY-130  specimens 

(3) 

One 

30" 

X 

30" 

X 

2” 

HY-130  specimen 

A. 4  Restrained  cracking  (window)  test  specimen 

A. 4.1  Gas  Metal  Arc  Welding 

(a)  One  15"  x  24"  x  2"  HY-130  specimen  on  a 
30"  x  30"  x  3"  SAE  1020  base  plate 

^These  specimens  had  H-type  slits  of  different  dimensions  to  simulate 
various  degrees  of  restraint.  The  degree  of  restraint  was  calculated 
numerically  (see  Section  2.7.1  of  Second  Technical  Progress  Report  to 
O.N.R.,  November  30,  L980) . 
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B.  Measurement  of  triaxial  distribution  of  residual  stresses 
B.l  Electron  Beam  welded  specimens 

(a)  One  24"  x  24"  x  1"  SAE  1020  specimen  (see  number  A. 1.2a) 

(b)  One  24  5/6"  x  26"  x  1"  HY-130  specimen  (see  number  A. 1.2c) 

B. 2  Gas  Metal  Arc  welded  specimens 

(a)  One  18"  x  32"  x  7/8"  HY-130  specimen 

(b)  One  18"  x  32"  x  1"  HY-130  specimen 

C.  Measurement  of  strain  energy  release  during  crack  propagation 

3 

C. l  Simply  restrained  specimens 

C.1.1  Gas  metal  Arc  Welding 


(a) 

One 

30" 

X 

30” 

X 

1/2 

"  SAE  1020  specimen 

(b) 

Two 

30" 

X 

30" 

X 

1" 

SAE  1020  specimens 

(O 

Two 

30" 

X 

30" 

X 

2" 

SAE  1020  specimens 

(d) 

Two 

30" 

X 

30" 

X 

7/8 

"  HY-130  specimens 

(e) 

One 

30" 

X 

30" 

X 

2" 

HY-130  specimen 

Details  of  the  experimental  procedures  used  and  recorded  results  are 
provided  in  the  three  technical  progress  reports  to  the  Office  of  Naval 
Research . 

2.1.3  Heat  Transfer  Analysis 

Efforts  to  improve  the  prediction  capabilities,  of  both  analytical 
and  numerical  solutions  of  the  heat  transfer  problem  during  welding,  were 
successfully  completed  during  the  first  half  of  the  fourth  year  of  this 
research  project,  and  are  summarized  in  this  section.  Details  on  this  part 

^ These  measurements  were  performed  using  the  Rosenthal-Norton  sectioning 
technique  (see  Section  2.6  of  First  Technical  Progress  Report  to  O.N.R., 
November  30,  1979  for  a  detailed  description  of  the  method). 

^ Same  specimens  as  in  A. 3.1. 
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of  the  research  program  can  be  found  in  the  third  technical  progress  report 
to  O.N.R.  and  the  following  thesis: 

Papazoglou,  V.J.,  "Analytical  Techniques  for  Determining  Temperatures, 
Thermal  Strains,  and  Residual  Stresses  During  Welding",  Ph.D.  Thesis, 

M.I.T. ,  May  1981. 

Closed-form  analytical  solutions  of  the  problem  of  heat  transfer 
during  welding,  though  not  as  accurate  as  the  numerical  ones,  are  very 
useful.  They  provide  for  the  establishment  of  the  general  laws  that  govern 
the  phenomenon,  thus  facilitating  a  better  understanding  of  it.  Moreover, 
they  allow  for  a  series  of  parametric  analyses  that  help  one  understand 
the  relative  importance  of  the  various  parameters. 

Based  on  these  considerations,  efforts  were  undertaken  during  the 
course  of  this  investigation  towards  devising  new  or  improving  already 
existing  analytical  solutions,  with  the  aim  of  obtaining  more  accurate 
predictions . 

Specifically,  a  new  three-dimensional  finite  heat  source  model  for 
solving  the  governing  partial  differential  equation  of  heat  transfer  was 
initially  developed,  taking  into  account  the  temperature  dependence  of  the 
material  properties  and  the  surface  heat  losses.  A  parametric  investiga¬ 
tion  of  the  effects  of  the  variation  of  various  model  parameters  was 
then  performed.  Results  were  finally  compared  with  other  closed  form 
solutions  and  with  experiments.  Agreement  was  generally  good  for  the 
case  of  thin  plates. 

In  addition  to  the  above  work,  the  conventional  point  heat  source 
4 

closed  form  solution  was  modified  so  as  to  more  accurately  simulate 
multipass  welding.  This  modification  permitted  positioning  of  the 
point  heat  source  at  the  center  of  each  welding  pass  (at  any  point  through 
the  plate's  thickness).  Comparison  of  the  results  obtained  was  again 
performed  with  experiments  and  the  previously  mentioned  finite  heat 
source  model. 

__ 

Rosenthal,  D.,  "Mathematical  Theory  of  Heat  Distribution  During  Welding 
and  Cutting",  Welding  Journal,  20  (5),  1941,  pp.  220s-234s. 
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Numerical  solutions  of  the  heat  transfer  problem  during  welding 
were  also  developed  in  Task  1.  Specifically,  using  the  multipurpose 
finite  element  program  ADINAT,  a  nonlinear  finite  element  model  was 
devised.  The  temperature  dependence  of  material  properties  and  the 
heat  losses  from  the  surface  were  again  taken  into  account.  The 
three-dimensional  character  of  the  heat  transfer  problem  was  simplified 
by  analyzing  only  a  cross-section  of  the  weldments,  located  in  the  mid¬ 
length  of  the  weld.  Arc  initiation  and  termination  effects  were 
neglected  and  the  welding  heat  source  was  assumed  to  move  at  a  constant 
high  speed  along  a  straight  line  in  a  planar  weld.  It  should  also  be 
noted  here  that  the  exact  heat  transfer  mechanism  inside  the  molten 
metal  pool  was  not  modeled.  Instead,  the  thermal  conductivity  of  the 
molten  metal  was  assumed  to  be  an  order  of  magnitude  higher  than  that 
of  the  material  in  the  solidus  temperature. 

Results  of  the  finite  element  modeling  are  presented  and  compared 
with  experiments  in  the  third  technical  progress  report  and  in  Papazoglou's 
thesis . 

2.1.4  Thermal  and  Residual  Stress  Analysis 

Using  the  predicted  temperature  distributions,  one  can  calculate  the 
transient  strains,  transient  stresses  and  residual  stresses  due  to 
welding.  The  calculation  of  strains  and  stresses,  however,  poses  a  much 
more  formidable  problem  than  the  one  encountered  in  the  heat  transfer 
analysis,  making  the  use  of  numerical  techniques  a  necessity.  These 
difficulties  stem  from  the  complicated  thermal-elastic-plastic  state 
developed  in  and  around  the  weld  metal  during  welding. 

Two  general  techniques  were  followed  by  the  M.I.T.  investigators  to 
solve  the  problem.  The  first  one  is  based  on  the  rather  simple  one- 

"Bathe,  K.J.,  "ADINAT-A  Finite  Element  Program  for  Automatic  Dynamic 
Incremental  Nonlinear  Analysis  of  Temperature,"  AVL  Report  82448-5, 

Mech.  Eng.  Dept.,  M.I.T. ,  May  1977  (revised  1981). 
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dimensional  analysis .6  Results  using  this  analysis  were  reported  in  the 
first  two  technical  progress  reports  to  O.N.R.  Some  details  of  this 
approach  are  also  given  in  Section  3.3.4  of  this  report,  where  the 
method  is  extended  so  as  to  calculate  strains  and  stresses  not  only 
during  welding  but  also  duiing  thermal  stress  relieving.  In  general, 
it  could  be  noted  that  the  one-dimensional  analysis  is  best  suited  for  the 
case  of  relatively  long  and  thin  (less  than  1/4"  thick)  mild  steel  plates. 

The  second  approach  is  based  on  the  finite  element  method,  and  can 
handle  thick  high  strength  steel  sections.  A  model  for  taking  into 
account  the  phase  transformation  strains  during  welding  of  quenched  and 
tempered  steels  was  developed  and  incorporated  into  the  finite  element 
program  ADINA.'7 

Specifically,  using  the  CCT  diagram  of  the  material  (either  known 
experimentally  or  derived  from  the  TTT  diagram)  and  the  temperature 
distribution  during  and  after  welding,  the  history  of  microstructure 
formation  during  the  cooling  states  is  predicted.  The  combined  trans¬ 
formation  and  thermal  strains  are  then  estimated  and  incorporated  into 
the  modified  finite  element  program  ADINA,  replacing  the  conventional 
thermal  strains. 

Comparison  of  the  finite  element  predictions  with  experimental  data 
yielded  a  good  correlation  in  general.  Further  details  on  this  part  of 
the  research  program  are  covered  in  the  third  technical  progress  report 
to  O.N.R.  and  the  following  thesis: 

Papazoglou,  V.J.,  "Analytical  Techniques  for  Determining  Temperatures, 
Thermal  Strains,  and  Residual  Stresses  During  Welding",  Ph.D.  thesis, 

M.I.T. ,  May  1981. 

An  overview  is  also  given  in  both  Appendices  A  and  B  of  this  report. 

^Papazoglou,  V.J.,  "Computer  Programs  for  the  One-dimensional  Analysis  of 
Thermal  Stresses  and  Thermal  Movement  during  Welding",  Manual  //  2  of 
Report  under  Contract  N00014-75-C-0469  to  0,N.R.,  M.I.T.,  1977. 

^  Bathe,  K..J.,  "ADINA-A  Finite  Element  Program  for  Automatic  Dynamic 
Incremental  Nonlinear  Analysis",  AVL  Report  82448-1,  Mech.  Eng.  Dept., 
M.I.T.,  September  1975  (revised  December  1978). 
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2.1.5  Fracture  Analysis 

Efforts  to  analyze  the  fracture  experiments  performed  in  Task  1.10 
were  completed  successfully  during  the  first  half  of  the  fourth  year  of 
this  research  program.  These  efforts  consisted  of  numerically  calculating 
the  stress  intensity  factor  and  the  crack  opening  displacement  using  the 
finite  element  method  and  based  on  a  newly  developed  hybrid  crack  element. 

Specifically,  a  circular  blade  saw  cut  was  used  in  the  experiments 
to  simulate  a  propagating  crack  in  the  weld  metal  of  the  H-type  slit 
specimens.  Crack  opening  displacement  (C.O.D.)  and  surface  strains 
were  measured  and  used  by  an  iterative  finite  element  code  in  order  to 
estimate  the  transverse  residual  stress  distribution  through  the  thickness, 
and  the  elastic  strain  energy  released  during  each  notch  size  increment. 

After  the  residual  stress  distribution  (through  the  thickness)  was 
obtained,  the  stress  intensity  factor  was  calculated  by  another  finite 
element  code.  This  last  program  automatically  considers  the  increasing 
size  of  the  crack.  A  special  hybrid  crack  element  was  developed,  in 
order  to  account  for  the  crack  tip  singularity  and  the  stress  distribution 
applied  on  the  crack  surface.  In  addition  to  the  stress  intensity  factois, 
the  crack  opening  displacement  was  also  calculated  by  the  program  and  was 
then  compared  to  the  experimentally  measured  one.  Good  agreement  between 
the  two  was  observed  for  the  case  of  small  cracks.  Both  finite  element 
programs  that  were  developed  for  the  fracture  analysis  were  based  on 

g 

the  general  finite  element  program  FEABL. 

Further  details  on  this  part  of  the  research  program  can  be  found 
in  the  third  technical  progress  report  and  in  the  following  thesis: 

Goncalves,  E. ,  "Fracture  Analysis  of  Welded  Structures",  Ph.I). 
thesis,  M.I.T.,  May  1981. 

A  review  is  also  provided  in  Appendix  B  of  this  report. 

8  _ 
Orringer,  0.,  French,  S.E.,  and  Winreich,  M. ,  "User's  Guide  for  the  Finite 

Element  Analysis  (FEABL  2,  4  and  5)  and  the  Element  Generation  Library 

(EGL)",  Aeroelastic  and  Structure  Research  Lab.,  Department  of  Aero,  and 

Astro.,  M.I.T.,  January  1978. 
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2.2  Work  Accomplished  on  Task  2  -  Cylindrical  Shells 

2.2.1  General  Status 

The  objective  of  Task  2  of  this  research  program  is  to  experimentally 
and  analytically  study  temperatures,  thermal  strains,  residual  stresses 
and  distortion  in  girth  welded  high  strength  steel  (specifically  HY-130) 
unstiffened  cylindrical  shells.  All  steps  of  this  task  (see  Table  1.1) 
have  been  completed  as  originally  planned. 

The  work  on  this  tack  has  been  completed  during  the  first  three 
years  of  the  research  program.  Details  of  this  work  can  be  found  in  the 
first  two  technical  progress  reports  to  O.N.R.  In  what  follows,  a  brief 
summary  of  the  work  performed  is  outlined. 

2.2.2  Experimental  Program 

An  initial  experiment  was  performed  outside  M.I.T.  on  butt  welding 
two  unstiffened  low  carbon  steel  cylindrical  shells.  This  exploratory 
experiment  was  of  great  value,  providing  the  M.I.T.  investigators  with 
experience  that  was  later  utilized  in  a  similar  experiment  conducted 
at  M.I.T.  using  HY-130  steel  cylinders. 

The  second  experiment,  performed  at  the  M.I.T.  facilities,  involved 
HY-130  cylindrical  shells.  The  shells  were  formed  from  60  in. x  63  in.x 
it  in. HY-130  flat  plates  by  first  bending  them  to  semi-circular  cylin¬ 
drical  shells  and  subsequently  seam-welding  the  two  halves  together 
using  the  GMA  welding  process.  After  completion  of  the  seam  weld,  the 
>>0  in. long,  18  in. 0D  cylinder  was  sawed  into  three  equal  sections,  and 
; 'repnrnt ions  of  the  grooves  were  made  for  the  joining  of  two  of  these 
-lions.  finally,  the  two  cylindrical  sections  were  girth  welded 
•  'ether  using  the  GMA  welding  process. 

i empe t a t uro  changes  and  transient  strains  in  the  axial  and  circum- 
N  teriti  il  directions  were  measured  at  preselected  locations  during  both 
' simp;  mt!  girth  welding  operations,  using  thermocouples  and  electric 
;  •  -s  i  ••  t  nice  st  r  tin  gages  ,  respectively.  The  radial  distortion  at  two 
}  •  i  n  t  was  also  i  out  i  mums  l  v  measured  during  the  girth  welding  of  the 
'  w  •.  1  i  nilt  r  s  . 
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Finally,  using  electric  resistance  strain  gages  and  the  strain 
relaxation  technique,  the  residual  stresses  on  the  outer  surfaces  of 
the  welded  cylinders  were  measured.  The  axial  and  tangential  stresses, 
as  distributed  in  both  the  axial  and  circumferential  directions,  were 
determined  and  compared  with  data  from  other  investigators. 

2.2.3  Analytical  Program 

The  problem  of  welding  two  unstiffened  cylindrical  shells  was 
analyzed  using  existing  computer  programs  developed  at  M.I.T.  during  the 
past  several  years.  A  program  based  on  the  analytical  solution  for  a 
moving  point  heat  source  was  used  for  the  heat  flow  analysis.  Then, 
strains,  stresses,  and  distortion  during  welding  were  analyzed  using  an 
axisymmet ric  finite  element  program. 

Good  correlation  was  observed  between  the  numerically  and  experi- 
mencally  obtained  temperature  distributions.  However,  this  was  not  the 
case  for  the  strain  and  stress  analysis.  It  is  believed  that  if 
adequate  funding  was  available,  better  results  would  have  been  obtained 
using  the  modified  finite  element  programs  ADINAT  and  ADINA. 

2 . 3  Work  Accomplished  on  Task  3  -  Bronze  Weldments 
2.3.1  General  Status 

The  objective  of  Task  3  of  the  current  research  project  is  to 
experimentally  generate  information  on  temperature  changes,  transient 
strains  and  distortion  on  bronze  weldments.  In  particular,  Nickel- 
Aluminum  (Ni-Al)  bronze  was  chosen  as  the  material  to  be  investigated, 
because  of  the  problems  encountered  in  controlling  distortion  during 
the  manufacturing  of  propeller  blades  using  this  material.  It  should 
be  pointed  out  that  this  investigation  was  prompted  from  the  fact  that, 
to  the  best  of  our  knowledge,  no  published  information  is  available 
on  the  subject. 

The  work  on  this  Task  was  initiated  and  completed  during  the  fourth 
year  of  this  research  program.  Details  on  this  work  can  be  found  in  the 
third  progress  report  to  O.N.R.  and  the  following  thesis: 
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McCord,  R.S.,  "An  Investigation  of  Strain,  Distortion,  and  Heat 
Flow  Distribution  During  Welding  of  Nickel-Aluminum  Bronze",  Ocean 
Engineer's  thesis,  M.I.T.,  June  1981. 

In  what  follows,  a  brief  summary  is  given  of  the  work  performed. 

2.3.2  Literature  Review  and  Experimental  Program 

A  literature  survey  was  first  conducted  to  obtain  information  on  the 
mechanical  and  physical  properties  of  Ni-Al  bronze  and  their  dependence 
on  temperature.  This  was  necessary  for  the  numerical  analysis  of  heat 
flow  and  thermal  stresses  during  welding. 

Two  series  of  experiments,  using  a  total  of  eight  specimens,  were 
conducted  to  measure  temperature,  thermal  strain,  and  distortion  changes 
during  the  welding  of  Ni-Al  bronze  plates.  Emphasis  was  placed  on  compar¬ 
ing  the  final  distortions  of  the  plates  after  welding. 

The  first  series  of  experiments  consisted  of  laying  a  bead-on-plate 
weld  in  two  passes.  Specimens  1  through  6,  all  measuring  24  in.  x  12  in.  x 
1/4  in.,  were  welded  in  this  manner,  with  the  heat  input  being  held 
approximately  constant  for  all  specimens.  Specimens  #1  and  # 2  were  welded 
in  the  unrestrained  condition  so  that  baseline  data  could  be  gathered  for 
comparison  purposes  with  other  conditions.  Specimen  #3  was  restrained 
during  both  weld  passes  using  C-clamps.  Specimen  // 4  was  block  welded 
in  four  blocks;  in  each  block,  measuring  6  in  in  length,  both  passes  were 
completed  before  moving  on  to  the  next  one.  Specimen  //5  was  welded  in  the 
unrestrained  condition;  it  was  preheated  and  postheated  using  oxyacetylene 
torches,  and  finally  the  weld  was  ground  down  at  preselected  locations. 
Specimen  # 6  was  welded  in  the  restrained  condition. 

The  second  series  of  experiments  consisted  of  laying,  in  two  passes, 
a  bead  on  the  edge  of  two  specimens,  H7  and  #8,  each  measuring  21  3/4  in.  x 
5  3/4  in.x  1/4  in.  Heat  input  as  well  as  all  other  parameters  were  kept 
constant  in  both  cases. 

Temperature,  strain,  and  distortion  changes  were  recorded  continuously 
during  the  welding  of  all  specimens.  Results  are  reported  in  the  third 
progress  report  to  O.N.R.  and  in  McCord's  thesis. 
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3.  WORK  ACCOMPLISHED  ON  TASK  4  -  THERMAL  STRESS 

RELIEVING 


3.1  General  Status 

The  objective  of  Task  4  of  this  research  project  is  to  study  the 
mechanisms  and  evaluate  the  effectiveness  of  thermal  stress  relieving 
in  weldments.  Task  4  consists  of  the  following  subtasks  (see  Table  1.2): 

Task  4.1:  Initial  Analysis  and  Verification  Tests 

Task  4.2:  Further  Analysis  and  Experiments 

Both  subtasks  are  now  completed  as  originally  proposed. 

An  extensive  literature  survey  was  initially  performed  covering 
thermal,  mechanical  and  vibratory  stress  relieving  treatments.  Special 
attention  was  given  to  the  applicable  military  and  industrial  codes. 

An  initial  analysis  of  the  effectiveness  of  mechanical  and  thermal 
stress  relieving  treatments  was  performed  for  the  cases  of  long,  thin 
plates  butt-  or  edge-welded.  The  one-dimensional  computer  program  earlier 
developed  at  M.I.T.  was  modified  to  calculate  strains  and  stresses  during 
and  after  uniform- temperature  heat  treatments.  Stress  relaxation,  due  to 
creep  at  high  temperatures,  was  also  included  in  the  analytical  model. 

Experiments  to  verify  the  analytical  results  were  then  conducted 
using  304  stainless  steel  plates.  This  was  decided  since  extensive 
creep  data  for  HY-130  steel  could  not  be  located  in  the  literature. 

Further  analysis  was  then  performed  in  order  to  investigate  the 
effectiveness  of  low-temperature  thermal  stress  relieving  treatments. 
Specifically,  the  one-dimensional  computer  program  was  further  modified, 
so  as  to  calculate  the  temperature,  strain,  and  stress  distributions 
during  and  after  line  heating  with  a  flame. 

Finally,  verification  experiments  were  conducted  using  edge-welded 
HY-130  steel  plates.  Both  uniform  and  line  heating  were  tested. 

A  summary  of  the  work  performed  is  presented  in  the  next  sections 
of  this  chapter.  Additional  data  is  in  the  following  theses: 
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Agapakis,  J.E.,  "Welding  of  High  Strength  and  Stainless  Steels: 

A  Study  on  Weld  Metal  Strength  and  Stress  Relieving",  S.M.  Thesis, 

M.I.T. ,  May  1982. 

Carpentier,  K.P.,  "Thermal  Stress  Relief  of  HY-130  Weldments", 

S.M.  thesis,  M.I.T. ,  May  1982. 

3. 2  Literature  Survey 

It  is  a  well  established  fact  that  during  the  fabrication  of  welded 
structures,  residual  stresses  are  developed  having  peak  values  that  reach 
the  yield  stress  of  the  material  being  welded.  This  fact  was  confirmed 
by  both  analysis  and  experiments  for  the  case  of  HY-130  steels  in  the 
course  of  the  present  research  program. 

Residual  stresses  must  be  a  major  cause  of  concern  since  they  usually 
are  detrimental,  directly  or  indirectly,  to  the  integrity  and  the  service 
behavior  of  a  welded  structure.  As  a  consequence,  several  methods,  aimed 
at  reducing  residual  stresses,  have  been  proposed  and  are  currently  in 
use.  These  methods  include  both  mechanical  and  thermal  treatments. 

In  the  next  few  parts  of  this  section,  a  brief  evaluation  of  the 
various  methods  is  attempted  and  some  important  results  are  presented. 

A  more  complete  review  is  given  in  the  previously  mentioned  theses. 

3.2.1  Thermal  Stress  Relieving  Treatments 


Thermal  stress  relieving  basically  involves  heating  of  the  part  to 

9 

a  subcritical  temperature,  below  A  ,  holding  it  at  that  temperature  to 

C1 

ensure  uniformity,  and  slow  cooling  to  room  temperature,  usually  in  air, 
to  prevent  the  reintroduction  of  stresses.  The  stress  relieving  temper¬ 
atures  usually  are  of  the  order  of  1100°F  to  1300°F  (590°C  to  700°C) , 
where  the  yield  strength  has  drastically  decreased  and  creep  occurs. 

The  welding  residual  stresses  can  no  longer  be  supported.  Thus  the 
stress  distribution  will  be  relatively  uniform  and  at  a  ver>  low  level. 


’By  A]  we  denote  the  lower  critical  temperature,  where  perlite  begins  to 
transform  into  austenite  (A  is  used  for  heating  and  A  for  cooling). 

C1  rl 
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Up  to  ,  the  higher  the  stress-relieving  temperature,  the  more  completely 
the  stress  is  removed. 

Full  relief  of  the  residual  stresses  can  only  be  ensured  through  an 
annealing  treatment.  This,  however,  would  be  more  costly  and  time 
consuming  and  would  cause  more  problems  due  to  the  higher  temperatures. 

Stress  relieving  is  always  followed  by  some  dimensional  changes. 

Even  warpage  and  distortion  can  result  when  the  residual  stresses  are 

high  enough.  Therefore,  it  may  be  necessary  to  straighten  the  part 

and  sometimes  to  stress  relieve  again  to  reduce  the  straightening  stresses. 

Stress  relief  heat  treatment  can  be  applied  either  to  the  structure 
as  a  whole,  or  to  specific  portions  of  it  if  the  size  is  prohibitively  large. 
In  the  former  case,  electric  and  fossil-fired  furnaces  are  generally  used, 
whereas  the  localized  heat  treatment  is  effected  by  means  of  electric 
resistance,  induction,  gas  flame,  and  exothermic  heating  techniques. 

For  loc.a.1.  stress  relieving  heat  treatments,  it  is  absolutely  necessary 
to  ensure  that  the  temperature  distribution  during  the  heat  treatment 
does  not  induce  new  thermal  stresses  which  exceed  the  material  yield 
stress  and  can  lead  to  the  development  of  new  residual  stresses  on  cooling. 
This  imposes  strict  requirements  on  the  level  and  uniformity  of  temperature, 
on  the  heating  and  cooling  rates,  and  on  the  width  of  the  heating  zone. 

The  latter  largely  determines  the  existing  temperature  gradient  through 
the  thickness  (between  the  heated  and  unheated  surfaces). 

For  butt  welded  plates,  it  was  experimentally  proven  that  satisfactory 
relief  of  residual  stresses  can  be  expected  if  uniform  heat  input  is 
applied  over  a  band-width  of  twice  the  length  of  the  weld,  as  shown  in 
Figure  3.1a.^  According  to  British  codes,  this  bandwidth  becomes  equal 
to  5 kRt  for  a  circumferential  girth  weld  on  a  cylinder  of  radius  R  and 
thickness  t  (Figure  3.1b).  Different  specifications,  however,  are  spelled 
out  in  different  welding  codes. 

^Burdekin,  F.M.,  "Heat  Treatment  of  Welded  Structures",  The  Welding 
Institute,  Abington,  England,  2nd  edition,  1969. 
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FIGURE  3.1  Band  widths  for  local  heat  treatments  to  achieve  stress 

relief  in:  (a)  flat  plate  situations,  (b)  circumferential 
butt  joints  in  cylinders,  and  (c)  butt  joints  in  spheres 


Nevertheless,  there  are  generally  three  main  requirements  that 
govern  the  selection  of  the  heat  treating  method: 

(a)  It  should  be  able  to  produce  the  required  temperature. 

(b)  The  required  temperature  should  be  controlled  within  specified 
limits  (e.g.,  +  10°  to  20°F  for  steels). 

(c)  It  should  be  possible  to  achieve  a  uniform  and  even  heating  and 
cooling  rate  throughout  the  heaviest  section  to  be  treated. 
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The  last  requirement  is  especially  important  for  the  case  of  joints  of 
complex  geometry  and  variable  thickness. 

3.2.2  Effects  of  Thermal  Stress  Relieving 

Stress  relieving  heat  treatments  are  usually  effective  in  reducing 
the  high  residual  stresses  present  in  a  weldment.  However,  they  also  have 
an  effect  on  the  microstructure  and  the  properties  of  both  the  base  plate 
and  the  weld  metal  since  they  are  carried  out  at  relatively  high  temper¬ 
atures.  These  effects  vary  with  the  material  under  consideration. 

During  the  last  decade,  the  desirability  of  post  weld  heat  treatments 
and  their  effects  were  extensively  investigated  by  the  Working  Group  on 
Thermal  Stress  Relief  of  the  Commission  X  of  the  International  Institute 
-if  Welding.  Information  from  the  series  of  documents^that  were  produced 
is  summarized  in  this  section. 

Specifically  for  non  work-hardened  base  metal  of  C-Mn  and  micro- 
all  eyed  steels,  it  was  generally  concluded  that  tensile  properties  are 
impaired  to  a  significant  extent,  especially  at  higher  temperatures. 
Resistance  to  brittle  fracture  is  affected  but  not  drastically.  Tem¬ 
perature  is  a  more  important  factor  than  soaking  time.  For  low  alloy  and 
creep  resistant  steels,  which  are  used  in  a  normalized  and  tempered  or 
a  quenched  and  tempered  condition,  the  effect  of  stress  relieving 
treatments  on  the  properties  will  depend  on  the  temperature.  The  effect 
will  be  minimal  if  the  treatment  is  carried  out  at  a  temperature  lower 
than  that  of  initial  tempering.  Additional  tempering  will  result  however, 
if  higher  temperature  treatment  is  performed.  This  is  usually  beneficial 
for  the  resistance  to  brittle  fracture,  but  detrimental  for  other 
properties  such  as  creep  resistance,  and  usually  is  not  recommended  by 
the  codes.  On  the  other  hand,  for  the  case  of  work-hardened  base  material, 
the  heat  treatments  usually  restore  the  base  properties  and  prevent  strain 
aging,  and  are  therefore  beneficial. 

1.1. W.  Working  Group  on  Stress  Relief  by  Heat  Treatment,  "Final  Report: 
Desirability  of  Post  Weld  Heat  Treatments  in  Welded  Construction,"  I.l.W. 
Document  X-913-78,  February  1979. 
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The  effect  of  stress  relieving  heat  treatments  on  the  properties  of 
the  heat  affected  zone  (HAZ)  will  not  only  depend  on  the  type  of  steel 
but  also  on  the  microstructural  state  of  the  HAZ.  Therefore,  welding 
procedure  and  conditions,  heat  input,  plate  thickness  and  distance  to  the 
fusion  line  are  important  parameters.  For  C-Mn  steels  the  heat  treatment 
will  in  general  soften  the  HAZ  structures,  except  if  carbide  precipi¬ 
tation  occurs.  The  yield  strength  of  the  HAZ  will  usually  be  higher 
than  that  of  the  base  material  after  the  same  treatment.  In  most  cases, 
heat  treatments  are  also  beneficial  for  resistance  to  brittle  fracture 
for  these  steels.  For  low  alloy  steels,  the  main  problem  is  to  retain 
satisfactory  toughness  in  the  HAZ.  The  effect  of  treatment  is  usually 
strongly  dependent  on  the  type  of  steel  and  the  exact  metallurgical 
changes  associated  with  the  welding  and  heat  treating  tempering. 

Although  not  much  work  has  been  done  on  the  effect  of  heat  treat¬ 
ments  on  the  properties  of  weld  metal,  it  appears  that  tensile  properties 
diminish  considerably  on  tempering.  Again,  temperature  seems  to  be  more 
important  than  the  soaking  time.  As  in  the  base  metal,  embrittlement 
can  occur  and  instances  of  stress  relief  cracking  have  been  reported. 

Stress-relief  cracking  is  defined  as  "intergranular  cracking  in  the 
heat  affected  zone  or  weld  metal  that  occurs  during  the  exposure  of 
welded  assemblies  to  elevated  temperatures  produced  by  post-weld  heat 
treatments  or  high-temperature  service".  It  has  also  been  referred 
to  in  the  literature  as  "post-v;eld-heat  cracking"  or  "reheat  cracking" 
and  in  general  is  caused  when  some  relief  of  stresses  by  creep  occurs. 

This  form  of  cracking  became  a  problem  with  the  austenitic  stainless 
steels  in  the  1950' s  and  with  the  low-alloy  constructional  steels  in  the 
1960's.  It  also  occurs  in  ferritic  creep-resisting  steels  and  nickel 
base  alloys  and  is  generally  related  to  precipitation  hardening.  Non¬ 
precipitation  hardening  materials,  such  as  plain  carbon  steels  and 
certain  nickel  alloys,  are  not  susceptible  to  reheat  cracking. 

The  cracks  can  be  positively  identified  by  metallographic  examination 
due  to  their  characteristic  branching  intergranular  morphology,  along 
the  coarse-grain  region  of  the  heat-affected  zone.  Cracking  usually 
occurs  at  high  temperature,  when  creep  ductility  is  insufficient  to 
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accomraodate  the  strains  required  for  the  relief  of  applied  or  residual 
stresses.  When  residual  stresses  are  high,  as  in  thicker  and  restrained 
sections,  reheat  cracking  is  most  likely  to  occur. 

Extensive  investigations  of  the  cracking  mechanism  and  of  possible 

remedies  have  been  undertaken  all  over  the  world,  and  are  in  detail 

12  12 

reviewed  by  Meitzner  and  Dhooge,  et  al .  In  1970  1.1. W.  established, 
in  Commission  X,  a  Working  Group  on  "Reheat  Cracking"  to  collect  and 
assimilate  information  on  the  subject. 

2.2.3  Stress  Relieving  of  HY-130  Steels 

The  very  high  residual  stresses  that  develop  during  welding  of  HY-130 
steels  make  stress  relief  treatments  very  attractive.  Substantial 
reduction  of  stresses  was  shown  to  occur  in  both  the  base  and  weld  metal 
at  temperatures  between  950°F  (510°C)  and  1050°F  (566°C) . ^  However,  a 
severe  degradation  of  notch  toughness  may  also  occur  at  these  temperatures. 
This  embrittlement,  which  is  phenomenologically  similar  to  temper 
embrittlement,  is  believed  to  be  influenced  by  the  soaking  temperature,  the 
time  at  temperature,  the  cooling  rates,  the  plastic  strain,  the  exact 
alloy  composition  and  prior  heat  treatment. 

Extensive  investigation  on  the  subject  was  performed  by  the  U.S.  Steel 
Corporation  and  the  U.S.  Navy  (N.S.R.D.C.  Materials  Laboratory).  Results 
by  Rosenstein^ indicated  that  stress  relief  cycles  are  cumulative  in  nature, 

^“Meitzner,  C.F.,  "Stress  Relief  Cracking  in  Steel  Weldments",  An  Inter¬ 
pretive  Report,  Welding  Research  Council,  Bulletin  211. 

^Dhooge,  A.,  et  al.,"A  Review  of  Work  Related  to  Reheat  Cracking  in  Nuclear 
Reactor  Pressure  Vessel  Steels",  International  Journal  of  Pressure  Vessels 
and  Pipeing  No.  (6)  1978,  pp.  329-409. 

14 

N.S.R.D.C.  reports  on  "Stress  Relief  Embrittlement  of  AX-140  and  E-11018 
Weld  Metals",  and  "Stress  Relief  Characteristics  of  a  5%  Ni  Weld  Metal", 
obtained  after  private  communication,  October  1981. 

Rosenstein,  A.H.,  "Phenomenological  Investigations  of  Stress  Relief 
Embrittlement,  Welding  Journal,  March  1970,  pp.  122s-131s. 
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by  comparing  the  tensile  and  Charpy  properties  resulting  from  consequent 
heat  treatments  of  short  duration  and  one  single  treatment  of  long  duration, 
both  having  the  same  soaking  times. 

Furthermore,  it  was  observed  that  the  degradation  of  notch  toughness 
was  maximum  when  stress  relieving  at  950°F  (510°C)  and  minimum  when  at 
1050°F  (566°C).  However,  investigations  of  the  effects  of  stress  relieving 
on  the  weld  metal  have  shown  that  the  higher  temperatures  may  substantially 
reduce  the  yield  and  tensile  strength  of  the  weld  metal,  and  thus  may 
cause  undermatching.  Therefore,  the  most  practical  temperature  range  for 
effective  stress  relieving  of  HY-130  weldments  is  between  1000°F  (538°C) 
and  1025°F  (552°C). 

Additionally,  it  was  established  that  the  degradation  of  toughness 
at  950°F  (510°C)  occurs  during  both  soaking  and  cooling.  The  isothermal 
degradation  was  shown  to  be  directly  dependent  upon  time  at  temperature. 

The  degradation  during  cooling,  on  the  other  hand,  is  inversely  related 
to  the  cooling  rate  and  does  not  depend  on  the  soaking  time.  Therefore, 
cooling  embrittlement  constitutes  the  major  portion  of  degradation  after 
short  times  at  temperature  and  only  a  minor  one  after  long  soaking 
periods . 

It  was  also  shown  that  stress  relief  at  the  tempering  temperature  of 
1050°F  (566°C)  results  in  softening  at  temperature  (accompanied  by 
increased  toughness)  and  embrittlement  on  cooling.  Therefore,  the 
resulting  properties  will  depend  both  on  the  cooling  rate  and  the  soaking 
time.  Furthermore,  it  should  be  mentioned  that  the  toughness  degradation 
due  to  stress-relief  can  in  general  be  recovered  by  retempering  to  a 
lower  strength  level. 

3.2. 4  A1 ternat i vo  Methods  of  Stress  Relieving 

When  the  metallurgical  effects  of  thermal  treatments  are  undesirable, 
it  is  possible  to  achieve  stress  relief  by  mechanical  means. 

One  such  method  is  the  well-known  overstressing  above  yield.  In 
general,  loading  the  welded  structure  to  a  level  higher  than  its  intended 
service  loading  will  reduce  the  effect  of  the  existing  residual  stresses 
on  its  service  performance.  Additionally,  a  first  overloading  succeeds 
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in  producing  compressive  residual  stresses  around  any  undetected  defects, 
resulting  in  what  are  believed  to  be  beneficial  effects  on  the  fatigue 
characteristics  of  the  weldment.  The  mechanical  treatment,  however, 
does  not  produce  any  improvement  in  the  fracture  toughness  of  the  HAZ  as 
opposed  to  the  thermal  treatments. 

An  alternative  to  overstressing  is  the  vibratory  stress  relief  (VSR) 
treatment.  Experimental  work  on  this  method  started  in  the  1930 's  and 
initially  gave  rise  to  very  optimistic  predictions.  During  the  last 
decade  both  experimental  and  analytical  work  was  intensified  in  an  effort 
to  investigate  the  mechanisms  and  effectiveness  of  the  VSR.  Although  this 
technique  is  currently  used  in  the  industry  with  some  success,  there 
still  remains  much  to  be  learned  about  its  mechanism,  or  even  if  it 
actually  does  relieve  stresses  as  claimed.  Nevertheless,  the  energy 
savings  accrued,  in  comparison  to  thermal  stress  relief  treatments,  is  a 
sufficient  enough  reason  to  exploit  its  use. 

The  impact  from  explosive  contact  charges  was  also  shown  to  be 
capable  t  redistribute  (rather  than  to  relieve)  welding  residual  stresses. 
Important  parameters  in  such  a  treatment  are  the  intensity  and  distribution 
of  the  explosive  load. 

For  the  purposes  of  this  study,  however,  no  further  consideration  was 
given  to  mechanical  or  explosive  stress  relieving  techniques.  Instead, 
the  method  of  thermal  stress  relief  treatment  was  examined  both  analytical¬ 
ly  and  experimentally. 

The  metallurgical  aspects  of  thermal  stress  relieving  w°re  also 
considered  to  be  outside  the  scope  of  this  investigation.  In  what  follows 
only  the  effectiveness  of  thermal  stress  relieving  treatments  in  terms 
of  the  attained  reduction  of  residual  stresses,  is  examined. 

3 .  1  Ana  lys is  of  Residua  1  St ress  Relaxation  Due  to  Heat  Treatments 

Residual  stress  changes  can  arise  during  all  three  stages  of  a  heat 
treatment,  Specifically,  referring  to  Figure  3.2  during  the  heating 
part  of  the  process,  residual  stresses  decrease  clue  to  the  temperature 
dependence  of  the  mechanical  properties,  mainly  through  a  reduction  o!  the 
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yield  strength  with  temperature.  During  the  holding  (or  soaking)  period, 

the  temperature  is  kept  constant  and  residual  stresses  are  reduced  due 

to  creep.  Experimental  evidence  suggests  that  the  major  portion  of  this 
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stress  reduction  occurs  in  the  first  part  of  this  stage.  Finally, 
during  the  cool-down  period,  residual  stresses  increase  due  to  the  tem¬ 
perature  dependence  of  the  mechanical  properties,  but  hopefully  (when 
the  treatment  is  successful)  not  to  their  initial  levels. 
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FIGURE  3.2  Stress  relieving  temperature  history 


IbX.S.R.D.U.  reports,  op.cit. 
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3.3.1  General  Considerations 

To  judge  the  effectiveness  of  a  stress  relief  treatment,  with  regards 
to  the  accomplished  reduction  of  residual  stresses,  it  should  be  necessary 
to  measure  the  maximum  residual  stresses  before  and  after  the  treatment. 

The  difference  would  be  a  realistic  measure  of  performance.  However, 
an  acceptable  alternative  to  the  time-consuming,  costly  and  usually 
destructive  residual  stress  measurements,  would  be  a  proper  analytical 
model.  Furthermore,  such  a  model  would  be  very  helpful  in  determining  an 
optimal  lower  temperature  heat  treatment,  where  a  properly  selected 
heating  pattern  would  most  effectively  reduce  stresses,  while  keeping 
the  metallurgical  changes  minimal. 

in  that  direction,  various  approaches  have  been  followed  in  the 
Literature  by  several  investigators.  Very  simple  uniform  residual  stress 
distributions  are  usually  assumed  for  the  weld  metal,  so  that  the  uni- 
dimensionil  stress-strain  curves  can  be  directly  employed.  Such  analytical 
results  arc  obtained  and  experimentally  verified  by  Tanaka.^  For  more 
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md  two-  or  three-dimensional  stress  states  numerical  models 


have  beer  proposed  to  handle  the  thermal-elastic-plastic  and  creep  analysis 
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required.  u'eda  and  Fukuda  present  a  finite  element  model  capable  of 

calculating  welding  residual  stresses  and  stress  relief  due  to  creep. 

I  9 

Fujita,  et  al  . ,  develop  a  thermo-visco-elastic-plastic  model  to  study 
tlu  median  ism  of  stress  relief  annealing.  Finally,  Cameron  and  Pember ton^ 


Tanaka,  .).,  "Decrease  in  Residual  Stresses,  Change  in  Mechanical  Properties 
and  CracKing  Due  to  Stress  Relieving  Heat  Treatment  of  HY-80  Steel", 

Welding  in  the  World,  J_0  (1/2),  1972. 

1 8 

I'eda,  H.,  and  Fukuda,  K.,  "Analysis  of  Welding  Stress  Relieving  by 
Annealing  Based  on  Finite  Element  Method",  Trans,  of  Japan  Weld.  Res. 

1 n st . ,  4  ( 1 ) ,  1975. 

1  9 

nijita,  Y . ,  Nomoto,  T.,  Aoyagi,  A.,  "A  Studv  on  Stress  Relaxation  due  to 
Heat  Treatment",  Dept,  of  Naval  Architecture,  University  of  Tokyo,  May 
1973,  IIW  Doe.  X-697-71. 

>  (  ) 

Cameron,  i.G.,  and  Pemberton,  C.S.,  "A  Theoretical  Sludv  o!  Thermal 
Stress  Relief  in  Thin  Shells  of  Revolution",  Inti.  J.  lor  Nnm.  Meth. 
in  Fng.,  II,  1977,  pn .  142 3- 11 37. 
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present  a  numerical  model  of  the  thermal  stress  relief  in  thin  shells  of 
revolution . 

For  the  purposes  of  the  current  study,  it  was  decided  that  the 

analysis  of  the  thermal  stress  relieving  operation  be  accomplished  using 

the  simple  one-dimensional  computer  program,  successfully  employed  in  the 

past  at  M.I.T.  for  the  prediction  of  residual  stresses  in  long,  thin,  butt 
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or  edge  welded  plates.  Some  modifications  had  to  be  made  in  this  program 
in  order  to  handle  the  case  of  thermal  stress  relieving.  These  modifi¬ 
cations  will  be  dealt  with  below  at  the  appropriate  places. 

3.3.2  Description  of  the  One-Dimensional  Model  -  Assumptions 

The  fundamental  assumptions  incorporated  in  this  model  are: 

(a)  The  plate  is  infinite  and  very  thin  (referring  to  Figure  3.3 
L  -*•  03  and  H  ->  0)  . 

(b)  The  welding  arc  is  modeled  as  a  line  heat  source,  and  there 

!.-  no  temperature  gradient  through  the  thickness  of  the  plate. 
(Two-dimensional  temperature  distribution). 

lc )  Furthermore,  the  temperature  distribution  is  stationary  if 
viewed  from  a  system  moving  with  the  heat  source.  (Quasi¬ 
stationary  state). 

(!)  Stress  is  non-zero  only  in  the  direction  parallel  to  the  weld 
centerline.  (One-dimensional  stress  distribution). 

(e)  These  stresses  are  a  function  of  the  transverse  distance  from 
the  weld  centerline  only. 

Additional  assumptions  for  the  analysis  of  thermal  stress  relief 
treatment  were  nude  as  follows: 

(f)  Anv  arbitrary  temperature  distribution  and  history  would  be 
input  to  the  modified  program.  However,  in  this  project  only  uniform  and 
line  heating  were  analyzed.  In  the  former  case,  a  uniform  temperature 
distribution  was  assumed  over  the  entire  plate,  changing  with  time  as  in 
Figure  3.2. 

1 

I’apaziu’ I  ou ,  V.  "Computer  Programs  for  the  One-Dimensional  Analysis  of 
Thermal  Sties.,-  -  and  Metal  Movement  during  Welding",  M.I.T.,  danuarv  1  9  7  7  . 
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In  the  case  of  line  heating,  however,  the  temperature  distribution 
during  the  treatment  was  calculated  in  a  way  analogous  to  the  temperature 
distribution  during  welding. 

Two  further  assumptions  were  introduced  fcr  the  case  of  uniform 
heating : 

(g)  Due  to  the  relatively  fast  heating  and  cooling  rates  it  was 
assumed  that  no  creep  occurs  during  the  heating  and  cooling 
stages  of  the  treatment. 

(h)  During  the  holding  or  soaking  period  at  the  stress  relieving 
temperature,  residual  stresses  can  only  decrease  due  to  creep. 

In  other  words,  if  creep  is  not  included  in  the  model,  no 
change  in  stresses  will  take  place  during  this  period. 

3.3.3  Temperature  Distribution 

Duri  welding,  the  non-uniform  and  changing  with  time  temperature 
distribution  is  estimated  in  the  one-dimensional  program  by  the  well-known 
Rosentha.L  solution.  The  exact  solution  for  a  line  heat  source  moving 
along  a  thin  and  infinite  plate  (as  in  Figure  3.3)  is: 


0 


= 

o  2-nkfl 


K  (~ ) 
o  2k 


(3.1) 


where  0  is  the  temperature  at  point  (x,  y)  at  time  t,  0^  is  the  initial 
temperature,  H  the  plate  thickness,  k  the  thermal  conductivity,  k  the 
thermal  diffusivity  (k  =  k/pc^),  p  is  the  density,  c^  the  specific  heat, 
Q  the  total  heat  input,  v  the  welding  speed  and  Ko(x)  is  the  modified 
Bessel  function  of  second  kind  and  zero  order. 

The  moving  coordinates  4  and  r  are: 


4  =  x  -  vt 


and 


r 


?  ? 
(■;“+y  ) 


1/2 


(3.2) 


(3.3) 


FIGURE  3.3  Weldment  configuration  (Butt  welding  of  plates) 


The  total  heat  input,  Q,  is 

Q  =  V.I.n  (3 

a 

where  V,  I  are  the  arc  voltage  and  current  respectively  and  n  the  arc 

3 

efficiency . 

Finally,  it  should  be  noted  that  Equation  (3.1)  was  also  modified 
to  account  for  heat  losses  due  to  ratiation  and  convection  from  the 
surfaces  of  the  plate,  becoming: 


where  h  =  average  surface  heat  loss  coefficient. 
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During  heat  treatment  in  a  furnace  the  temperature  distribution  can 
be  assumed  to  be  uniform  along  the  entire  plate.  For  the  case  of  a 
localized  treatment,  however,  by  flame  heating  for  example  the  temperature 
distribution  can  be  calculated  modifying  the  solution  for  a  point  heat 
source  moving  on  an  semi-infinite  body: 


(3.5) 


where 


(3.6) 


and  all  the  other  variables  same  as  in  Equation  (3.1). 

This  solution  of  the  governing  differential  equations  satisfies  the 
following  boundary  conditions: 

Um  2r  R2k  —  )•  =  Q  (3.7) 

R-*  O'  ' 

and  also 


?.  im  0  =  0  (3.8) 

R-^oo  0 


Condition  (3.7)  is  the  heat  equilibrium  equation  on  an  infinitely 
small  hemisphere  drawn  around  the  point  heat  source.  Condition  (3.8) 
merely  states  that  the  temperature  must  not  change  at  great  distances 
from  the  source. 

In  order  to  more  accurately  describe  the  temperature  distribution, 
however,  for  the  case  of  a  finite  width  and  finite  thickness  plate  heated 
with  a  traveling  heat  source,  as  in  Figure  3.4  two  more  boundary  conditions 
have  to  be  satisfied: 


±0. 

"!Z 


0  at  z  =  0  and  z=H 


(3.9) 


<  K 


TRAVELLING  HEAT  SOURCE 


EXISTING  EDGE  WELD 


SECTION  A-A 


FIGURE  3.4  Geometry  of  plate  undergoing  strip  heating 


— —  =Oaty=G  and  y  =  -F 

I  y 


(3.10) 


To  satisfy  these  boundary  conditions  an  infinite  number  of  heat  source 
images  with  respect  to  the  plate  bounaaries  must  be  used,  as  in  Figure 
3.5.  Then  the  solution  becomes: 
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I-  ,  dTL  .  dTR  ,  dBL  ,  nBR- 

^  -\vR  -AvR  -AvR  -\vR 

■:  y.  e _ nm  e _ nm  e _ nm  e  nm 

n= 1  m=l  _TL  _TR  _BL  BR 


(3.11) 


Rn  V  +  y  +  (2nli_z> 
Rn  =V?2  +  y2  +  (2nH+z)2 


A  finite  number  of  images  was  used  in  the  computer  implementation  of 
Equation  (3.11) . 

J . 3.4  Stress  Analysis 

To  calculate  the  transient  and  residual  stresses  during  and  after 

welding  and  subsequent  heat  treatments,  the  method  of  successive  elastic 

,  ,  22 

solutions  is  employed.  The  procedure,  outlined  by  Mendelson,  was  first 

2 

"Mendel son,  A,  "Plasticity:  Theory  and  Applications",  Macmillan  Publ.  C 
New  York,  1968. 
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FIGURE  3.5  Arrangement  of  heat  source  images  for  a 
finite  width  and  thickness  plate 

23  24 

used  in  the  solution  of  welding  problems  by  Tall  ,  and  later  by  Masubuchi  . 

To  analyze  the  stress  state  at  the  center  cross  section  of  the  plate 
(Figure  3.3)  due  to  an  arbitrary,  and  changing  with  time,  temperature 
distribution,  0(y,t),  it  is  assumed  that  at  time  t,  the  section  is  a  part 
of  an  infinitely  long  plate  subject  to  the  same  temperature  distribution 

23 

Tall,  L.,  "Residual  Stresses  in  Welded  Plates,  A  Theoretical  Study", 
Welding  Journal,  4_3  (1),  1964  ,  pp.  10s-23s. 

24 

Masubuchi,  K. ,  Simmons,  F.B.,  and  Monroe,  R.E.,  Analysis  of  Thermal 
Stresses  and  Metal  Movement  During  Welding",  RSIC-820,  Redstone  Scientific 
Information  Center,  Redstone  Arsenal,  Alabama,  July,  1968. 
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over  its  entire  length,  as  in  Figure  3.6.  This  temperature  profile  will 
remain  the  same  during  the  current  time  increment,  At. 
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or 


o  =  E(e  -  a*A0  -  eln  -  Aein) 

XX  XX 


(3.16) 


where 


a  / E  =  Elastic  part  of  strain,  e 


el 


a'A0  =  Thermal  strain,  e 


th 


A6  =  0  - 


m 


pi 


=  Accumulated  (during  the  previous  time  increments) 

pi  c 

=  e  +  e 


inelastic  strain 
Plastic  strain 


Ac 


x 

in 


=  Creep  strain 


=  Change  in  inelastic  strain  during  the  time  increment  At 


From  global  equilibrium  (no  external  forces  and  moments  acting  on 
the  plate) 

*+c 


L 

r: 


o  dy  =  0 
x 


a  y  dy  =  0 

x 


(3.17a) 


(3.17b) 


Substituting  Eqs.  (3.14)  and  (3.15)  into  (3.17),  a  set  of  linear 
equations  is  obtained  for  the  determination  of  the  unknown  coefficients  c^ 
and  Solving  this  system  and  substituting  back  into  Eq .  (3.14),  the 

following  expression  is  obtained  for  the  total  strain: 


\(y)  =  (A1 


(A, 


y a2) 
yA3) 


E(»-A0 


E(u • A9 


+tln  +  At 1 y  dy 

X  X 


m  .  in. 

+  c  +  A  c  )y  dy 
X  x  y 


(3. 18) 
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where 


(3.19) 


Equations  (3.16)  and  (3.19)  are  not  enough  to  solve  the  problem. 
What  is  still  needed  is  a  stress-strain  law  and  a  relation  between 
stress  and  creep  strain  increments. 


To  proceed  further,  the  assumption  is  made  that  creep  will  only 
take  place  during  the  soaking  stage  of  the  temperature  history.  Thus, 
during  this  period  the  accumulated  plastic  strain,  ,  will  remain 
constant.  The  heating  and  cooling  stages  where  no  creep  occurs  are 
treated  in  exactly  the  same  way  as  the  welding  problem. 

(A)  Specifically  during  welding  when  creep  is  not  taken  into  account: 


in 

£ 

X 


(y) 


■  ex1<y) 


(3.20) 


At  each  time  step,  the  total  strain  is  first  calculated  along 
the  cross  section  from  Eq.  (3.18)  assuming  that  no  plastic  strain  exists. 

t^(y)  =  0  (3.21) 

The  mechanical  strain,  em,  then  is: 

£>)  =  r-x(y)  -  exh(y)  =  ex(y)  ~  a-AO(v)  (3.22) 


Now  assuming  a  bilinear  stress-strain  law,  a  first  approximation  of  the 
plastic  strain  along  the  cross  section  can  be  obtained,  as  in  Figure 
3.7.  This  value  can  be  used  again  in  Eq .  (3.18)  to  obtain  a  second 


FIGURE  3.7  Bilinear  stress  strain  law 
used  in  1-D  model 


FIGURE  3.8  Uniaxial  creep  curve 
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approximation  of  the  total  strain  and  the  process  can  be  repeated 
until  convergence  is  reached. 


Further  details  of  this  iterative  procedure,  which  can  also 
be  applied  during  the  heating  and  cooling  stages  of  a  heat  treatment, 
can  be  found  in  the  previously  mentioned  references  on  the  successive 
elastic  solutions  method.  However,  it  should  be  noted  that  during 
the  calculation  of  the  total  strains  at  each  time  step,  the  accumulated 
plastic  strains  from  previous  time  steps  should  be  included  to  account 
for  possible  elastic  unloading  or  reverse  yielding. 

(B)  During  the  holding  period,  when  creep  is  taken  into  account:,  this 
procedure  has  to  be  slightly  modified. 

At  the  start  of  the  first  time  increment  the  inelastic  strain  is: 


in.  v  pi,  . 

E*  <y>  ■  <  <*> 


(3.23) 


<  1 


where 

is  the  total  accumulated  plastic  strain  up  to  that  instant 
(due  to  welding  and  heating) 

As  mentioned  before,  will  remain  constant  during  the  whole 
soaking  period. 

To  get  a  first  approximation  of  the  total  strain,  e  ,  after  At, 
from  Eq .  (3.18)  it  is  now  assumed  that 


Ln  -  a  c  -  n 
—  At  -  U 


(3.24) 


This  first  approximation  of  the  total  strain,  e  ,  is  then  substituted 
in  Eq .  (.3.16)  to  obtain  a  first  approximation  of  stress  a^.  Using 

this  stress  approximation  and  the  appropriate  creep  law  (Section  3.3.5), 

,  .  .  r  ,  ■  •  „  c  ,  in 

a  second  approximation  for  the  creep  strain  inclement,  Ac^  =  , 

is  obtained.  This  value  is  now  again  substituted  in  Eq .  (3.18)  for 

a  new  approxim.it  ion  of  the  total  strain  c  and  the  process  is  repeated 

until  convergence  is  reached. 
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At  the  start  of  the  second  and  any  subsequent  time  increments, 
the  total  strains  will  be  equal  to  the  initial  plastic  strain  plus 
the  accumulated  creep  strain  during  the  previous  time  increments: 

n-1 

<n(y)  =  ex1(y)  +  Aex(y)  (3.25) 

.  c 

Assuming  again  that  Ae^  is  zero  during  this  time  step,  the  process 
outlined  above  can  be  repeated. 

3.3.5  Creep  Laws 

Creep,  the  time  dependent  deformation  and  fracture  of  materials, 
is  probably  the  most  general  type  of  material  behavior.  The  current 
state  of  the  art  in  inelastic  high-temperature  analysis  requires  that 
plasticity  and  creep  constitutive  equations  be  formulated  largely  on 
independent  bases.  However,  elevated-temperature  deformation  is, 
essentially,  the  result  of  time-dependent  processes  where  uoth  plastic 
.is.c  e ;  . -v  I  citaviors  are  present  simultaneously.  Prior  creep  deform¬ 
ations  influence  subsequent  plastic  behavior  and  vice  versa.  However, 

onlv  limited  information  is  available  on  these  mutual  interactions  as 

25 

outlined  by  Corum,  et  al .  For  the  purposed  of  this  study  the  two 
behaviors  were  modeled  separately,  as  already  noted  in  the  previous 
section . 

Very  limited  information  is  available  in  the  literature  on  the 

creep  behavior  of  high-strength,  quenched  and  tempered  steels,  as  HY-80 

and  HY-L30.  Only  some  data  on  the  minimum  creep  rupture  time  are 

26 

reported  (as  by  Domis  ). 

5 

Corum,  J.M.,  et  al..  Interim  Guidelines  for  Detailed  Inelastic 
Analysis  of  High-Temperature  Reactor  System  Components",  Oak  Ridge 
National  Laboratory  Report  ORNL-5014,  Dec.  1974. 

6 

Domis,  W.F.,  'Creep  and  Creep  Rupture  Properties  of  HY-80  and  HY-130 
(T)  Steels",  l' . S .  Steel  Applied  Research  Laboratory,  Report  No. 
19.012-006  ( 1 ) ,  July  15,  1963. 
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For  stainless  steels,  on  the  other  hand,  numerous  studies  have 

been  performed  to  investigate  their  elevated-temperature  inelastic 

27 

behavior.  As  reported  by  Clinard,  et  al.  ,  uniaxial  creep  behavior 
(Figure  3.8)  of  stainless  steels  during  the  primary  and  secondary 
stage  can  be  very  well  modeled  by  an  equation  of  the  form: 

e^o.t.T)  =  f(cr,T)  ri-e_r(ax,T)tl  +  g(a.T)  .  t  (3.26) 

A  A  A  I  I  A 

where  eL  is  the  uniaxial  creep  strain,  a  the  applied  stress,  T  the 

X  X 

test  temperature  and  t  the  time. 

The  functions  f (a  ,T),  r(a  ,T)  and  g(o  ,T)  can  be  deduced  from 
xx  x  27 

creep  test  data  by  curve  fitting.  Clinard,  et  al.  ,  report  for  304 
stainless  steels,  the  following  representation  at  T  -  1100°F  (594°C). 

v  ,  ..-5  1.843 

i  (ij  )  =  5 .436x10  o 
x  x 

r ( 0  )  =  5.929xl0~5  exp(0.2029o  )  (3.27) 

.3.0 

sinh  (0.14790^)1 

where 

a  is  expressed  in  Ksi,  t  in  hours  and  the  creep  strain 

and  in  inches/inch, 
x 

Creep  response  to  constant  uniaxial  compression  is  usually  assumed 
identical  to  that  in  tension  (actually  a  reflection  of  it  with  respect 
to  the  time  axis). 

For  tiie  complete  description  of  the  time  dependent  behavior  of 
the  material,  a  "hardening  rule"  is  needed  in  order  to  predict  the 
creep  response  when  the  stress  levels  are  changing.  The  two  most 
eommonlv  used  rules,  time  hardening  and  strain  hardening,  are  schemat¬ 
ically  shown  in  Figure  1.9.  It  should  be  pointed  out,  however,  that 

'7 

Clinard,  3.A.,  et  al.,  "Ver i i icat ion  by  Comparison  ol  Independent 
Computer  Program  Solutions",  Pressure  Vessels  and  Piping  Computer 
Program  evaluation  and  Qua  1 i i i cat  ion,  Fnergv  Technology  Conference , 
Houston,  Texas,  September  1977,  A.S.M.P.,  pp.  27-49. 


g(c  )  =  6.73x10 
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these  rules  are  strictly  correct  only  in  step  changes  of  stress  that 
are  of  long  duration. 
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FIGURE  3.9  Strain-hardening  and  time-hardening  models 

of  creep  response  under  a  stepwise  varying  load 


Experimental  evidence  tends  to  support  a  strain  hardening 
formulation  for  the  case  of  304  and  316  stainless  steels.  In  our 
problem,  however,  where  small  changes  of  stress  take  place  at  each 
infinitesimal  time  step,  it  was  decided,  for  computational  efficiency, 
to  adopt  t  ime-hardening,.  Specifically  the  creep  strain  increment 
Ae^(v)  accumulated  between  time  t  and  t+  At  (at  each  point  along  the 
cross  section)  is: 


-41- 


e 


i 


i 

L 


Ae^  (y)  =  e^[ox(y),  (t+At),  T(y)j  -  [0x(y)  » t  »T(y)  ]  (3.28) 

where 

T(y)  =  the  temperature  distribution,  and 
a  (y)  =  the  current  approximation  of  the  stress  distribution 

In  Eq.  (3.28)  the  creep  strain  rate  Ae^(y)/At  is  a  function 
of  stress,  time  and  temperature.  Time  hardening  was  assumed  in 
order  to  avoid  the  added  computations  of  solving  Eq.  (3.26)  for  time, 
as  a  strain  hardening  formulation  would  require. 

3.3.6  Computer  Implementation  and  Results 

It  should  be  noted  at  this  point  that  in  the  calculation  of  the 
temperature  distribution  during  welding  and  non-uniform  heat  treating, 
a  special  iterative  scheme  is  employed  in  order  to  account  for  the 
temperature  dependence  of  the  material  properties. 

In  the  stress  calculation,  all  variables  and  formulas  are  non- 
dimensionalized  and  all  the  integrations  are  performed  numerically. 
Finally,  it  is  worth  pointing  out  that  the  previously  derived  equations 
are  considerably  simplified  in  special  cases.  So,  during  butt  or 
bead-on-plate  welding,  the  temperature,  and  thus  the  strain  and  stress, 
distribution  is  symmetrical  with  respect  to  the  weld  line.  Further¬ 
more,  during  uniform  heat  treating  the  temperature  is  assumed  to  be 
constant  over  the  whole  plate  during  each  time  step.  Further  details 
of  the  above  analysis,  a  complete  listing  of  the  modified  one¬ 
dimensional  program,  and  results  of  various  sample  cases  are  given  in 
Agapakis*  thesis. 

The  predictions  of  the  one-dimensional  program  for  the  case  of 
edge  welding  and  subsequent  stress  relieving  of  304  stainless  steel 
plates  are  presented  next  in  this  section.  Welding  conditions  were 
assumed  exactly  the  same  as  in  the  experiments  described  in  Section 
3.4.]  of  tiiis  report.  Temperatures,  strains  and  stresses  were  cal¬ 
culated  across  a  center  strip  of  the  specimen  throughout  welding  and 
stress  relieving  operations. 
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A  range  of  values  was  found  in  the  literature  for  the  arc  efficiency 
and  the  surface  heat  loss  coefficient.  Since,  no  experimental  measure¬ 
ments  of  these  parameters  were  made  in  this  study,  the  actually  selected 
values,  within  that  range,  were  such  as  to  minimize  the  deviation  of  the 
predicted  temperature  history  from  the  experimentally  measured  one. 

The  predicted  temperatures,  mechanical  strains  and  stresses  during 
welding  are  plotted  in  Figures  3.10,  3.11,  3.12.  For  ease  of  com¬ 
parison  with  experimental  data,  the  same  locations  (0.5,  1.0,  1.5,  2.0 
and  3.0  inches  from  the  weld  line)  are  selected. 

During  uniform  stress  relieving  at  1000°F,  the  assumed  temperature 
history  is  shown  in  Figure  3.13  and  the  predicted  variations  in 
mechanical  strain  and  stress  are  plotted,  versus  time,  in  Figures 
3.14  and  3.15.  Comparisons  of  the  predicted  residual  stress  distribu¬ 
tion  after  welding  and  after  stress  relieving  are  given  in  Figure  3.16. 

Analytical  results  for  HY-130  steel  are  given  in  Section  3.4.2, 
together  with  experimental  measurements. 

3. 4  Experimental  Study 

To  verify  the  analytical  predictions,  experiments  were  initially 
performed  using  304  stainless  steel  plates.  As  already  mentioned, 
this  was  decided  since  extensive  creep  data  for  HY-130  steel  could 
not  be  located.  Only  uniform-temperature  heat  treatments  were  tested 
in  this  first  phase. 

Subsequently,  further  experiments  were  conducted  using  HY-130 
steel  plates.  Soth  uniform  and  line  heating  treatments  were  tested. 

A  summary  of  the  experimental  procedures  and  results  will  be 
presented  in  the  next  few  parts  of  this  section. 


i 


MECHANICAL  STRAINS  (X  1000)  WELDING  TEMPERATURES  (°F) 


420 


FIGURE  3.10  Temperatures  during  edge  welding,  as 

predicted  by  the  one-dimensional  progra 


FIGURE  3.11  Mechanical  strains  during  edge  welding,; 

predicted  by  the  one-dimensional  prograr 


STRESSES  (ksi ) 


REMAINING  STRESSES  (ksi) 


TRANSVERSE  DISTANCE  (in.) 


FIGURE  3.16  Comparison  of  residual  stresses  before, 
during  and  after  stress  relieving  at 
1100°F,  as  predicted  by  the  one-dimensional 
program 
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3.4.1  Initial  Verification  Tests  on  304  Stainless  Steel  Plates 

28 

The  geometry  and  the  dimensions  of  the  four  main  specimens  are 
shown  in  Figure  3.17.  All  plates  were  edge  welded  and  all  but  one 
were  subsequently  stress  relieved  in  a  furnace  at  different  holding 
temperatures.  The  plates  were  finally  sectioned  and  residual  stresses 
were  measured  by  stress  relaxation. 

The  temperature  and  strain  changes  during  welding  and  subsequent 
heat  treatment  were  measured  by  thermocouples  and  electrical- 
resistance  strain  gages  attached  on  the  plates.  The  thermocouple  and 
strain  gage  locations  are  also  depicted  in  Figure  3.17.  The  total 
number  of  strain  gages  and  their  configuration  is  given  in  Table  3.1 
for  all  specimens. 


□  Strain  qage  locations 
,  Thermocouple  locations 


FIGURE  3.17  Specimen  geometry  and  instrumentation 

Five  more  specimens  were  used  during  this  phase,  in  order  (a)  to 
determine  the  appropriate  welding  conditions,  (b)  to  evaluate  the 
performance  of  the  data  acquisition  system,  and  (c)  to  test  the  heating 
uniformity  of  the  stress-relieving  furnace. 


Table  3.1:  Arrangement  of  Strain  Gages, and  Thermocouples 

for  the  Verification  Tests 


Specimen 

Strain  Gages 

Thermocouples 

Specimen  No. 

Dimensions 

Number 

Configuration 

Number 

Configuration 

1 

22"x6"x3/8" 

5 

A 

5 

C 

2 

22"x6"x3/8" 

A 

5 

C 

3 

22"x6"x3/8" 

10 

B 

5 

C 

4 

22"x6"x3/8" 

102 

B 

5 

C 

Key:  (A) : 

Single  gages  in  the  longitudinal  direction  on  the  one  side 

of  the  plate  only 

E 

(B)  : 

Single  gages  in  the  longitudinal  direction  on  both  sides 

of  the  plate  (to  determine  possible  bending) 

(C): 

Thermocouples  located  at  the  surface  of  the  specimen  on 

one  side  only 

I 

X,  1 

Notes : 

Five  additional  gages  were  added  in  the  transverse  direction 

before  cutting. 

5 

JOriginal  strain  gages  were  destroyed  during  stress  relieving. 
Ten  new  gages  were  installed  in  both  the  longitudinal  and 
transverse  directions  before  cutting. 


* 


I 


L 


J 
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Straight  polarity  short  circuiting,  gas  metal  arc  welding  was 
performed  on  all  specimens  utilizing  a  308  stainless  steel  wire 
(0.035"  in  diameter)  and  a  90%  He  -  7  1/2%  A  -  2  1/2%  C02  gas  mixture. 
The  finally  selected  welding  conditions  were  the  same  for  all 
specimens : 


Arc  Voltage 
Arc  Current 
Wire  Feed  Speed 
Weld  Travel  Speed 


23  volts 
90  Amperes 
175  inches/min 
16  inches/min 


Very  similar  temperature  and  strain  distributions  were  observed  on 
all  specimens.  This  is  evidenced  in  Figures  3.18  to  3.21,  which  depict 
the  temperature  and  strain  history  at  various  distances  from  the  weld 
line  in  specimens  #2  and  //4.  Based  on  this  similarity  in  the  distribu¬ 
tion  of  transient  temperatures  and  strains,  it  is  reasonable  to  assume 
that  very  similar  residual  stress  distributions  existed  in  all  the 
specimens  alter  welding. 

After  welding  and  cooling,  specimen  it 2  was  directly  cut  in  order 
to  determine  the  residual  stresses  in  the  as-welded  condition  whereas 
specimens  if  1,  #3  and  if 4  were  stress  relieved  at  different  holding 
temperatures  before  cutting.  Details  on  the  stress  relieving  conditions 
are  given  in  Table  3.2. 

The  temperature  history  during  stress  relieving  is  depicted  in 
Figures  3.22  and  3.23  for  specimens  it  1  and  it3  respectively.  Strain 
measurements  during  stress  relieving  were  only  taken  for  specimens  if  1 
and  v 3  since  the  gages  were  destroyed  at  temperatures  above  550°F. 
Strains  for  specimen  it  3  are  shown  in  Figures  3.24  and  3.25. 


Finally,  after  stress  relieving,  the  plates  were  sectioned  and 
the  obtained  distribution  of  residual  stresses  was  compared  with  the 
one  in  the  as-welded  case  (Figure  3.26).  More  details  on  these 
experiments  and  the  results  can  be  found  in  Agapakis'  thesis. 
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TIME  IN  SECONDS 


WELDING  OF  SPECIMEN  #2 

FIGURE  3. i,j  Thermocouple  readings  during  welding  of  specimen  #2 


WELDING  OF  SPECIMEN  =2 

FIGURE  3.19  Strains  during  welding  of  specimen  1‘2  ("corrected  for 
temperature- induced  apparent  strain  and  gage  factor 
vnriat ions. 
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/triucation  xests 


Specimen  No. 

Holding  Temperature  °F  (°C) 

Time  in  Furnace  hr 

1 

500°F  (260°C) 

6 

2 

not  stress  relieved 

- 

3 

approx.  370°F  (188°C) 

4 

4 

1100°F  (593°C) 

7 

6 

500°F  (260°C) 

2 

370°F  (188°C) 

I 


Note:  All  stress  relieved  specimens  were  cooled  in  air  after 

heating . 


THERMOCOUPLE  READINGS  IN  DEGREES  F  (front  side)  ~  THERMOCOUPLE  READINGS  IN  DEGREES  F 


iCTTO 


0.0  5000.  10C00.  15000.  20000.  25000. 

TIME  IN  SECONDS 

STRESS  RELIEVING  OF  SPECIMEN  #1 

UPC  3.22  Thermocouple  readings  during  stress- 
relieving  of  specimen  #1 


0.0  5000.  10000.  15000.  20000. 

TIME  IN  SECONDS 


STRESS  RELIEVING  OF  SPECIMEN  =3 


FIGURE  3.23  Thermocouple  readings  during  stress 
relieving  of  specimen  A 3 


TIME  IN  SECONDS 

STRESS  RELIEVING  OF  SPECIMEN  #3 


FIGURE  3.24  Strains  during  stress-relieving  of  specimen  #3,  front  side 
(corrected  for  temperature  induced  apparent  strain  and  gage 
factor  variations) 


FIGURE  3.25  Strains  during  stress  relieving  of  specimen  f/3,  back  side 
(corrected  for  temperature  induced  apparent  strain  and 
ga re  factor  variations) 

_ J 
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It  should  also  be  mentioned  at  this  point  that  considerable  effort 
has  been  recently  devoted  in  upgrading  our  experimental  data  acquisition 
systems.  A  microcomputer  (MINC-23,  manufactured  by  the  Digital 
Equipment  Corporation)  was  purchased  for  this  purpose,  using  funds  not 
related  to  this  research  project,  and  was  interfaced  to  the  signal 
conditioning  equipment  for  strain  gages  and  thermocouples.  Further 
details  on  the  data  acquisition  system  are  considered  outside  the 
scope  of  this  report,  however,  and  can  be  found  in: 

Agapakis,  J.E.,  "Fundamentals  of  Computer  Aided  Experimentation 
for  Welding",  Department  of  Ocean  Engineering,  M.I.T.,  December  1982. 

3.4.2  Further  Experiments  on  HY-130  Steel  Plates 

Four  HY-130  steel  specimens  (22"x6"xl/2")  were  tested  in  this 
phase.  Their  geometry  and  configuration  was  the  same  as  that  of  the 
stainless  steel  specimens  (Figure  3.17).  All  plates  were  edge  welded, 
and  all  bur  one  were  subsequently  stress  relieved  (one  in  a  uniform 
temperature  furnace  and  the  other  two  by  line  heating  with  an  oxy- 
acetvlene  flame).  The  plates  were  then  cut,  ard  residual  stresses  were 
measured  by  stress  relaxation. 

Only  temperatures  were  measured  during  welding  and  stress  relieving 
operations.  Strain  gages  were  installed  on  all  specimens  (in  both 
the  longitudinal  and  the  transverse  direction)  only  before  cutting. 

Their  location  is  again  depicted  in  Figure  3.17.  Spray  transfer, 
gas  metal  arc  welding  was  performed  on  all  specimens  utilizing  an 
AX-140  (0.04S"  in  diameter)  wire  and  a  98%  A  -  2%  0^  gas  mixture.  The 
selected  welding  conditions  were  (same  for  all  specimens): 

Arc  Voltage  27  Volts 

Arc  Current  225  Amperes 

Wire  Feed  Speed  285  inches/min 

Weld  Travel  Speed  21  inches'/min 

The  temperature  distribution  during  welding  for  plate  #2  is  given 
in  Figure  3.27.  Very  similar  temperature  distributions  were  observed 
on  all  plates. 


I 
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After  welding,  specimen  //I  was  cut  in  order  to  determine  the  as- 
welded  residual  stress  distribution.  Specimen  if 2  was  heated  in  a 
furnace  at  1000°F'*'  and  air  cooled  to  room  temperature.  Specimens  #3 
and  #4,  on  the  other  hand,  were  flame  heated  using  an  oxyacetylene 
torch.  The  torch  was  moved  at  a  constant  speed  of  1  inch/min  along  a 
line  1.64  inches  (in  spec ’’men  if  3)  and  4.30  inches  (in  specimen  #4) 
away  from  the  weld  line.  The  temperature  distribution  during  stress 
relieving  is  given  in  Figures  3.28  and  3.29  for  specimens  if 2  and  if 4, 
respec  t i velv . 

To-  selection  of  the  torch  velocity  and  location  was  based  on 
computer  predictions  for  maximum  stress  relaxation.  The  resulting 
residual  stress  distributions  are  given  together  with  the  computer 
predictions  in  Figures  3.30  to  3.33. 

further  details  on  these  experiments  and  their  results  can  be 
found  in  (arpentier's  thesis.  It  should  be  noted  however,  that  due  to 
t!  iiod  number  of  experiments  no  definite  conclusions  can  be 

drawn  on  the  eff  'tiveness  of  line  heating  as  a  stress  relieving 
technique.  At  this  point,  further  investigation  is  considered  necessary. 


Since  c ret  p  ■  IP'.  i  tool  mid  not  he  handled  bv  our  rank  1  the 
plates  wo!,  I'.:'  in  tin.  ’a:  cot  oniv  tor  about  an  hour. 
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4.  WORK  ACCOMPLISHED  ON  TASK  5  -  IMPROVED  COMPUTER  PROGRAMS 

4.1  General  Status 

The  objective  of  Task  5  of  the  current  research  program  is  to 
develop  manuals  of  improved  computer  programs  for  analyzing  (1)  heat 
flow,  (2)  transient  thermal  stresses  and  metal  movement,  and  (3) 
residual  stresses  and  distortion  of  weldments.  Although  Task  5  includes 
some  efforts  at  improving  existing  computer  programs,  it  does  not 
include  major  improvements  of  old  computer  programs  or  the  development 
of  new  ones. 

A  unique  feature  of  Task  5  is  the  revision  of  the  programs  in 
such  a  way  that  they  can  become  more  useful  to  practicing  engineers. 
Moreover,  a  series  of  examples  to  aid  engineers  in  understanding  their 
use,  capabilities,  and  limitations  are  also  to  be  included. 

Two  subtasks  to  be  performed  under  Task  5  have  been  identified  in 
the  proposal  to  O.N.R.  dated  November  1980,  as  follows: 

Task  5.1:  Development  of  Improved  Manuals  of  Independent 
Welding  Programs 

Task  5.2:  Development  of  Manuals  of  Programs  which  are 
Compatible  with  the  ADINA  Program 

Details  on  the  work  done  are  included  in  the  next  two  sections. 
Note,  however,  that  all  manuals  are  published  separately  from  this 
final  report. 

4 . 2  Manila  is  of  Independent  Welding  Programs 

During  the  past  decade  M.I.T.  researchers,  under  the  direction  of 
Professor  K.  Masubuchi,  have  developed  a  series  of  computer  programs 

t 

for  analyzing  various  problems  related  to  heat  flow,  transient  thermal 
strains,  residual  stresses,  and  distortion  in  weldments.  On  some 
computer  programs,  manuals  have  already  been  prepared,  some  of  which 
were  included  in  the  monograph  prepared  under  Task  1  of  a  previous 
O.N.R.  Sponsored  Research  Program.  Since  then,  largo  amounts  of 
experimental  data  have  been  obtained  and  compared  with  analytical 
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predictions.  Considerable  efforts  have  been  made  at  the  same  time  to 
improve  these  programs. 

Among  these  programs,  the  one-dimensional  program  was  proven  to 
be  very  useful  for  the  case  of  welding  thin  plates  and  T-shaped  built- 
up  beams.  As  a  consequence,  it  was  decided  that  a  new  manual  for 
this  program  needs  to  be  written.  This  new  manual  is  published 
separately : 

Imakita,  A.,  Papazoglou,  V.J.,  and  Masubuchi,  K. ,  "Consolidated 
Manual  of  One-Dimensional  Computer  Programs  for  Analyzing  Heat  Flow, 
Transient  Thermal  Strains,  Residual  Stresses,  and  Distortion  in 
Weldments",  prepared  under  Contract  N00014-75-C-0469  (M.I.T.  OSP 
#  82558)  for  the  Office  of  Naval  Research,  M.I.T. ,  July  1982 

In  what  follows  a  summary  of  the  modifications  made  in  these 
programs  is  described.  For  more  details,  one  is  referred  to  the  afore¬ 
mentioned  manual. 

One-Dimensional  Program.  This  computer  program  is  capable  of 
calculating  the  temperatures,  thermal  strains,  and  stresses  developed 
in  a  thin  plate  during  and  after  welding,  assuming  a  one-dimensional 
state  of  stress.  The  following  modif ications, as  compared  to  the  earlier 
version  of  the  program, were  made  prior  to  developing  the  revised  manual: 

(1)  The  numerical  integration  method  was  modified 

(2)  The  number  of  integration  points  was  increased  to  improve 
accuracy 

(3)  The  temperature  dependence  of  the  surface  heat  loss 
coefficient  was  taken  into  account,  and 

(4)  The  output  format  was  changed  to  facilitate  easier  handling 
of  the  results 

A  detailed  sample  analysis,  accompanied  by  a  number  of  comments, 


is  also  included  in  the  manual. 
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One-Dimensional  Program  for  T-shaped  Beams.  This  computer  program 
is  an  extension  of  the  previous  one  and  is  aimed  at  calculating  tem¬ 
peratures,  thermal  strains,  residual  stresses,  and  longitudinal  bending 
distortion  during  the  welding  fabrication  of  built-up  T-beams.  The 
following  modifications  were  made  in  this  program: 

(1)  The  temperature  dependence  of  material  properties  was  used 
in  conjunction  with  a  parabolic  interpolation  instead  of  the 
linear  one  previously  used 

(2)  The  heat  flow  between  the  web  and  the  flange  plates  was 
taken  into  account 

(3)  The  strain  calculation  procedure  was  changed  to  conform 
with  the  one  used  in  the  simple  one-dimensional  program 

(4)  The  program  is  now  capable  to  also  perform  calculations  for 
the  case  of  edge  welds,  and 

(3)  The  output  data  format  was  changed,  while  at  the  same  time 
more  results  can  be  printed  out 

The  manual  includes  detailed  sample  analyses  of  both  edge  and 
T-beam  welds. 

No  efforts  have  been  or  will  be  undertaken  to  develop  manuals 
for  the  two-dimensional  finite  element  programs  that  have  been  developed 
independently  in  the  past  by  M.I.T.  investigators.  This  decision  was 
based  on  the  availability  of  the  much  more  sophisticated  and  efficient 
programs  ADINAT  and  ADINA. 

4.3  Manuals  of  Programs  Compatible  to  ADINA 

During  the  initial  stages  of  the  current  research  project,  it  was 
decided  that  all  computer  programs  developed  should  be  compatible  with 
a  general  purpose  finite  element  computer  program,  to  facilitate  easier 
integration  with  other  kinds  of  numerical  analyses,  e.g.,  structural 
static  and  dynamic,  fracture,  etc.  Upon  consultation  with  represent¬ 
atives  of  the  Office  of  Naval  Research,  the  general  purpose  finite 
element  codes  ADINAT  and  ADINA,  capable  of  performing  nonlinear  heat 
transfer  and  the rma 1-e last ie-p last ic  stress  analyses  respectively  ,  were 
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chosen.  Both  of  these  codes  were  developed  by  Professor  K.  J.  Bathe 
and  co-workers  at  the  Department  of  Mechanical  Engineering,  M.I.T. 

To  make  these  programs  capable  of  analyzing  the  welding  problem, 
i.e. ,  heat  flow  and  stress  calculations  during  welding,  several 
modifications  had  to  be  made.  Details  of  these  modifications  were 
presented  in  the  third  progress  report  to  O.N.R.  Efforts  were  made 
to  publish  manuals  that  explain  in  detail  the  programs  developed  and 
modified.  These  manuals  also  contain  directions  for  the  more  efficient 
use  of  the  programs  as  well  as  some  sample  analyses. 

More  specifically,  the  following  manual  is  completed  and  published 
independently : 

Imakita,  A.,  "Heat  Source  Models  Using  Finite  Element  Method  for 
Arc  Welding",  Dept,  of  Ocean  Engineering,  M.I.T. ,  July  1982. 

This  manual  gives  details  on  the  use  of  the  general  purpose  finite 
element  program  ADINAT,  for  the  calculation  of  temperature  distribution 
during  welding.  A  review  of  past  work  in  the  area  is  also  included 
togeLnor  with  some  sample  cases  and  discussion  on  the  effects  of 
various  model  parameters. 

Two  more  manuals  are  currently  in  the  final  stages  of  preparation: 

Papazoglou,  V.J.,  "Computer  Programs  for  the  Calculation  of 
Schematic  Continuous  Cooling  Transformation  Diagrams",  report  to  the 
Office  of  Naval  Research  (under  preparation). 

Papazoglou,  V..J. ,  "Numerical  Calculation  of  Microstructure  History 
and  Phase  Transformation  Strains",  report  to  the  Office  of  Naval 
Research  (under  preparation). 

The  first  one  describes  the  programs  developed  to  calculate  schematic 
continuous  cooling  transformation  (CCT)  diagrams  from  isothermal  data, 
i.e.,  from  experimentally  obtained  time-temperature-transformation 
(TTT)  diagrams. 

The  second  describes  the  program  developed  to  calculate  the 
microstructure  history  of  a  material  resulting  from  a  given  temperature 
historv  and  based  on  an  idealized  CCT  diagram  of  the  material.  It  also 


5.  PUBLICATIONS  AND  DEGREES  GRANTED 


5.1  Publications 


Work  done  under  the  present  contract  from  the  Office  of  Naval 
Research  has  resulted  in  the  following  publications: 

(1)  Masubuchi,  K.,  "Models  of  Stresses  and  Deformation  Due  to 
Welding:  A  Review",  presented  at  the  Conference  on  Modeling  of 
Casting  and  Welding  Processes,  Franklin  Pierce  College,  Rindge, 

New  Hampshire,  August  4-8,  1980. 

(2)  Masubuchi,  K.,  and  Papazoglou,  V.J. ,  "Thermal  Strains  and 
Residual  Stresses  in  Heavy  HY-130  Butt  Welds",  Proceedings  of 
1980  Fall  Meeting  of  SESA,  Ft.  Lauderdale,  Florida,  October 
12-15,  1980. 

(3)  Papazoglou,  V. J. ,  and  Masubuchi,  K. ,  "Analytical  Methods  for 
Determining  Temperatures,  Thermal  Strains,  and  Residual  Stresses 
Due  to  Welding",  presented  at  the  1980  ASM  Materials  and  Processes 
Congress ,  Cleveland,  Ohio,  October  28-30,  1980. 

(4)  Papazoglou,  V.J.,  and  Masubuchi,  K. ,  "Analysis  of  Thermal 
Strains  and  Residual  Stresses  in  High  Strength  Steel  Weldments", 
presented  at  the  62nd  Annual  Meeting  of  AWS ,  Cleveland,  Ohio, 

April  6-9,  1981;  submitted  for  publication  in  the  Welding  Journal. 

(5)  Masubuchi,  K. ,  "Welding  Stresses",  presented  at  the  28th 
Sagamore  Army  Materials  Research  Conference  on  Residual  Stresses 
and  Stress  Relaxation,  Lake  Placid,  New  York,  July  13-17,  1981. 

(6)  Masubuchi,  K.,  and  Agapakis,  J.,  "Analysis  and  Control  of 
Residual  Stress,  Distortion  and  their  Consequences  in  Welded 
Structures",  Proceedings  of  Conference  on  Trends  in  Welding 
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Numerical  Analysis  of  Thermal 
v;i22£  Stresses  During  Welding  Including 
Phase  Transformation  Effects 

K.  Masubuchi 

Professor,  The  problem  of  determining  temperature  distributions,  transient  thermal  strains, 
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Department  of  Ocean  Engineering  process  is  solved  using  the  finite  element  method.  First,  a  nonlinear  heat  transfer 

Massachusetts  intitule  of  Technology!  analysis  is  performed  taking  into  account  the  temperature  dependence  of  the 

Cambrige.  Mass  02139  material  properties,  and  convection  and  radiation  surface  heat  losses.  This  is 

followed  by  a  thermo-elastic-plastic  stress  analysis  that  incorporates  phase  trans¬ 
formation  strains.  Finally,  the  theoretical  predictions  are  compared  with  ex¬ 
perimentally  obtained  data  showing  good  correlation. 


Introduction 

The  welding  processes  have  been  and  are  extensively  used 
for  the  fabrication  of  various  structures  ranging  from  bridges 
and  machinery  to  all  kinds  of  sea-going  vessels  to  nuclear 
reactors  and  space  vehicles.  This  is  the  case  because  of  the 
many  advantages  they  offer  compared  to  other  fabrication 
techniques  excellent  mechanical  properties,  air  and  water 
tightness,  and  good  joining  efficiency,  to  name  but  a  few. 

At  the  same  time,  however,  welding  creates  various 
problems  of  its  own  that  have  to  be  solved.  In  the  past  these 
problems  have  been  tackled  through  experimental  in¬ 
vestigations  resulting  in  the  collection  of  extensive  data  on 
low -carbon  steel  using  conventional  arc  welding  techniques. 
Because  of  the  large  amount  of  manpower  and  the  high  cost 
required  for  such  investigations,  however,  this  method  is 
rather  inefficient  when  different  materials  (high-strength 
quenched  and  tempered  steels,  aluminum  and  titanium  alloys, 
etc.)  or  modern  welding  techniques  (electron  beam,  laser,  etc.) 
need  to  be  studied.  Fortunately  enough,  since  the  advent  of 
the  computer,  a  new  tool  became  available  which,  combined 
with  limited  experimental  efforts,  could  give  solutions  to 
many  of  the  welding  problems. 

One  of  the  problems  suitable  for  computer  analysis  is  the 
one  involving  residual  stresses  and  distortion  during  welding. 
The  uneven  temperature  distribution  produced  during  welding 
gives  rise  to  incompatible  strains  which  in  turn  result  in  self- 
equilibrating  residual  stresses  that  remain  in  the  structure 
after  it  has  cooled  down  to  ambient  temperature.  The 
phenomena  involved  are  very  complex  necessitating  the  use  of 
numerical  techniques  for  their  solution.  Moreover,  the 
prediction  of  residual  stresses  is  of  paramount  importance 
since  these  stresses  adversely  affect  the  service  behavior  of 
welded  structures  including  their  brittle  fracture,  stress 
corrosion  cracking,  fatigue,  and  buckling  characteristics. 
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Several  investigators  have  contributed  towards  the  solution 
of  the  welding  residual  stress  prediction  problem  in  the  past 
ten  or  so  years  [1-12J.  In  most  of  these  studies  use  is  made  of 
the  powerful  finite  element  method  in  two  dimensions  (plane 
stress,  plane  strain,  and  axisymmetric  analyses).  Moreover, 
the  thermomechanical  problem  is  assumed  to  be  uncoupled, 
with  the  thermal  and  mechanical  parts  being  treated 
separately  [2],  Good  correlation  is  generally  observed  between 
the  obtained  results  and  experimental  data.  In  cases  where 
poor  agreement  is  found,  the  discrepancy  is  attributed  to  the 
fact  that  no  consideration  is  given  to  the  phenomenon  of 
allotropic  phase  transformations  [2,  5). 

This  paper  describes  a  technique  for  analyzing  tem¬ 
peratures,  thermal  stresses,  and  residual  stresses  using  the 
finite  element  method.  Particular  emphasis  is  placed  on  the 
rational  determination  of  the  phase  transformation  strains. 
The  computed  results  are  then  compared  with  experimental 
data  obtained  during  butt  welding  25.4-mm-thick  HY-130 
plates  using  the  multipass  Gas  Metal  Arc  Welding  (GMAW) 
process. 

Heal  Transfer  Analysis 

The  problem  of  heat  How  during  welding  was  solved  using 
the  multipurpose  finite  element  program  ADINAT  developed 
by  Bathe  and  co-workers  at  the  Department  of  Mechanical 
Engineering,  M.I.T.  [13],  Following  is  a  discussion  of  the 
particular  features  of  the  program  pertaining  to  the  welding 
analysis. 

Finite  Element  Formulation.  The  governing  isoparametric 
finite  element  equations  for  the  nonlinear  heat  transfer 
problem  can  be  written  in  matrix  form  as  follows: 

tK  AlT1  =  t  +  aATq  ( i -  1)  (I) 

with 

r  +  ,*±t,ii)  =r  +  «Ar,( I  -  \  )  +  ±8"'  (2) 

Equation  (I)  represents  the  heat  How  equilibrium  at  time 
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Fig.  1  Weldment  configuration 


/  +  <*•  At,  where  0<a<  1  and  a  is  chosen  to  obtain  optimum 
stability  and  accuracy  in  the  solution.  Furthermore,  the 
solution  using  these  equations  corresponds  to  the  modified 
Newton- Raphson  iteration  scheme. 

All  boundary  conditions  relevant  to  the  welding  problem 
are  incorporated  in  the  matrices  of  equation  (1).  They  include: 

1  Convective  heat  losses  from  the  plates’  surfaces,  modeled 
according  to  Newton’s  law. 

2  Radiation  heat  losses  which  are  very  significant  in  the 
vicinity  of  the  weld  metal  due  to  the  large  difference  between 
the  surface  and  environmental  temperatures  and  which  are 
modeled  according  to  the  well-known  quartic  Stefan- 
Boltzman  law. 

3  The  heat  input  during  welding  is  modeled  using  the 
concept  of  arc  efficiency.  For  the  space  distribution  of  the 
heat  input,  a  consistent  formulation  is  adopted  assuming  a 
uniform  distribution  over  the  top  of  each  weld  bead.  Finally, 
the  time  distribution  of  the  heat  input  is  modeled  in  such  a 
way  that  the  passing  of  the  arc  over  the  cross  section  examined 
can  be  simulated  (a  linear  increase  as  the  arc  approaches, 
uniform  as  the  arc  travels  over  the  cross  section,  followed  by  a 
linear  decrease). 

A  family  of  one-step  methods  is  considered  for  the  time 
integration  using  the  parameter  a  [14j.  The  scheme  was  found 
to  be  unconditionally  stable  for  a  >  1/2  and  to  generally  give 
better  solution  accuracy  when  a  =  1  (Euler  backward 
method).  Finaly,  the  modified  Newton  iteration  is  guaranteed 
to  converge  provided  the  time  step  At  is  small  enough. 

Two-Dimensionality  of  Problem.  Assuming  that  the 
welding  heat  source  (arc)  is  moving  at  a  constant  speed  along 
a  regular  path  (e.g.,  a  straight  line  in  a  planar  weld),  that  the 


weld  speed  is  sufficiently  high  relative  to  the  material’s 
characteristic  diffusion  rate,  and  that  end  effects  resulting 
from  either  initiation  or  termination  of  the  heat  source  can  be 
neglected,  the  three-dimensional  character  of  the  heat  transfer 
problem  can  be  simplified  [3,  6,  15).  This  can  be  achieved  by 
analyzing  a  cross  section  of  the  weldment  of  unit  thickness 
and  located  in  the  mid-length  region  of  the  weld  (see  Fig.  1 ). 

Modeling  of  Weld  Metal.  The  heat  transfer  mechanism  in 
the  weld  metal,  when  molten,  is  extremely  complex  and  its 
physics  are  not  well  understood  as  of  today.  These  com¬ 
plexities  arise  from  the  difficulties  in  adequately  modeling  the 
convective  motion  of  the  molten  metal,  the  electric  heating 
due  to  the  current  flow  in  the  base  metal,  the  boundary 
conditions  for  heat  losses,  etc. 

Although  some  progress  has  recently  been  made  towards  a 
proper  modeling  of  the  fluid  flow  and  the  convective  heat 
transfer  [16],  these  efforts  are  still  directed  towards  simple 
cases  and  are  rather  inconclusive.  As  a  consequence,  most 
investigators  try  to  simulate  the  convective  heat  transfer 
mechanism  in  the  molten  metal  by  assuming  a  higher  [3,  6)  or 
lower  [2]  value  for  the  thermal  conductivity  than  that  for  the 
material  at  the  solidus  temperature.  In  this  study  a  higher 
value  is  used,  since  it  is  believed  that  it  better  approximates 
the  actual  phenomena. 

Finally,  since  the  cases  analyzed  involve  multipass  welding, 
it  is  necessary  to  find  a  way  to  model  the  laying  of  the  various 
beads  during  the  welding  cycle.  This  is  made  possible  by  the 
element  birth-and-death  capabilities  of  the  code  used. 

Transient  Strain  and  Residual  Stress  Analysis 

Using  the  predicted  temperature  distributions  obtained 
from  the  heat  transfer  analysis,  the  transient  stresses,  and 
residual  stresses  can  be  calculated.  In  the  present  study  these 
calculations  were  performed  using  the  multipurpose  finite 
element  nonlinear  stress  analysis  computer  program  ADINA 
(17),  properly  modified  to  take  into  account  the  trans¬ 
formation  strains  [15).  Details  of  the  method  used  are 
prt  e  ited  in  the  forthcoming. 

Finite  Element  Formulation.  The  thermal-elastic-plastic 
constitutive  model  (with  the  phase  transformation 
modifications  discussed  later)  is  used  to  study  the  mechanical 
part  of  the  welding  problem  for  the  case  of  plane  strain.  The 
governing  incremental  finite  element  equations  can  be  written 
as 

tkAU1'*  =  r  +  Ar„ -t  +  Aty  (»-  1)  (3) 

The  solution  using  equation  (3)  corresponds  to  the  modified 
Newton-Raphson  iteration  procedure  which  is  helpful  in 
improving  the  solution  accuracy  and  in  many  cases  in 
preventing  the  development  of  numerical  instabilities.  The 
convergence  of  the  iteration  can  be  accelerated  using  the 


Nomenclature 


B  =  bainite  transformation 
temperature 

<?  -  effective  plastic  strain 
e,  =  thermal  strain,  mixture  of 
phases 

ez  =  transformation  strain 
e,  =  thermal  strain,  single  phase 
/  =  fraction  transformed  at  8 
fB  =  fraction  transformed  at  9, 

F  =  vector  of  nodal  point  forces 
h  =  heat  convection  coefficient 
K  =  stiffness  matrix 
M  =  martensite  transformation 
temperature 


Q  =  effective  heat  flow  vector 
R  =  external  force  vector 
I  =  time 

=  bainite  transformation 
starting  time 

U  =  displacement  vector 
v  =  welding  speed 
x,y,z  =  Cartesian  coordinates 
a  =  time  integration  parameter 
AS  =  temperature  increment 
AU  =  displacement  increment 
At,,  =  bainite  transformation  time 
interval 

=  Kroneckcr  delta 


=  thermal  and  transformation 
strain 

i)„  =  arc  efficiency 
S  =  temperature  vector 
0  =  temperature 

Subscripts 

5  =  start 
/  =  finish 

Superscripts 

left  =  time  at  which  a  quantity 
occurs 

right  =  iteration  counter 
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Fig.  3  Typical  dialatatlonal  curve  Including  phase  transformation 
effects 

Ailken  method  or,  in  complex  material  nonlinear  cases  (like 
the  welding  problem),  improved  using  the  BFGS  (Broyden- 
Fletcher-Cioldfarb-Shanno)  matrix  updating  method. 

The  usual  small  strain  assumption  of  decomposing  the  total 
strain  at  any  time  instance  in  its  elastic ,  plastic,  and  thermal 
components  is  used.  Creep  is  neglected,  given  that  in  the 
welding  problem  the  time  intervals  at  high  temperature  are 
short.  The  plastic  strain  component  is  calculated  using  the  von 
Mises  yield  criterion,  modified  for  temperature  dependence, 
and  the  associated  flow  rule  Kinematic  hardening,  thought 
by  many  to  better  model  the  phenomena  involved,  is  used;  it  is 
based  on  the  assumption  that  the  size  of  the  yield  surface 
depends  on  temperature  only,  whereas  its  translation  rate  in 
the  stress  space  depends  on  the  plastic  strain  rate. 

Phase  Transformation  F.f feels.  In  most  previous  stress 
analyses  of  the  welding  problem  no  consideration  is  given  to 
the  phenomenon  of  phase  transformations.  The  limited  ef¬ 
forts  reported  in  the  literature  are  based  on  the  use  of  a  small 
number  (usually  three)  of  experimentally  obtained 
dilatational  curves  for  specific  cooling  rates  and  do  not 
concentrate  on  the  fundamentals  of  the  phenomena  involved 
(I0|  A  rational  model  developed  to  calculate  the  phase 
transformation  strains  is  discussed  in  this  section. 

The  model  starts  from  the  Continuous  Cooling  irans 
formation  (CCT)  diagram  of  the  material  under  consideration 
which  is  either  known  experimentally  or  can  be  derived  from 
isothermal  data  (T  IT  diagram).  Methods  for  the  latter  case 
are  reported  in  the  literature  by  considering  linear  cooling 
rates  [18.  19).  Based  on  the  tact,  however,  that  the  welding 
cooling  rates  resemble  log-linear  curves  some  modifications  to 
the  available  methods  had  to  he  made  to  take  this  f.ul  into 
account  [  1 5 1 

Csing  the  l  Cl  diagram  (see  Fig  3  tor  an  idealized  ver-ion 
ot  it)  and  the  temperature  distribution  during  welding,  the 


history  of  microstructure  formation  during  the  cooling  stages 
can  be  predicted  for  each  point  in  the  FIAZ  and  the  weld 
metal. 

To  predict  the  microstructures  during  the  cooling  stage  of 
the  welding  process  an  incremental  strategy  is  involved  so  that 
the  model  can  be  compatible  with  the  step-by-step  solution  of 
the  nonlinear  stress  analysis  using  ADINA.  A  procedure  has 
therefore  to  be  established  that  will  enable  one  to  calculate  the 
proceeding  of  each  transformation.  The  following  equation 
has  been  used  for  this  purpose: 

*3)’] 

Equation  (4)  is  applied  separately  for  each  allotropic  trans¬ 
formation,  namely  for  the  martensite  to  austenite  one  upon 
heating  and  for  the  austenite  to  bainite  and  austenite  to 
martensite  ones  upon  cooling  [15]. 

Based  on  the  microstructure  predictions  calculated  by  the 
methodology  outlined  in  the  foregoing,  the  combined  thermal 
and  transformation  strains,  e™,  developed  during  the  weld¬ 
ing  thermal  cycle  can  then  be  found  from 

e™=  (e,  +e2  +  <?3)*$„  (5) 

The  estimation  of  the  thermal  strains  e,  and  e,  requires 
knowledge  of  the  microstructures  present  and  of  the  average 
thermal  expansion  coefficients  of  austenite,  bainite,  and 
martensite.  Similarly,  to  calculate  the  transformation  strain, 
t>:,  one  has  to  also  know  the  transformation  strains  for  each 
separate  phase  change,  namely  from  martensite  to  austenite, 
from  austenite  to  bainite,  and  from  austenite  to  martensite. 

To  show  the  impact  the  inclusion  of  phase  transformation 
effects  has  on  the  thermal  strain  of  a  section  which  is  free  to 
expand  or  contact.  Fig  3  is  constructed.  The  solid  line  depicts 
for  a  specific  temperature  history  the  thermal  strain  during 
the  heating  stage  of  this  history,  including  the  phase  trans¬ 
formation  effects  (phase  change  from  martensite  to 
austenite).  This  is  the  usual  curve  assumed  as  input  in 
previous  stress  analyses  using  the  finite  element  method. 
Upon  eooling,  the  same  curve  was  usually  assumed,  thus 
neglecting  the  phase  changes  from  austenite  to  bainite  and/or 
martensite.  By  contrast,  the  intermittent  line  shows  clearly  the 
effect  the  expansion  accompanying  these  transformations  has 
on  the  thermal  strain. 

Welding  Model.  A  cross  section  of  the  weldment  at  its 
midlength  is  used  in  this  study  under  plane  strain  conditions. 
This  is  rationalized  by  the  fact  that  for  relatively  long  plates 
the  maximum  stresses  are  developed  in  this  region. 

I  he  boundary  conditions  used  in  the  analysis  allow  free 
expansion  of  the  weldment  in  the  transverse  direction  as  well 
as  bending  At  the  same  time  the  structure  should  be  properly 
restrained  to  eliminate  all  possible  modes  of  rigid  body 
motion;  otherwise  the  stiffness  matrix  will  not  be  positive 
definite 

■\n  important  point  is  the  accumulation  of  plastic  strains  in 
the  regions  that  become  molten  during  the  welding  cycle. 
When  the  temperature  reaches  the  liquids  these  plastic  straips 
are  physically  relieved,  starting  to  accumulate  again  when  the 
metal  solidifies.  Flu*  presence  ot  nonzero  material  properties 
above  the  liquidus,  however,  would  cause  the  plastic  strains 
not  onlv  to  continue  accumulating  but  also  to  reach  ar- 
tifk  i all v  high  values  owing  to  the  very  low  magnitude  of  the 
mechanical  properties  at  these  temperatures.  It  was, 
therefore.  necessary  to  modify  ADINA  by  imposing  a  total 
relief  of  plastic  strains  when  the  material  melts. 

I  malls,  one  of  the  most  important  decisions  an  analyst  has 
to  make  when  perloinnng  a  nonlinear  incremental  -tress 
ana!’- -iv  w  tin.  -oluiion  strategy  to  be  followed,  because  the 
i..  :u,v  and  1 K  .onvei pence  characleiislics  of  the  solution 
depend  vc-\  mu. n  r  1  his  n  espeoulh  true  lot  complex 
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Fig.  5(a)  Finite  element  mesh  used  in  analysis  of  welding  problem 


Fig.  5<6)  F  nite  element  mesh,  continued  (12.7  mm  to  the  right  of  weid 
centerline) 


Fig.  6  Comparison  of  finite  element  results  with  experimental  data 
(weld  pass  No.  1) 


Fig.  7  Comparison  of  finite  element  results  with  experimental  data 
(pass  No.  2) 


situations  involving  highly  nonlinear  material  behavior  like 
the  one  encountered  in  the  welding  problem.  In  this  study 
these  parameters  were  chosen  based  on  previous  experience 
and  on  a  number  of  numerical  experiments. 

(experimental  Program 

A  series  of  experiments  consisting  of  butt  welding  HY-130 
steel  plates  using  the  GMAW  process  were  performed,  (each 
plate  measured  305  mm  x  610  mm  x  25.4  mm  resulting, 
after  welding,  in  the  configuration  shown  in  Fig.  4.  The  weld 
joint  consisted  of  a  double- V  groove  with  a  60-deg  included 
angle  in  accordance  with  U.S.  Navy  specifications.  All 
welding  experiments  were  performed  with  the  plates  resting 
on  knife  edge  supports  in  order  to  better  simulate  the  com¬ 
plete!'  unrestrained  condition.  F  or  more  details  on  the  ex 
pcrimental  procedures  used  see  [15]. 

I ernperatures  and  strains  were  measured  throughout  the 
'■  Mine  operation  using  Chrome!/  Alumcl  adhesively  bonded 
•  -oiiplo  and  9()-dcg  Rosette  electric  resistance  strain 
.  i.  counted  on  the  plates'  surface  respectively,  and 
i  M  at  various  distances  from  the  weld  centerline  (15). 


Comparison  of  Numerical  and  experimental  Results 

The  finite  element  mesh  used  for  the  numerical  calculations 
is  shown  in  Fig.  5,  where  only  the  left  half  of  the  plate  and  the 
weld  metal  is  shown,  the  right  half  being  exactly  symmetric. 
This  mesh  represents  the  unit  thickness  cross  section  of  the 
plate  at  its  midlength.  Four  to  six-node  isoparametric 
quadrilateral  elements  were  used,  taking  special  care  not  io 
include  any  triangular  elements  in  the  mesh  since  they  afe 
considered  by  many  as  inefficient.  A  total  of  77  nodes  and  47 
elements  were  utilized  with  the  mesh  being  finer  near  the  weld 
centerline  because  of  the  large  gradients  of  temperature  and 
stresses  presert  in  this  region.  The  temperature  dependence  of 
all  material  properties  was  considered  in  the  analyses  (1 5], 

Temperature  Distributions.  The  heat  transfer  analysis  was 
performed  for  the  first  three  welding  passes  only.  A  total  of 
four  analyses  (A I  to  A4)  were  made  by  varying  the  arc  ef¬ 
ficiency  and  the  temperature  variation  of  the  heat  convection 
coefficient  ( 1 5 [ . 

Cases  A2,  A3,  and  A4  are  compared  with  the  ex¬ 
perimentally  obtained  results  in  F  ig.  6  for  a  point  on  the 
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Fig.  8  Accumulated  effective  plastic  strain  history  at  various 
distances  trom  the  weld  centerline 


plate's  top  surface  12.7  mm  away  from  the  weld  centerline 
and  for  the  first  welding  pass.  Considering  the  same  heat 
input  ( =  0.60),  little  difference  is  found  between  cases  A2 
and  A3.  This  is  due  to  the  fact  that  although  higher  values  for 
the  convection  coefficent  were  chosen  in  analysis  A3,  these 
values  were  not  high  enough  to  significantly  alter  the  heat 
losses  from  the  top  and  bottom  surfaces  of  the  plate  and 
consequently  the  temperature  distribution.  A  substantial 
increase  in  the  convection  coefficient  was  therefore  chosen  for 
case  A4;  at  the  same  time,  however,  an  increase  in  the  arc 
efficiency  was  made  (>/„  -  0.70)  to  partially  compensate  for 
the  higher  h,  and  thus  to  obtain  a  good  estimate  of  the 
maximum  temperature  rachcd.  As  seen  in  F  ig.  6  the  com¬ 
bination  of  values  used  in  ihis  latter  analysis  succeeded  in 
bringing  the  cooling  rate  much  closer  to  the  experimental  one. 
At  the  same  time  a  4-percent  overprediction  of  the  maximum 
temperature  is  observed.  The  difference,  however,  was  very 
small  so  that  any  further  analysis  was  not  felt  necessary. 

Similar  remits  were  found  at  other  points  of  the  plate,  and 
for  other  passes.  C  ase  A4  came  always  closest  at  matching  the 
experimentally  obtained  results  (see  F  ig.  7). 


DC 


Fig.  10  Comparison  of  experimentally  measured  transverse  transient 
strains  on  the  plate's  top  surface  with  numerical  predictions 

Strain  and  Stress  Distributions.  Due  to  the  high  costs  in¬ 
volved  in  performing  the  numerical  stress  analysis,  only  the 
first  welding  pass  was  thoroughly  analyzed.  The  results 
obtained  are  d.scussed  below. 

Figure  8  shows  the  accumulation  of  the  effective  plastic 
strain,  e~p,  for  the  first  welding  pass  at  three  points  located 
on  the  plate’s  midthickness  and  at  various  distances  away 
from  the  weld  centerline.  In  the  weld  metal  zone,  plastic  strain 
is  built  up  rapidly  as  the  material  yields  in  compression.  At  the 
melting  region  this  plastic  strain  is  completely  relieved,  to 
rapidly  reach  again  a  level  in  close  proximity  to  that  existing 
prior  to  melting  upon  solidification.  As  cooling  continues  and 
yielding  in  tension  progresses,  the  plastic  strain  steadily  in¬ 
creases  but  at  a  smaller  rate  than  before.  In  the  heat-affected 
zone,  plastic  strain  is  accumulated  during  heat-up  in  much  the 
same  way  as  in  the  weld  metal  zone  but  not  quite  as  rapidly. 
As  the  material  cools,  it  yields  in  tension  and  further  plastic 
strain  is  accumulated. 

Figure  9  shows  calculated  transient  longitudinal  stress 
distributions.  Compressive  stresses  exist  in  the  weld  metal 
prior  to  melting.  When  the  metal  is  in  its  molten  stage, 
negligible  compressive  or  tensile  stresses  were  calculated.  As 
cooling  commences,  tensile  stresses  start  appearing  in  the  weld 
metal.  These  stresses  then  build  up  to  the  residual  stress 
pattern  when  ambient  temperature  is  reached.  For  self- 
equilibrating  purposes,  compressive  stresses  exist  in  areas 
removed  from  the  weld  centerline.  Note  in  Fig.  9  the  effect  of 
phase  transformation,  causing  a  sudden  decrease  in  the  tensile 
stresses. 

Comparison  of  the  experimentally  measured  transverse 
transient  strain  history  at  a  point  located  on  the  top  surface  of 
the  plate  is  made  with  the  numerically  obtained  results  in  Fig. 
10.  The  correlation  is  relatively  good  if  one  takes  into  account 
the  various  assumptions  involved  in  modeling  the  complex 
welding  problem.  A  delay  in  the  transition  from  tensile  to 
compressive  and  again  from  compressive  to  tensile  strains  is 
observed  in  the  analysis.  The  same  delay  as  far  as  the  oc¬ 
currence  of  the  maximum  strain  is  also  exhibited.  It  is  believ  ed 
that  this  phenomenon  is  primarily  due  to  the  relative 
coarseness  of  the  finite  element  mesh,  and  the  complex 
loading  history  present  in  the  welding  problem. 

Conclusions 

This  paper  presented  a  methodology  for  predicting  tem¬ 
peratures,  transient  strains,  and  residual  stresses  due  to 
welding  using  the  finite  element  method.  Particular  emphasis 
was  placed  on  the  incorporation  of  phase  transformation 
effects  in  the  analysis,  an  important  parameter  when  welding 
quenched  and  tempered  steels.  Comparison  of  the  computed 
results  with  experimentally  obtained  data  showed  the  im¬ 
portance  of  considering  these  effects. 

Finally,  it  should  be  noted  that  the  analysis  presented  is  a 
rather  expensive  one.  It  is  thus  recognized  that  it  cannot  be 
used  in  everyday  practice.  The  authors  believe,  however,  that 
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benefits  can  accrue  in  the  case  of  manufacturing  critical 
structures  when  an  accurate  estimation  of  the  metallurgical 
microstructures  and  residual  stresses  that  result  from  the 
welding  operation  is  very  essential. 
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APPENDIX  B 

Residual  Stresses  Due  to  Welding:  Computer-Aided  Analysis  of 
Their  Formation  and  Consequences 

Vassilios  J.  Papazoglou,1  Associate  Member,  Koichi  Masubuchi,2  Member,  Edison  Gonsalves,3 
Visitor,  and  Akihiko  Imakita,4  Visitor 


Results  of  recent  investigations  at  M.l.T.  on  the  subject  of  residual  stresses  due  to  welding  are  sum¬ 
marized.  Part  1  deals  with  the  computer-aided  prediction  of  residual  stress  distributions.  Both 
closed-form  analytical  and  numerical  solutions,  the  latter  based  on  the  finite-element  method,  to  the 
problems  of  determining  temperature,  transient  strain,  and  residual  stress  distributions  are  present¬ 
ed.  Experimentally  obtained  data  are  then  compared  with  predictions  obtained  by  the  developed 
computer  programs  to  test  their  validity.  Guidelines  for  the  applicability  of  each  model  are  also  in¬ 
cluded.  Part  2  discusses  the  computer-aided  analysis  of  the  effects  of  residual  stresses  on  the  ser¬ 
vice  behavior  of  welded  structures.  A  finite-element  analysis  of  the  fracture  characteristics  of  high- 
strength  steel  weldments  is  described  in  some  detail  as  a  case  study.  The  obtained  results  are  also 
compared  with  experimental  data. 


Introduction 

The  PROBLEM  of  predicting  residual  stresses  due  to  welding 
has  long  been  recognized  by  ship  designers  and  fabricators  as 
very  important  but  at  the  same  time  as  a  very  difficult  one  to 
analyze.  This  difficulty  has  its  origin  in  the  complex  mecha¬ 
nism  of  residual  stress  formation  which  starts  from  the  uneven 
temperature  distribution  caused  by  the  intense,  concentrated 
heat  source  associated  with  all  fusion  welding  processes.  The 
incompatible  strains  that  are  formed  as  a  consequence  give  rise 
in  turn  to  self -equilibrating  stresses  that  remain  in  the  welded 
structure  after  it  has  cooled  down  to  ambient  temperature,  thus 
producing  the  so-called  residual  stresses 

To  effectively  analyze  the  problem  and  account  for  the 
nonlinear  phenomena  associated  with  it,  one  has  to  draw  ujxin 
knowledge  from  a  variety  of  scientific  disciplines,  including 
heat  transfer,  applied  mechanics,  numerical  analysis,  and 
materials  science.  On  top  of  that,  the  lengthy  calculations  re¬ 
quired  make  the  solution  prohibitive  unless  one  has  access  to 
a  fast  electronic  digital  computer.  It  is  therefore  no  coincidence 
that  researchers  became  actively  involved  in  the  field  only 
during  the  past  decade  or  so.  during  which  time  computers  and 
the  powerful  numerical  techniques  developed  in  association 
w  ith  them  have  started  to  lx1  more  widely  used 

A  logical  question  that  could  come  to  one's  mind  at  this  [xiint 
is  why  should  one  l>e  concerned  w  ith  the  prediction  of  residual 
stresses  in  welded  structures,  especially  since  it  is  such  a  for- 
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midable  problem.  The  answer  is  that  the  effect  they  have  on 
the  service  behavior  of  these  structures  can  be  detrimental. 
Brittle  fracture  can  occur  earlier  because  the  presence  of  re¬ 
sidual  stresses,  combined  with  any  external  loading,  can  sub¬ 
stantially  decrease  the  critical  flaw  size  of  the  weldment 
Compressive  residual  stresses  located  in  the  regions  some  dis¬ 
tance  away  from  the  weld  line  can  substantially  decrease  the 
critical  buckling  stress  of  a  structure,  especially  if  it  is  composed 
of  thin  plates,  causing  structural  instabilities.  Similar  effects 
occur  in  the  areas  of  fatigue  fracture,  stress-corrosion  cracking, 
hydrogen  embrittlement,  and  others. 

ft  thus  becomes  evident  that  the  practicing  engineer  would 
like  to  have  a  design  tool  to  account  for  residual  stresses  and 
their  effects.  More  specifically,  what  he  wishes  to  do  is  to 
change  design  and  fabrication  parameters,  such  as  plate 
thickness,  joint  design,  welding  conditions,  and  welding  se¬ 
quence.  so  that  the  adverse  effects  of  residual  stresses  (and  the 
associated  distortion)  can  he  reduced  to  acceptable  levels  from 
the  point  of  view  of  reliability.  It  is  generally  much  Irettor  to 
achieve  this  goal  during  an  early  stage  of  design  and  fabrication 
rather  than  confronting  the  problem  at  later  stages  of  fabrica¬ 
tion  This  is  especially  true  when  critical  structures,  such  as 
submarines  and  other  types  of  deejr-sea  submersibles.  are  to  !>e 
built  using  new  materials  (for  example,  HY-150  steel,  titanium 
alloys)  or  unconventional  welding  processes  (for  example, 
electron  Iream,  laser  beam)  or  both. 

One  of  the  most  comprehensive  sources  of  information  re¬ 
garding  the  problem  of  residual  stresses  and  their  consequences 
is  a  recently  published  book  b\  Masubuchi  j  1  i  ’  This  book, 
w  hich  can  lx-  of  great  value  to  both  designers  and  fabricators 
of  welded  structures,  represents  a  systematic  presentation  of 
know  ledge  on  the  subject  accumulated  over  the  past  '50  or  so 
years  mainlv  on  the  basis  of  extensive  parametric  exix'rimental 
investigations,  which  can  collectively  lx-  thought  of  as  '  first 
generation"  or  "pre-computer  age  methixlologies  These 
metlxxlologies  have  served  and  still  serve  their  piirjxisc  in  eases 
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Fig.  1  Systematic  approach  for  predicting  residual  stresses  due  to 
welding  and  their  consequences 


that  have  already  been  extensively  investigated,  for  example, 
joining  of  mild  steel  structures  using  conventional  arc  welding 
techniques.  Time  and  cost  requirements  in  today  s  competitive 
world,  however,  do  not  allow  them  to  be  used  w  hen  new  ma¬ 
terials  and  processes  have  to  Ik-  effectively  used  in  the  con¬ 
struction  of  seagoing  vessels  and  offshore  structures.  The  need 
therefore  exists  for  a  new,  more  versatile  and  efficient  com¬ 
puter-age  methodology. 

This  paper  discusses  such  a  methodology  based  on  the  ex¬ 
tensive  use  of  computers  in  a  very  sy  stematic  way  as  shown  in 
Fig.  1.  Part  1  of  the  paper  discusses  state-of-the-art  methods 
for  predicting  residual  stresses  tand  distortion)  due  to  welding; 
it  thus  corresponds  to  Analysis  l  between  Steps  0  and  1  in  Fig. 
1  The  methods  presented  range  from  relatively  simple  one- 
dimensional  analyses  to  more  complicated  ones  based  on  the 
finite-element  method  with  nonlinearities  Comparing  the 
obtained  results  w  ith  data  collected  through  series  of  experi¬ 
ments,  the  predictive  capabilities  of  the  various  methods  are 
finally  assessed. 

In  Part  2  of  the  paper  the  subject  of  analyzing  the  conse¬ 
quences  of  residual  stresses  on  the  service  behavior  of  w  elded 
structures  is  addressed  (Analysis  2  Ix-twc-en  Steps  I  and  2  in  Fig 
1 ).  To  demonstrate  how  the  computer  can  help  in  this  respect, 
details  of  the  special  case  of  analyzing  the  fracture  character¬ 
istics  of  weldments  are  presented  The  numerical  results  ob¬ 
tained  using  the  finite-element  method  are  finally  compared 
w  ith  exerimentally  obtained  data 

It  is  not  the  intent  of  the  authors  to  present  in  detail  all  the 
mathematics  involved  in  the  discussed  methodologies.  The 
interested  reader  is  referred  tor  this  pur|x>se  to  other  sources 
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The  emphasis  will  rather  lx-  placed  on  how  the  ship 


designers  and  fabricators  can  effectively  use  the  state-of  -the-art 
methodologies  available  today  lor  solving  the  residual  stress 
problem 


Pari  I 

Prediction  of  welding  residual  stresses 


In  this  part  of  the  pafier  the  problem  ol  predicting  residual 
stresses  due  to  welding  is  addressed  After  a  brief  introduction 
on  the  physical  as(iects  of  residual  stress  formation,  the  various 


methodologies  available  for  calculating  residual  stress  distri¬ 
butions  are  discussed.  Particular  emphasis  is  placed  on  the 
method  based  on  a  computer  simulation  of  the  welding  process 
as  it  relates  to  residual  stress  formation.  This  consists  of  first 
solving  the  heat-transfer  problem  and  then  utilizing  the  ob¬ 
tained  results  to  perform  the  stress  analysis. 

Throughout  the  discussion  special  reference  is  made  to  the 
prediction  capabilities  of  the  analyses  by  comparing  analytical 
results  with  experimental  data  collected  from  series  of  experi¬ 
ments  conducted  over  the  past  several  years  by  Massachusetts 
Institute  of  Technology  (M.I.T.)  investigators. 

Genera!  information  on  residual  stresses 

Residual  stresses  are  defined  as  these  stresses  that  exist  in  a 
body  if  all  extern; '  loads  are  removed.  Various  alternative 
technical  terms  have  been  used  and  can  still  be  found  in  the 
literature  that  refer  to  residual  stresses,  including  internal 
stresses,  initial  stresses,  inherent  stresses,  reaction  stresses,  and 
locked-in  stresses. 

Areas  in  which  residual  stresses  can  exist  vary  greatly  in  scale 
from  a  large  [xirtion  of  a  metal  structure  down  to  areas  mea¬ 
surable  only  on  the  atomic  scale  1 1 1.  Regarding  the  former  case, 
macroscopic  residual  stresses  can  occur  in  scales  ranging  from 
ships  heated  by  solar  radiation  from  one  side  to  a  small  area  on 
the  top  of  a  plate  where  grinding  took  place.  On  the  other 
hand,  microscopic  residua!  stresses  can  Ik-  pnxluced  on  the 
atomic  scale  during  a  phase  transformation  or  near  a  dislocation. 
Furthermore,  residual  stresses  can  also  lx-  classified,  according 
to  the  mechanism  of  formation,  to  those  pnxluced  by  structural 
mismatching  and  those  pnxluced  by  an  uneven  distribution  of 
nonelastic  strains  (including  plastic  and  thermal  strains). 

In  the  case  of  welding,  with  which  this  paper  is  concerned, 
residual  stresses  can  lx-  classified  as  lx-ing  macroscopic  in  scale 
and  as  lx-ing  produced  by  uneven  distributions  of  nonelastic 
strains  Finally .  it  should  lx-  noted  that  since  residual  stresses 
exist  without  any  external  loads  they  should  always  satisfy  force 
and  moment  equilibrium,  dial  is.  they  should  Ik-  self-equili¬ 
brating. 

Welding  residual  stress  formation 

To  physically  understand  how  residual  stresses  are  formed 
during  welding,  the  simple  case  of  a  straight  liead-nn-plate  weld 
w  ill  lx-  described  in  some  detail  j  1 1  Figure  2  shows  schemat¬ 
ically  the  changes  of  temperature  and  stresses  that  occur  during 
such  a  process  The  welding  are.  which  is  moving  at  a  speed 
e.  is  presently  located  at  the  origin  ().  as  shown  in  Fig.  2(a). 

Figure  2 (/<)  shows  the  tein|x-ralure  distribution  along  several 
cross  sections.  Along  Section  A- A,  which  is  ahead  of  the 
welding  arc  the  temperature  change  due  to  welding.  Ad,  is 
almost  zero.  Along  Section  R-R,  which  crosses  the  welding  are, 
the  tenqx-rature  change  is  extremely  rapid  and  the  distribution 
is  very  uneven  Along  Section  (-().  which  is  some  distance 
Ix-hind  the  welding  arc,  the  distribution  of  teni|x-rature  change 
is  as  shown  in  Fig  2i/»i  -V  Along  Section  I)-D,  w  hich  is  very 
tar  from  the  welding  arc,  the  temperature  change  due  to 
welding  again  diminishes 

Figure  2(r>  shows  the  distribution  of  stresses  along  these 
sections  in  the  x -direction,  n,.  Stress  in  the  (/-direction.  n,t.  and 
shearing  stress.  also  evist  in  a  two-dimensional  stress 
held 

Along  Sxx  lion  A  A.  thermal  stresses  due  to  welding  are  almost 
zero  The  stress  distribution  along  Section  R-R  is  shown  in  Fig. 
2(c'  2  because  molten  metal  w  ill  not  sup|x>rt  a  load,  stress 
underneath  the  welding  arc  is  close  to  zero  Stresses  in  regions 
a  short  distance  from  the  are  are  compressive.  Iiecause  the  ex¬ 
pansion  of  these  areas  is  restrained  bv  the  surrounding  metal 
w here  the  teiiqx-rat ui es  are  lower  Since  the  temperatures  of 
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Fig.  2  Schematic  representation  of  changes  in  temperature  and  stresses  during 

welding 


these  areas  are  high  and  the  yield  strength  of  the  material  low, 
stresses  in  these  areas  are  as  high  as  the  yield  strength  of  the 
material  at  corresponding  temperatures.  The  magnitude  of 
compressive  stress  passes  through  a  maximum  with  increasing 
distance  from  the  weld  or  with  decreasing  temperature. 
However,  stresses  in  areas  away  from  the  weld  aTe  tensile  and 
balance  with  compressive  stresses  in  areas  near  the  weld.  In 
other  words 

f  axdy  =  0 

across  Section  B-B.'1  Thus,  the  stress  distribution  along  Section 
B-B  is  as  shown  in  Fig.  2(c)-2. 

Stresses  are  distributed  along  Section  C-C  as  shown  in  Fig. 
2(r)-3.  Since  the  weld  metal  and  base  metal  regions  near  the 
weld  have  cooled,  they  contract  and  cause  tensile  stresses  in 
regions  close  to  the  weld.  As  the  distance  from  the  weld  in¬ 
creases.  the  stresses  first  change  to  compressive  and  then  become 
tensile 

Figure  2(c)-4  shows  the  stress  distribution  along  Section  I)-I). 
High  tensile  stresses  are  produced  in  regions  near  the  weld, 
while  compressive  stresses  are  produced  in  regions  away  from 
the  weld.  This  is  the  usual  distribution  of  residual  stresses  that 
remain  after  welding  is  completed 

The  cross-hatched  area,  M-VT,  in  Fig  2 (a )  shows  the  region 
where  plastic  deformation  occurs  during  the  welding  thermal 
cycle.  The  egg-shaped  region  near  the  origin  ()  indicates  the 
region  where  the  metal  is  melted. 

The  region  outside  the  cross-hatched  area  remains  elastic 
during  the  entire  welding  thermal  cycle 

Sources  of  residual  stresses — In  the  previous  section  the 
difference  in  shrinkage  of  differently  heated  and  cooled  areas 
of  a  welder!  joint  was  identified  as  the  primary  cause  of  residual 
stress  formation,  resulting  in  high  longitudinal  stresses,  rrt,  in 
the  weld  metal.  Similar  tensile  stresses,  at/,  arise  in  the 

This  neglects  the  effect  of  eu  and  r,„  on  the  equilibrium 
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transverse  direction,  too,  but  of  smaller  magnitude. 

This  situation  is  typical  for  the  case  examined.  If,  however, 
one  is  also  interested  in  the  through-thickness  distribution  of 
residual  stresses  or  in  materials  exhibiting  phase  transformations, 
or  both,  two  additional  residual  stress  sources  can  be  identified 
16j. 

One  source  is  the  uneven  cooling  in  the  thickness  direction 
of  the  weld  Surface  layers  of  the  weld  and  the  highly  heated 
areas  close  to  it  usually  cool  more  rapidly  than  the  interior, 
especially  in  the  case  of  thick  plates.  Thermal  stresses  thus  arise 
over  a  cross  section  which  can  lead  to  heterogeneous  plastic 
deformation  and  hence  to  residual  stresses.  These  “quenching" 
residual  stresses  are  expected  to  be  compressive  at  the  surface 
of  the  highly  heated  areas  and  to  self-equilibrate  with  the  tensile 
ones  in  the  inner  regions. 

The  other  source  of  residual  stresses  comes  from  the  phase 
transformations  that  occur  during  cooling;  in  the  case  of  steel, 
for  example,  austenite  is  transformed  into  ferrite,  bainite,  or 
martensite,  or  a  combination  of  them.  These  transformations 
an*  accompanied  by  an  increase  in  specific  volume,  causing  the 
material  Ireing  transformed  (in  the  weld  metal  and  the  heat- 
affected  zone)  to  tend  to  expand.  This  expansion,  however, 
is  hindered  by  the  cooler  material  which  is  not  living  trans¬ 
formed.  inducing  compressive  stresses  in  the  transformed 
material  and  tensile  stresses  in  the  other  regions 

The  total  residual  stresses  due  to  welding  can  thus  lie  found 
bv  combining  the  effects  of  the  aforementioned  three 
sources 

Methods  for  predicting  residual  stresses 

The  incompatible  nonelastic  strains  produced  in  the  weld¬ 
ment  as  a  result  of  the  nonuniform  tem|>eraturc  distribution 
are  formed  in  a  very  complex  manner  as  discussed  in  previous 
sections  The  computational  efforts  required  to  analyze  the 
phenomena  involved,  including  high-temperature  plasticity 
calculations  in  a  multidimensional  stress  space,  have  limited 
investigators  efforts  in  the  past  to  such  simple  cases  as 
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Fig.  3  Weldment  configuration 

1.  Spot  welding  in  which  temperature  and  strain  changes 
are  axisymmetric. 

2.  Instantaneous  heating  along  the  edge  of  a  strip  in  which 
temperature  and  stress  changes  occur  in  one  dimension  only. 

In  more  complex  cases,  however,  and  before  the  advent  of 
the  era  of  powerful  digital  computers,  a  method  had  to  be  de¬ 
viser!  for  the  prediction  of  residual  stresses.  Such  a  method  was 
developed  by  Masubuchi  in  the  late  1950’s  |7,8|  on  the  basis  of 
Moriguchi's  theory  of  incompatible  strains  |9|.  This  method 
consists  of  assuming  an  incompatible  strain  distribution  in  the 
weldment  and  calculating  the  residual  stresses  on  the  basis  of 
the  theory  of  elasticity,  recognizing  that  this  incompatibility 
is  a  form  of  singularity  in  an  elastic  stress  field.  Analyses  using 
singular  points,  such  as  a  concentrated  load  or  a  center  of  di¬ 
lation,  for  solving  stress  problems  under  various  boundary 
conditions  are  well  established  and  can  thus  lie  very  useful  for 
analyzing  residual  stresses  under  similar  boundary  condi¬ 
tions. 

One  major  drawback  of  this  method,  however,  is  that  it  is 
very  difficult  to  estimate  the  distribution  of  the  incompatible 
strains.  In  some  representative  cases  the  distribution  has  been 
found  experimentally  1 10|.  On  that  basis,  estimates  can  then 
l>e  made  for  other  similar  situations. 

It  is  nevertheless  felt  that  since  the  advent  of  the  modern 
digital  computer  a  more  accurate  prediction  of  the  residual 
stress  distribution  can  be  made  by  simulating  the  phenomena 
that  occur  during  the  welding  process. 

Computer  simulation — To  accurately  simulate  the 
thermo-mechanical  behavior  of  a  weld,  or  for  this  matter  any 
phenomenon  that  involves  Ixith  thermal  and  stress  analysis,  one 
should  start  from  first  principles,  which  in  this  case  is  the  first 
law  of  thermodvnamics.  This  would  mean  that  one  would 
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have  to  solve  a  problem  containing  mechanical  and  thermal 
coupling,  a  fact  that  makes  the  analysis  extremely  complicated, 
if  not  impossible,  on  the  basis  of  the  present  state  of  knowledge, 
except  in  cases  where  thermoelastic  modeling  is  sufficient. 
Initial  attempts  at  formulating  and  investigating  the  more 
complex  problem  of  coupled  thermoplasticity  have  been  un¬ 
dertaken  by  Mroz  1 1 1,12],  but  the  whole  subject  area  is  still  at 
its  initial  stages  of  development. 

It  becomes  necessary,  therefore,  to  uncouple  the  thermal  and 
mechanical  parts  of  the  welding  problem  and  solve  each  one 
separately.  The  assumptions  required  for  this  uncoupling  have 
been  examined  by  Hibbitt  and  Marcal  1 13);  the  most  critical 
ones  are  the  neglect  of  dimensional  changes  and  the  neglect  of 
cross-coupling  between  thermal  and  mechanical  work. 

Based  on  these  assumptions  the  problem  can  be  solved  in  two 
steps.  First,  the  heat  flow  during  welding  is  analyzed.  The 
obtained  temporal  and  spatial  temperature  distributions  are 
then  used  as  one  of  the  inputs  for  the  subsequent  strain  and  stress 
analysis.  Various  analytical  and  numerical  methods  that  can 
be  used  in  each  step  are  examined  in  the  next  sections  of  this 
part  of  the  paper. 

Problem  characteristics — Figure  3  shows  the  physical  phe¬ 
nomenon  to  he  examined  A  welding  arc  is  traveling  with  a 
speed  v  between  two  plates,  causing  them  fo  coalesce  by  pro¬ 
viding  filler  metal. 

The  four  parts  of  a  weld  that  are  subjected  to  different 
thermal  histories  are  shown  in  Fig.  4.  Part  1  constitutes  the 
filler  metal  which  is  deposited  molten  and  later  solidifies  as 
cooling  begins.  Part  2  is  the  part  of  the  joint  that  melts  and  later 
resolidifies  during  cooling.  Both  these  parts  define  the  weld 
metal.  Part  3  is  the  heat-affected  zone  (HAZ),  defined  for  steel 
as  that  part  of  the  joint  in  which  the  maximum  temperature 
reached  is  above  the  Ai  but  below  the  solidus  temperature  7 
Finally.  Part  4  is  the  base  metal 

These  four  thermal  histories  and  the  related  changes  in 
physical  properties  are  shown  in  Fig  5  for  the  case  of  a  single¬ 
pass  weld  1 1 4 1. 

Analysis  of  heat  transfer  during  welding 

The  importance  of  accurately  predicting  the  tem|x*rature 
distribution  during  welding  has  lxx>n  recognized  for  many  years 
by  Ixith  scientists  and  engineers  working  with  welding  prob¬ 
lems.  This  importance  stems  from  the  fact  that  most  of  the 
phenomena  subsequently  encountered,  such  as  residual  stress**;, 
distortion,  and  metallurgical  changes,  have  their  origin  m  the 
uneven  temix-rature  distribution  and  the  fast  heating  and 
cxxiling  rates  that  occur  during  the  welding  o|X“ration. 

All  of  the  early  attempts  at  solving  the  problem  ol  heat  flow 
during  welding  were  analytical  in  nature  since  they  w  ere  |x-r- 
formed  Ixfore  the  advent  ol  the  computer.  As  a  consequence 
several  simplifying  assumptions  had  to  lx-  made  to  allow  for  the 
solution  of  the  highly  nonlinear  governing  partial  differential 
equation  and  the  accompany  ing  Ixiundarv  conditions  Starting 
in  the  mid-sixties,  however,  several  investigators  1  rum  around 
the  w  orld  have  used  the  computer  to  numerically  predict  the 
thermal  history  during  welding  with  much  greater  accuracy. 

As  will;  any  complex  engineering  problems,  the  choice  of 
w  hether  to  use  an  analy  tical  or  a  numerical  solution  has  to  lx- 
based  on  a  cost  versus  accuracy  tradeoff.  Analytical  solutions 
are  much  more  mex|x-nsive  and.  though  not  as  accurate  as  the 
numerical  ones,  provide  nevertheless  for  the  establishment  of 
the  general  laws  and  thus  facilitate  a  gixxl  understanding  of  the 
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phenomena  involved.  On  the  other  hand,  the  more  expensive 
numerical  solutions  are  necessary  whenever  accuracy  is  of 
paramount  importance,  as  for  example  when  a  metallurgical 
characterization  of  the  weld  metal  and  the  HAZ  is  needed,  or 
when  a  subsequent  stress  analysis  to  determine  the  transient 
strains  and  residual  stresses  is  required. 

Closed-form  analytical  solutions 

Fundamental  solutions — The  first  exact  analytical  solutions 
of  the  problem  of  heat  flow  during  welding  were  obtained  by 
Rosenthal  [  15— 17|  in  the  late  1930’s  and  early  1940's,  although 
a  particular  case  was  considered  independently  at  around  the 
same  time  by  Boulton  and  Lance  Martin  [18].  Rosenthal  solved 
the  conventional  heat  conduction  differential  equation  for 
constant  point,  line,  and  plane  heat  sources  moving  at  a  constant 
speed  with  respect  to  a  fixed  Cartesian  coordinate  system.  To 
facilitate  easier  handling  of  the  problem,  he  assumed  that 
welding  was  performed  over  a  sufficient  length  so  that  the 
temperature  distribution  around  the  heat  source  would  not 
change  if  viewed  from  a  coordinate  system  moving  with  the 
heat  source.  This  phenomenon  is  called  quasi-stationary  or 
quasi-steady  state.  Additional  assumptions  were  made  as  fol¬ 
lows: 


1.  The  physical  properties  of  the  conducting  medium  are 
constant. 

2.  The  heat  losses  through  the  surface  of  the  conducting 
medium  to  the  surrounding  atmosphere  are  neglected. 

3.  Heat  created  in  electric  welding  by  the  Joule  effect  is 
negligible. 

4.  The  phase  changes  and  the  accompanying  absorption 
or  release  of  latent  heat  in  the  conducting  medium  are  ne¬ 
glected. 

5  The  conducting  medium  is  infinitely  large  in  the  two- 
dimensional  case  (line  heat  source)  and  semi-infinitely  large 
in  the  three-dimensional  case  (point  source). 


Based  on  the  foregoing  assumptions  Rosenthal  developed  the 
following  exact  solutions  for  the  two-  and  three-dimensional 
cases  respectively  (see  Fig.  3): 
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(2) 


w  here 

II a  =  initial  temperature 
=  total  heat  input 
H  =  plate  thickness 
£  =  x  —  vt 

r  =  welding  (or  arc  travel)  S[x*ed 
t  -  tirne 

ix.y.z)  =  fixed  Cartesian  coordinate  system 
r  =(£-  +  !/-)'  ’ 

R  =  (£2  +  (/-  +  z-)'  2 
k  =  thermal  conductivity 
'/> A  =  k  i>r  =  k  =  thermal  diffusivity 
/i  =  density 
c  =  specific  heat 

K„{x)  =  modified  Bessel  function  of  sis  mid  kind  and  zero 
order 

For  the  case  of  thin  plates  (two-dimensional  solution) 
Rosenthal !  I(i|  and  other  investigators  have  suggested  that  heat 
losses  through  the  surface  to  the  surrounding  atmosphere  might 
have  to  fie  taken  into  account  bv  replacing  in  the  Bessel  function 
of  equation  ( I  I  the  factor  An  by 
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Flfl.5  Thermal  histories  and  related  changes  of  physical  properties 
of  the  four  weld  parts 


(An)2  + 


h\  +  h-z  1  2 
kH 


(3) 


where  hi  and  h2  are  the  heat-transfer  coefficients,  assumed 
constant,  at  the  top  and  bottom  of  the  plate,  respectively. 

Furthermore,  by  using  the  so-called  "method  of  images"  or 
“image  source  method”  [16],  one  can  get  solutions  for  the  cases 
of  large  but  finite  thickness  and/or  finite  breadth  plates.  For 
example,  a  three-dimensional  solution  for  laying  a  weld  bead 
on  the  top  of  a  finite-thickness  plate  with  adiabatic  boundary 
conditions  can  be  obtained  by  adding  an  infinite  series  to 
equation  (2),  yielding 
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H.I  =!£2  +  y-  +  (2nH  -z)-\'  2 

K„  =l£2  +  y2  +  (2  nil  +  z)2)1  2 

Following  these  initial  developments  many  investigators 
tested  the  validity  of  the  proposed  equations  experimentally . 
A  thorough  exjxisition  of  these  works  can  lx1  found  in  Myers  et 
al  |  !9|.  Comments  on  this  subject  are  reported  in  a  later  sei  - 

(ion 

Modifications  to  fundamental  solutions  Looking  at  data 
on  thermal  conductivity,  specific  heat,  and  density  one  finds 
that  all  these  parameters  are  highly  dependent  upon  tenqx'r- 
ature.  thus  making  the  constant- properties  assumption  of  the 
analytical  solutions  unrealistic,  especially  for  the  regions  close 
to  the  heat  source  where  the  material  exhibits  very  high  tem- 
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Table  1  Value*  of  arc  efficiency  for  various  processes 


Christensen  [23] 

Rykalin  [25| 

Tsai  [26) 

GMAW 

0.65  to  0.85 

mild  steel 

0.66  to  0.70 

0.80  to  0.90 

aluminum 

0.70  to  0.85 

SAW 

0.90  to  0.99 

0.90  to  0.99 

0.85  to  0.98 

SMAW 

mild  steel,  ac 

0.66  to  0.85 

0.65  to  0.85 

0.55  to  0.90 

GTAW 

mild  steel,  ac 

0.22  to  0.48 

0.20  to  0.50 

mild  steel,  dc 

0.36  to  0.46 

0.45  to  0.75 

aluminum,  ac 

0.21  to  0.43 

0.20  to  0.50 

GMAW  =  gas  metal  arc  welding. 
SAW  =  submerged  arc  welding. 
SMAW  =  shielded  metal  arc  welding. 
(’■TAW  =  gas  tungsten  arc  welding. 


peratures.  It  is  desirable  to  predict  the  high-temperature  region 
as  accurately  as  possible,  however,  since  it  is  this  region  that  is 
directly  related  to  the  size  of  the  plastic  zone  and  the  accom¬ 
panying  residual  stresses  and  distortion. 

Researchers  at  M.I.T.  (20-22]  adopted  the  iteration  method 
to  take  into  account  the  temperature  dependence  of  material 
properties.  The  fundamental  heat  source  solution  with  ma¬ 
terial  properties  at  some  temperature,  say  300°C  (572° F), 
provides  the  first  approximate  solution  at  a  particular  point. 
This  temperature  is  compared  with  the  initial  guess  and  if  the 
two  temperatures  disagree  by  more  than  0.5°C,  new  properties 
are  found  for  a  temperature  halfway  in  between.  The  >  new 
values  are  used  to  obtain  a  new  temperature  estimate.  The 
process  is  repeated  until  convergence  is  reached.  It  should  be 
pointed  out  that  although  this  iteration  method  generally  gives 
good  predictions  outside  the  fusion  zone,  there  is  no  guarantee 
that  it  will  converge  to  the  correct  solution  since  the  approxi¬ 
mation  used  may  not  satisfy  the  energy  conservation  law. 

In  addition,  the  conventional  point  heat  source  closed  from 
solution  fails  to  give  good  results  in  the  case  of  multipass 
welding.  This  is  due  to  the  fact  that  the  solution  is  based  on  the 
point  source  being  located  at  the  ten  surface  of  the  plates  being 
welded. 

To  accommodate  the  multipass  welding  case,  a  modification 
of  the  solution  was  made  (2|  enabling  one  to  locate  the  heat 
source  at  any  point  through  the  plate's  thickness.  It  thus  be¬ 
comes  possible  to  simulate  each  welding  pass  by  positioning  the 
(X)int  source  at  the  center  of  the  pass. 

The  basic  assumptions  of  the  conventional  solution  were  kept 
the  same.  Furthermore,  the  adiabatic  boundary  conditions  on 
the  top  and  Imttom  surfaces  of  the  plate  were  satisfied  by  using 
the  method  of  images.  The  obtained  solution  can  then  lx-  ex¬ 
pressed  by  the  following  equation: 


"l  «  +  ..“I  H„  +  R«  |i 

where 

Rn  —  \  £"  +  4-  (OTn  e  c )- 

R„  =  \  i2  ±  ij-  f  tOfl„  —  jP 

or,,  =  v  ofvTTTr 

ORn  =  v  Of^\ 

w  ith  F  and  G  being  the  distances  of  the  point  heat  source  from 
the  top  and  bottom  surfaces  of  the  plate,  respectively .  and  ( >7'„ 
and  ()B„  the  distances  of  the  point  heat  source  from  the  nth 


imaginary  ones  with  respect  to  the  top  and  bottom  surfaces, 
respectively.  All  other  quantities  have  been  previously  de¬ 
fined. 

Three-dimensional  finite  heat  source  model — Most  inves¬ 
tigators  agree  that  perhaps  the  most  critical  input  required  for 
the  welding  thermal  analysis  is  the  power,  Q,  that  enters  the 
plate  or  section  being  welded.  It  is  customary  to  express  this 
total  heat  input  by  the  formula 

Q  =  r)a  •  V  •  I  (6) 

where  V  and  /  are  the  arc  voltage  and  current,  respc  ctively, 
their  product  giving  the  electric  arc  power.  The  other  pa¬ 
rameter  in  the  equation  (6),  rja,  is  called  arc  efficiency;  it  rep¬ 
resents  the  ratio  of  the  power  introduced  by  the  arc  into  the 
metal  to  the  total  electric  arc  power.  In  other  words  it  provides 
a  semi-empirical  way  of  taking  into  account  the  various  heat 
losses  that  occur  through  electrode  tip  heating,  radiation  to  the 
surrounding  atmosphere,  metal  spatter,  etc. 

Arc  efficiency,  r?„,  is  heavily  dependent  on  the  welding 
process  used,  the  penetration  achieved,  the  shielding  gas  and 
many  other  factors  that  make  it  very  difficult  to  predict.  It  is 
therefore  estimated  either  through  experimental  measurements 
using  the  calorimetry  method  [23—25]  or  through  semiempirical 
correlations  with  the  weld  width  in  the  case  of  single-pass 
welding  1 26|.  Values  of  the  arc  efficiency  for  various  welding 
processes  as  proposed  by  several  investigators  are  given  in  Table 
1. 

Of  equal  importance  to  the  magnitude  of  the  total  heat  input 
is  its  distribution.  At  the  solutions  presented  earlier  a  point  or 
line  heat  source  was  assumed  As  Rykalin  [25 1  and  other  in¬ 
vestigators  report,  however,  a  more  realistic  approach  is  to  as¬ 
sume  a  Gaussian  radial  heat  flux  distribution  of  the  form 

q(r)  =  q,,-e~(:r2  (7) 

where 

</o  =  maximum  heat  flux  at  center  of  heat  spot 
=  CQ/tt,  |W/cm*| 

C  =  heat  flux  concentration  coefficient  =  3/rj,  [cm_2| 
r  =  radial  distance  from  center  of  heat  spot,  |cm| 

Th  =  radius  of  heat  spot,  [cm] 

Equation  (7)  is  valid  for  a  stationary  arc.  High-speed  cine¬ 
matography  reveals,  however,  that  during  welding,  when  the 
arc  is  moving,  the  arc  column  is  not  radially  symmetric  but 
rather  distorted  backwards.  This  observation  led  Tsai  |2f>|  to 
propose  the  following  equation  for  the  arc  heat  distribution 
instead  of  equation  (7): 

(/tr^l  =  q„- v,  i  (HI 

where  all  the  symlxils  have  Ixvn  previously  defined. 

Rased  on  the  preceding  discussion  it  is  ev  ident  that  in  order 
to  more  accurately  describe  the  tenqierature  distribution  and 
cr Milinc.  rates  in  the  region  close  to  the  weld  more  realistic  as¬ 
sumptions  should  be  employed.  A  three-dimensional  finite 
beat  source  model  for  sob  ing  the  governing  partial  differential 
equation  of  heat  transfer  was  thus  develop'd  [2  under  the 
follow  ing  assumptions 

1  (^uasi-slationary  state,  that  is.  steady -state  conditions 
ns  it  1 1  respect  to  a  coordinate  system  moving  with  the  heat 
source. 

2  The  heat  input  is  provided  l>\  a  i  loving  three-dimen¬ 
sional  skewed  normally  distributed  heat  source  moving  on  the 
surface  of  the  plate  and  given  by  equation  (8). 

3  The  thermal  conductivity  of  the  material,  k.  is  assumed 
to  lie  a  linear  function  of  temperature.  0.  given  In 

k\H)  =  k„-W  +  -yt«-«„:,  (9) 
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where  Ar0  is  the  value  of  the  thermal  conductivity  at  the  initial 
plate  temperature  and  y  the  proportionality  coefficient 

4.  The  thermal  diffusivity  of  the  material,  k,  is  assumed  to 
be  constant. 

5.  Convective  and  radiation  boundary  heat  losses  from  the 
plate's  surface  are  taken  into  account  through  a  constant  av¬ 
erage  '  effective"  heat-transfer  coefficient  h.  which  can  be 
different  for  the  top  and  bottom  surfaces  t  plate. 

6.  The  initial  temperature  of  the  plate,  f\,  in  be  different 
from  the  environmental  (ambient)  temperature,  6e,  to  allow 
for  preheating 

7.  Phase  transformation  and  Joule  heating  effects  can  be 
neglected 

The  final  result  obtained  from  this  model  is  given  in  Ap¬ 
pendix  1. 

Numerical  solutions 

Since  the  advent  of  the  electronic  digital  computer  several 
serious  efforts  have  been  made  to  numerically  solve  field,  and 
in  particular,  heat-transfer,  problems.  And  although  the  finite 
difference  method  had  initially  the  edge,  the  advantages  of  the 
finite-element  method,  especially  if  coupled  with  thermal  stress 
analysis,  are  more  and  more  recognized  today. 

Over  the  years  many  finite-element  programs  have  thus  been 
developed  that  are  capable  of  p<  rforming  heat-transfer  anal¬ 
yses.  Several  of  these  codes  can  take  various  nonlinearities  into 
account  in  a  more  <'r  less  sophisticated  manner.  One  of  the 
most  sophisticated  ones  is  ADINAT  (Automatic  Dynamic  In¬ 
cremental  Nonlinear  Analy  sis  for  Temperatures)  developed  by 
Bathe  and  co-workers  1 27,28)  over  a  period  of  years.  Some 
details  of  this  program  will  be  discussed  later.  It  suffices  here 
to  mention  that  ADINAT  can  take  into  account  temperature- 
dependent  material  properties  as  well  as  nonlinear  convection 
and  radiation  boundary  conditions. 

The  discussion  so  far  dealt  with  multipurpose  heat-transfer 
finite-element  method  (FF,M)  programs.  Concurrently, 
however,  several  investigators  concerned  with  the  welding 
probk  have  developed  similar  programs  fiibbit  and  Marcal 
[  Id]  developed  such  a  program  in  1973  which  was  later  used 
by  other  investigators  too  At  M  l  T  a  team  headed  by  Ma- 
subuchi  also  developed  similar  programs  in  the  early  1970's  |20, 
29,30!  Friedman  131  33 1  has  also  made  substantial  contri¬ 
butions  in  the  case  of  (IT  A  welding  working  at  the  Westing- 
house  F.leetric  Corporation's  Bettis  Atomic  Power  Labora¬ 
tory 

The  way  each  of  the  aforementioned  individuals  approached 
the  various  as|xvls  of  the  welding  problem  are  discussed  at  the 
appropriate  places  in  the  next  subsections,  where  some  details 
on  the  finite-element  formulation  of  the  welding  heat  transfer 
problem  are  presented 

Governing  equation  -  The  governing  incremental  isopara¬ 
metric  finite-element  equations  for  the  nonlinear  heat-transfer 
problem  can  !>c  written  in  matrix  form  as  follows  [27,28] 

.'K*  +■  'K  f  'K-  MG  ='•  +  '*  " 

+  '  ■  ±'Q'  (10) 

or 

'KAH  •  (I0n) 

where  'K  is  the  eff  ‘olive  eonduetivitx  matrix  at  tune  t  con¬ 
sisting  of  the  condiK  ti\  it v.  nonlinear  convection,  and  radiation 
matrices,  '  '  "  Q is  the  heat  flow  vector  including  the  effects 
of  surface  heat  flow  inputs,  internal  heat  generation  and  tern 
perature-dependeiil  heat  capacitv,  '  ‘  '  1  is  the  effective 

heat  flow  vector,  and  Ad 1 '  vs  the  increment  in  mxlal-point 
temperatures  in  iteration  i 


Tabt*  2  Temperature  dependence  ot  heat  convection  coefficient 


O'-0e[R\ 

/t[Btu/sec  in.2  R  \ 

Case  1 

Case  2 

Case  3 

0 

0 

0 

0 

100 

0.2  X  10~5 

0.2  X  10-5 

0.1  x  10-4 

500 

0.5  X  10~5 

0.1  X  10-4 

0.6  X  10~4 

1000 

0.1  X  10-" 

0.2  X  10-4 

0.1  X  10-3 

5000 

0.58  X  10"3 

0.58  X  10-" 

0.6  x  10-  > 

50000 

0,12  X  10-2 

0.12  X  10-2 

0.12  X  10-2 

Note:  0*  =  surface  temperature. 

8P  =  environmental  temperature. 


I  +  i<A!0(i)  —  I  +  (a) 

Equation  (10)  represents  the  heat  flow  equilibrium  at  time 
f-f-o-Af,  where  0<a<l  and  «  is  chosen  to  obtain  optimum 
stability  and  accuracy  in  the  solution.  Furthermore,  the  so¬ 
lution  using  these  equations  corresponds  to  the  modified 
Newton-Raphson  iteration  scheme. 

Boundary  conditions — Convection  and  radiation  boundary 
conditions  are  taken  into  account  by  including  the  matrices  ‘Kc 
and  ‘Kr  and  the  vectors  * +  <>A'  Q?  and  ,  +  i'%lQr  in  equation 
(10).  Additional  external  heat  now  input  on  tne  boundary  is 
specified  in  1  f  ,"it  Q  as  surface  heat  flow  input.  Prescribed 
temperature  conditions  can  also  be  specified. 

.  With  respect  to  the  welding  problem  these  boundary  con¬ 
ditions  can  be  stated  in  more  detail  as  follows: 

1 .  Convection  heat  losses  from  the  plates’  surfaces  can  be 
modeled  according  to  Newton’s  linear  law 

q'  =  h-(0e-0')  (12) 

where  h  is  the  temperature-dependent  convection  coefficient, 
tlP  the  environmental  temperature,  and  0s  the  surface  tem¬ 
perature  at  the  point  under  consideration.  There  is  some  dif¬ 
ficulty  in  estimating  the  temperature  dependence  of  h.  Al 
though  most  previous  investigators  considered  /i  to  be  constant, 
efforts  were  made  at  M.I.T.  to  rationally  estimate  it  on  the  basis 
of  past  semi-empirical  studies.  Table  2  provides  the  three 
different  sets  of  values  tested.  Discussion  on  the  obtained  re¬ 
sults  are  made  later. 

2.  Radiation  heat  losses  are  very  significant  in  the  vicinity 
of  the  weld  metal  because  of  the  large  difference  between  the 
surface  and  environmental  temperatures.  These  losses  are 
modeled  according  to  the  quartic  Stefan-Boltzman  law 

q'  =  <nA0}  ~  O'4)  (13) 

where 

a  -  Stefan-Boltzman  constant 

(  =  emissivity  of  surface 

A  =  shape  factor 

0r  =  sink  temperature 

0 *  =  surface  temperature 

In  the  M.I.T.  studies  the  shape  factor  was  taken  to  Ire  un  >  and 
the  emissivity  cot'fficient  0.8 

3  Temperatures  can  lx-  prescribed  at  any  point  and/or 
surface  of  a  weldment.  Such  a  case,  however,  seldom  arises  in 
wel..  .ig  analyses. 

4  The  heat  input  during  welding  is  modeled  using  the 
concept  of  arc  efficiency .  For  the  space  distribution  of  the  heat 
input,  a  consistent  formulation  should  be  adopted  assuming  a 
uniform  distribution  over  the  top  of  each  weld  bead  Finally, 
the  tmie  distribution  of  the  heat  input  is  modeled  in  such  a  way 
that  the  passing  of  the  arc  over  the  cross  section  examined  couid 
lx-  simulated  (a  linear  increase  as  the  arc  approaches,  uniform 
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as  the  are  travels  over  the  cross  section,  followed  by  a  linear 
decrease) 

Time  integration  schemes— A  family  of  one-step  methods 
1 34 1  is  considered  for  the  time  integration  using  the  parameter 
a.  The  scheme  is  fo>\  to  be  unconditionally  stable  for  <v  > 
[A  and  to  generally  give  better  solution  accuracy  when  <»  =  1 
(Euler  backward  method).  Finally  the  modified  Newton  it¬ 
eration  is  guaranteed  to  converge  provided  the  time  step  A t  is 
small  enough. 

Material  properties — One  of  the  significant  features  of  the 
finite-element  method  is  that  it  can  take  into  account  any 
nonlinear  dependence  of  the  material  physical  properties  with 
temperature,  something  that  is  very  important  in  the  welding 
analy  sis  Furthermore,  it  can  incorporate  the  latent  heat  of 
fusion  or  of  any  solid-state  material  transformation,  which 
cannot  be  furnished  bv  closed-form  analytical  solutions  1 13, 
:W|. 

Dimensionality  of  problem  —  At  first  glance  even  the  simpler 
welding  heat- transfer  analysis  looks  as  a  three-dimensional  one. 
It  the  assumptions  are  made,  however,  that  the*  welding  heat 
source  (arc)  is  moving  at  a  constant  speed  along  a  regular  path 
llor  example,  a  straight  line  in  a  planar  weld),  that  the  weld 
speed  is  sufficiently  high  relative  to  the  material’s  characteristic 
diffusion  rate,  and  that  end  effects  resulting  from  either  initi¬ 
ation  or  termination  of  the  heat  source  can  In-  neglected,  the 
three-dimensional  character  of  the  heat-transfer  problem  can 
lie  simplified  This  can  lie  achieved  by  analyzing  a  cross  section 
of  the  weldment  ol  unit  thickness  and  located  in  the  midlength 
region  of  the  weld  |2,'51  i. 

Molten-pool  modeling  The  heat-transfer  mechanism  in 
the  weld  metal,  when  molten,  is  extremely  complex  and  its 
phy  sics  are  not  well  understood  as  of  today  These  complexities 
arise  not  only  from  the  difficulty  involved  in  modeling  the 
welding  arc  heat  flux  correctly,  hut  also  from  the  liehaviorol 
the  convective  motion  of  the  molten  metal,  the  thermal  prop¬ 
erties  of  the  molten  metal  i  including  the  phase  transformations 
that  take  placcduring  melting  and  solidification),  the  electric 
heating  due  to  the  current  flow  in  the  I  use  metal,  the  Ixnindary 
conditions  for  heat  losses,  etc  Mam  investigators  have  rec¬ 
ognized  these  difficulties  and  have  tried  to  find  ways  to  ap¬ 
proximate  the  phenomena  involved 

A  look  at  the  available  literature  reveals  that  there  are  gen¬ 
erally  three  way  s  tor  handling  the  problem  The  first  one.  still 
at  its  dt  ■veloping  stages  tries  to  understand  and  subsequently 
mathematically  model  the  physical  phenomena  involved,  that 
is.  the  fluid  flow  .  the  convective  heat  transfer,  etc  '.'V>i  Since 
no  conclusive  general  results  are  available  vet.  this  approach 
cannot  be  generally  used 

In  the  second  wav  the  problem  is  divided  into  two  parts 
I'irst  the  shape  ol  the  molten  pool  is  semi -empiric  ally  deter 
mined,  then  the  heat  f  low  equations  are  solved  numericall'  m 
the  solid  metal  onlv  using  the  melting  isotherm  as  a  boundary 
condition  2ti  This  method  cannot  lx-  used,  however.  Ixvause 
the  temperature  distribution  in  the  weld  metal  has  to  be  cal¬ 
culated  too  d  a  stress  analysis  is  to  follow  imost  plastic  defor¬ 
mation  that  causes  the  tin  mat  ion  of  residual  stresses  takes  place 
in  this  region i 

I'lie  third  and  filial  approach  tries  to  simulate  the  convective 
heat  transfer  nice  linnism  in  the  molten  metal  bv  using  a  v  .din- 
tor  the  thermal  comlm  tiv  itv  of  the  molten  metal  an  order  ol 
magnitude  higher  than  that  of  the  material  at  the  solidus  tel u 
perature  This  approai  h.  although  strictly  not  physically 
correct .  was  used  in  the  M  I  I  investigation 

Finally  since  the  cases  analv  zed  involved  multipass  welding, 
it  is  necessary  to  I  mil  a  wav  to  model  t  fie  lav  ing  ol  the  various 
I  leads  during  the  welding  evde  Ibis  was  made  |xixsiblc  bv  the 
element  birth  and  death  <  a| labilities  ol  the  rode  used  In  oilier 
w  i  ids.  (he  program  \  I  >  I  N  \  I  is  <  i  pa  hie  ol  giv  mg  birth  to  a 


predetermined  number  of  elements  at  predetermined  time 
instances,  thus  enabling  one  to  model  the  laying  of  a  bead  by 
specifying  the  appearance  of  the  elements  representing  it  at  the 
time  it  physically  appears  |27|. 

Analysis  of  transient  strains  and  residual 
stresses  due  to  welding 

Using  the  temperature  distributions  predicted  on  the  basis 
of  the  techniques  described  in  the  previous  section,  one  can 
calculate  the  transient  strains,  transient  stresses  and  residual 
stresses  due  to  welding  since  the  problem  is  assumed  to  Ire  un¬ 
coupled  The  calculation  of  strains  and  stresses,  however,  poses 
a  much  more  formidable  problem  than  the  one  encountered 
in  the  heat-transfer  analysis,  making  the  use  of  numerical 
techniques  a  necessity.  These  difficulties  stem  from  the 
complicated  thermal-elastic-plastic  state  developed  in  and 
around  the  weld  metal  during  welding 

Two  general  techniques  have  been  developed  to  solve  the 
problem  One  is  a  simple  one-dimensional  analysis  and  the 
other  a  more  sophisticated  one  based  on  the  finite-element 
method.  Both  will  lx-  considered  here,  although  emphasis  will 
lx>  placed  on  the  latter 

One-dimensional  analysis 

The  one-dimensional  model  for  calculating  stresses  parallel 
to  the  weld  line  only  was  first  develo|>ed  in  1 904  by  Tall  |3fi|. 
Using  this  model,  Masubuchi.  Simons,  and  Monroe  |37|  devel¬ 
oped  a  computer  program  which  was  later  modified  and  im¬ 
proved  both  at  M  I  T.  and  elsewhere  1 20-22,38,39]. 

The  basic  assumption  inherent  in  the  one-dimensional  stress 
analy  sis  is  that  the  only  stress  present  is  the  one  parallel  to  Un¬ 
weld  line,  (i,,  and  that  this  stress  is  a  function  of  tin-  transverse 
distance  from  the  weld  centerline  only  As  a  consequence,  the 
equilibrium  conditions  are  not  satisfied.  Despite  this,  however, 
it  ap|x-ars  that  the  obtained  solutions  correlate  reasonably  well 
vv  it li  experimental  data  in  certain  cases,  as  will  be  further  ex¬ 
plored 

The  analy  sis  lends  basically  the  procedure  originally  promised 
bv  Tall  ].%|;  the  computer  e<xle  implementing  the  solution  with 
the  I;  esl  modifications  can  be  found  in  1 22 1 

The  algorithm  for  solv  ing  the  problem  is  based  on  the  method 
of  successive  elastic  solutions  as  pro|X>sed  by  Mendclson  |-l()j 
The  program  can  take  into  account  the  temperature  depen¬ 
dence  ol  all  material  properties,  any  tv  pe  of  strain  hardening 
and  can  solve  Ix-ad  on  plate,  edge,  and  butt  we  Ids  of  flat  plates 
with  finite  width  One  of  die  input  requirements,  the  tem¬ 
perature  distribution,  can  also  he  ealeulated  if  desired  In  the 
same  program  usinj'  the  line  heat  source  solution.  The  output 
at  each  time  step  consists  of  tin-  temperature,  total  strain,  me 
chanu  al  strain,  plastic*  stram.  and  stress  at  each  of  the  prede¬ 
termined  |>oints  located  at  various  transverse-  distances  from  the 
vs  eld  centerline 

I’rex  ions  applications  i  »l  the  program  have  show  n  that  it  i»ives 
Hood  Jesuits  m  the  ease  of  thin  plates  This  happens  because 
in  thin  plates  all  stresses,  except  a,  are  vein  small  sometimes 
an  order  of  maumtude  smallei  than  ct,  l  unher  discussion  is 
delaved  until  the  next  section  when  experimentally  obtained 
i  exults  are  compared  with  analv  heal  solutions 

l  inallv.  it  should  In*  mentioned  that  tlieone-dimension.il 
analv  s|s  ran  ho  used  with  pi  ope  r  minor  modifications  for  the 
analv  sis  of  simple  st  met  in  a!  forms  othei  than  plates  \n  ex¬ 
ample  is  the  welding  analysis  of  built-up  I  -beams  ' 20.22! 

Finite-element  analysis 

The  complex  hehav  ior  of  ,i  weldment  and  in  particular  the 
hnjlilv  nonlmeai  material  response  and  the  material  loading 
and  unloading  lh.it  «»e< m  in  the  multidimensional  stress  space 
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can  be  handled  more  accurately  using  numerical  techniques 
such  as  the  finite-element  method. 

One  of  the  first  applications  of  FEM  to  weld  problems  was 
presented  by  Hibbitt  and  Marcal  (13|,  who  considered  a 
thermo-mechanical  model  for  the  welding  and  subsequent 
loading  of  a  fabricated  structure.  Their  model  simulates  GMA 
welding  processes  and  accounts  for  temperature-dependent 
material  properties.  Friedman  1 3 1  —33 )  also  developed  fi¬ 
nite-element  analysis  procedures  for  calculating  stresses  and 
distortions  in  longitudinal  butt  welds.  These  procedures  are 
applicable  to  planar  or  axisymmetric  welds.  Rybicki  and  co- 
workers  (41 1  have  developed  similar  procedures.  At  M.I.T.  a 
team  headed  by  Masubuchi  has  also  developed  two-dimensional 
finite-element  programs  capable  of  performing  plane-strain 
and  plane-stress  analyses  ( 1,20, 29,30 1. 

Japanese  investigators  have  also  been  very  active  in  the  field. 
Ueda  and  co-workers  [14],  Satoh  et  al  [42|,  and  more  recently 
Fujita  and  Nonioto  (43|  have  all  developed  models  to  calculate 
transient  strains  and  residual  stresses  due  to  welding  based 
primarily  on  the  initial  strain  method. 

A  recent  report  [ 44 1  contains  an  annotated  bibliography  of 
research  efforts  in  the  area  that  have  been  published  since 
1977. 

In  this  subsection  some  details  on  the  finite-element  formu¬ 
lation  of  the  welding  stress  problem  are  presented.  Although 
the  discussion  pertains  to  the  multipurpose  nonlinear  stress  fi¬ 
nite-element  program  ADINA  [45,46],  it  is  felt  that  what  is 
discussed  can  be  equally  well  applied  to  any  similar  code. 

Finite-element  formulation — Some  basic  considerations 
regarding  the  thermo-elastic-plastic  and  creep  constitutive 
model  used  in  conjunction  with  ADINA  are  discussed  in  the 
following 

The  governu.g  incremental  finite-element  equations  for  the 
problem  can  be  written  as  |45| 

'KALI1*'  =  (14) 

where  K  is  the  tangent  stiffness  matrix  corresponding  to  time 
r;  '  •  is  the  nodal-point  external  force  vector  at  time  t  +  At: 

1 1  Atf  ii  —  1 1  jsa  vector  of  nodal-point  forces  that  are  equivalent, 
in  the  virtual  work  sense,  to  the  internal  element  stresses  at  time 
t  +  At  and  iteration  i  -  1 

(15) 

and  AV 11 '  is  the  increment  in  nodal-point  displacement  in  it¬ 
eration  i 

,  .  ,  =  / ,  Aip.-n  +  AD1"  (16) 

The  solution  using  equation  ( 14)  corresponds  to  the  modified 
Newton -Raphson  iteration  procedure  which  is  helpful  in  im¬ 
proving  the  solution  accuracy  and  in  many  cases  in  preventing 
the  development  of  numerical  instabilities.  The  convergence 
of  the  iteration  can  be  accelerated  using  the  Aitken  method  or. 
in  complex  material  nonlinear  cases  (like  the  welding  problem), 
improved  using  the  BFGS  (Bros  den- Fleteher-Goldfarb- 
Shanno)  matrix  updating  method  '  461 

In  the  thermo-elastic  -plastic  and  creep  model,  and  assuming 
infinitesimal  strains,  the  total  strain  at  time  r.  is  assumed 
to  lie  given  by 


w  here 

e,;  =  elastic  straii 
=  plastic  strain 
(’[,  =  creep  strain 
eJi  =  thermal  strain 


so  that  at  any  time  r  during  the  response  the  stress  is  given  by 
the  constitutive  law  for  an  isotropic  thermo-elastic  material 

rti  =  ^Cfir^ers  -  'e?  -  'et  -  ’ e ™)  (18) 

with  TCfjrs  denoting  a  component  of  the  elastic  constitutive 
tensor. 

The  thermal  strains  are 

Te™  =  ram(r6  -  0H)5rs  (19) 

where  Tam  is  the  average  thermal  expansion  coefficient,  0R  is 
the  reference  temperature,  and  5rs  the  Kronecker  delta.  This 
term  can  be  modified  to  include,  in  addition  to  the  thermal 
strain,  the  strains  that  are  induced  from  the  solid  phase  trans¬ 
formations  occurring  during  the  heating  and  cooling  stages  of 
the  thermal  history  [2,47 1. 

The  creep  strains  can  be  determined  by  any  of  a  number  of 
different  approaches  proposed  in  the  literature.  Given  that 
in  the  welding  problem  the  time  intervals  at  high  temperature 
are  short,  however,  creep  can  be  neglected. 

For  the  plastic  strains,  Tep,  the  situation  is  more  complicated. 
Although  the  classical  theory  of  isothermal  plasticity  is  a  well- 
tested  one,  extension  of  the  theory  to  nonisothermal  cases  is 
difficult  to  substantiate  experimentally.  Several  investigators 
have  proposed  modifications  but  very  few  experiments  have 
been  performed.  Relatively  good  agreement  between  theory 
and  experiments  has  generally  been  reported  but  for  temper¬ 
atures  up  to  about  538°C  (1000°F)  only.  During  welding, 
though,  the  temperatures  rise  to  above  the  A  \  temperature 
|654°C  (1210°F)  for  steel]  inside  the  HAZ-base  metal  boundary 
and  above  the  liquidus  temperature  in  the  weld  metal.  For 
lack  of  any  alternative,  however,  the  same  nonisothermal  theory 
of  plasticity  can  be  used  throughout  the  temperature  range 
encountered  in  welding  problems. 

The  general  form  of  the  yield  or  loading  function  for  mul- 
tiaxial  stress  conditions  ; 

F(’(rtj,’otij,’(ry)  (20) 

where  r«.,  and  r(T,  are  functions  of  the  history  of  plastic  de¬ 
formation  and  temperature.  For  elastic  behavior,  rF  <  0.  and 
for  plastic  behavior  ’F  =0  As  a  consequence  of  Drucker’s 
postulate  for  stable  plastic  materials.  rF  defines  a  convex  surface 
in  the  stress-temperature  space.  It  is  also  assumed  that  the 
isothermal  normality  condition  remains  valid,  so  that 


where  '  \  is  a  positive  scalar  The  selection  of  a  hardening  rule 
is  also  required  for  the  calculation  of  7\.  In  ADINA  either 
isotropic  or  kinematic  hardening  can  be  assumed.  Because 
cyclic  plasticity  is  expected  in  the  welding  problem,  the  kine¬ 
matic  hardening  mechanism,  thought  by  many  to  better  model 
the  phenomena  involved,  should  be  chosen.  The  assumptions 
involved  in  this  mechanism  are  that  the  size  of  the  yield  surface 
dc|>ends  on  the  temperature  only,  whereas  the  translation  rate 
of  the  yield  surface  in  the  stress  space  depends  on  the  plastic 
strain  rate. 

Further  details  on  the  numerical  asjrects  of  the  problem  can 
lx-  found  in  |46|. 

Geometry  and  boundary  conditions  A  cross  section  of  the 
weldment  in  its  midlength  can  lx-  used  to  calculate  the  transient 
strains,  transient  stresses,  and  residual  stresses  due  to  welding 
This  is  rationalized  by  the  fact  that  for  relatively  long  plates  the 
maximum  stresses  are  developed  in  this  region.  Furthermore, 
the  plane-strain  assumption  can  lx-  used  (that  is.  all  plane  sec¬ 
tions  normal  to  the  weld  line  remain  plane  during  the  entire 
welding  process) 

The  boundary  conditions  to  lx-  list'd  in  the  analy  sis  should  be 
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EXPERIMENT  No  3  (25  mm  FROM  WELD  LINE) 
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Fig.  6  Temperature  distribution — analysis  versus  experiment  (bead-on-plate  weld,  line  heat 

source) 


such  so  as  to  allow  free  expansion  of  the  weldment  in  the 
transverse  direction  as  well  as  trending.  At  the  same  time  the 
structure  should  be  properly  restrained  to  eliminate  all  possible 
modes  of  rigid-body  motion;  otherwise  the  stiffness  matrix  will 
not  be  positive  definite. 

The  foregoing  holds  for  a  relatively  thick  weldment  In  the 
case  of  a  thin  plate,  a  plane-stress  analysis  can  alternatively  la- 
performed  over  the  top  surface  of  the  plate  (assuming  negligible 
through-thickness  stress  gradients) 

Material  properties—  The  temperature  dependence  of 
Young's  modulus  E,  Poisson  s  ratio.  /'.  virgin  yield  stress,  crvll. 
and  strain  hardening  modulus,  /vj ,  are  all  required  input  for 
the  thermo-elastic-plastic  model  used 

A  word  of  caution  is  appropriate  al  this  point  as  tar  as  the 
material  properties  alxive  the  iiquidus  tenijx-rature  are  con¬ 
cerned.  The  material  does  not  have  am  strength  when  molten 
since  all  its  mechanical  pro|x-rties  are  zero  But  due  to  nu¬ 
merical  considerations,  zero  properties  cannot  generally  be 
entered  as  input  to  the  program  Hence  to  avoid  any  insta¬ 
bilities  (or  even  divisions  by  zero!  \er\  small  values  for  /•..  a,  u. 
anti  Er  should  lx-  used  above  the  Iiquidus  temperature 

Another  point  is  the  accumulation  of  plastic  -.trains  in  the 
regions  that  become  molten  during  the  welding  cycle.  When 
the  temperature  reaches  the  Iiquidus  these  plastic  strains  are 
physically  relieved,  starting  to  accumulate  again  when  the 
metal  solidifies.  The  presence  ol  nonzero  material  pro|x-rties 
above  the  Iiquidus.  however  would  cause  the  plastic  strains  not 
only  to  continue  accumulating  but  also  to  reach  artificiallv  high 
values  owing  to  the  very  low  magnitude  of  the  mechanical 
properties  It  was  therefore  necessarv  to  moddv  A1)I\A  !  v 
imposing  a  total  relief  of  plastic  strains  when  the  material 
melts. 

Solution  strategy  Oncol  the  most  unpoitant  decisions  an 
analyst  has  to  make  when  [K-rformmg  a  nonlinear  incremental 
stress  analysis  is  the  solution  strategy  to  lx-  followed.  Ix-eaiise 
the  accuracy  and  the  convergence  characterislu  s  of  the  solution 
depend  very  much  on  it  This  is  especiallv  true  lor  complex 


situations  involving  highly  nonlinear  material  behavior  like  the 
one  encountered  in  the  welding  problem.  It  is  therefore  nec¬ 
essary  to  perform  several  numerical  experiments  prior  to  a  full 
analysis  of  the  problem. 

Experimental  verification  of  models 

In  the  previous  two  sections  various  analy  tical  and  numerical 
methixls  for  the  prediction  of  temperature,  transient  strain,  and 
residual  stress  distributions  due  to  welding  were  descrilx-d.  It 
w  as  suggested  that  these  methods  could  be  very  useful  to  de¬ 
signers  and  fabricators  for  the  prediction  of  residual  stresses 
Before  the  methixls  can  lx-  widely  used,  however,  the  level  of 
confidence  in  their  predictive  capabilities  should  lx-  well  er 
tablished  This  can  lx-  done  by  comparing  their  predictions 
w  it  1 1  results  obtained  from  systematic  series  of  experiments  that 
cover  different  cases  encountered  in  practice.  Through  these 
comparisons  the  analy  tical  and  numerical  models  can  lx-  tested 
and  mixlified  where  necessary 

During  the  past  ten  or  so  soars  efforts  towards  this  goal  have 
Ix-en  continuously  made  by  several  M  I  T  investigators.  Their 
results  have  Ix-en  published  in  a  numlx-r  of  theses,  re|xirts,  and 
pa|x-rs  12,4.5.20  22.2fi.29.30.37.38.47j.  In  this  section  some 
ol  the  most  recent  investigations  will  lx-  rojxirted  At  the  same 
time  the  capabilities  and  limitations  of  the  several  methrxlol- 
ogies  will  lx-  discussed 

Temperature  distributions 

l  ine  heat  source  To  investigate  the  predictive  capabilities 
ol  the  line  heal  source  solution,  ev|x-riments  performed  on 
9  .95  nun  thick  i 1  y  in  l  9091  TO  aluminum  alloy  plates  -4S  were 
analv  zed  Head-on-plate  welds  using  tlx-  (Al  AW  process  wore 
performed  with  an  arc  voltage  ol  19  V.  arc  current  of  2-4(1  A. 
and  w elding  speed  o|  I  s  (i  mill  sec  it)  54  in  seel  I  he  plates 
were  79  2  cm  lot)  in  )  long  and  45  7  cm  i  1  S  in  I  vv  ide  Figure 
9  shows  how  the  experimentally  obtained  data  compare  with 
the  analv  t  leal  predictions  22  at  two  points  located  25  and  79 
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mm  (1  and  3  in.)  away  from  the  weld  centerline  and  on  the 
midsection  of  the  plate  (quasi-stationary  state).  It  can  be  seen 
that  the  correlation  is  very  favorable,  Note  that  this  would  not 
be  the  case  if  locations  closer  to  the  weld  centerline  were  con¬ 
sidered  because  of  the  singularity  (heat  source)  present  on  this 
line. 

Similar  good  predictions  have  been  obtained  when  plates  up 
to  about  12.7  mm  (V2  in.)  in  thickness  and  bead-on-plate  or  butt 
welded  using  conventional  processes  were  analyzed.  In  ad¬ 
dition,  analyses  of  plates  up  to  25.4  mm  (1  in.)  thick  and  welded 
using  high-energy-density  welding  processes  (electron  beam 
and  laser)  seem  to  accurately  predict  temperature  distributions 
(49). 

A  common  characteristic  of  all  the  cases  for  which  the  line 
heat  source  solution  gives  reasonably  accurate  predictions  at 
points  some  distance  away  from  the  weld  centerline  is  the 
presence  of  negligible  temperature  gradients  in  the  plate's 
thickness  direction  (two-dimensional  heat  transfer).  This 
characteristic  can  therefore  be  used  as  a  criterion  for  when  such 
a  solution  can  be  utilized. 

Finite  heat  source — As  previously  mentioned,  the  line  heat 
source  solution  does  not  give  accurate  results  at  points  close  to 
the  weld  centerline,  even  for  thin  plates,  the  reason  being  that 
the  heat  source  introduces  a  singularity  which  does  not  model 
the  phenomena  in  this  area  properly  This  difficulty,  however, 
can  be  overcome  by  the  three-dimensional  point  heat  source 
model  discussed  in  an  earlier  section. 

To  test  this  proposition,  experimental  results  previously  ob¬ 
tained  j 26 1  when  welding  a  3.2-mm-thick  (Vs  in. )  low-carbon 
steel  plate  were  analyzed.  The  welding  speed  was  3.80 
mm/sec  (0. 15  in. /sec)  and  the  arc  power  5000  J/sec.  Figure 
7  shows  the  comparison  between  experiments  and  analysis  at 
two  points  located  10.2  and  17.8  mm  (0.4  and  0.7  in.)  away  from 
the  weld  centerline  As  postulated,  the  finite  heat  source  model 
gives  good  results  in  this  thin-plate  case. 

Point  heat  source — The  point  heat  source  solution  is  more 
appropriate  when  relatively  thick  welded  plates,  say  more  than 
12  7  mm  (Vi  in),  are  to  be  analyzed.  This  is  the  case  because 
in  these  instances  the  heat-transfer  problem  is  a  three-dimen¬ 
sional  one  due  to  the  presence  of  temperature  gradients  in  the 
thickness  direction  too.  When  joining  thicker  plates,  however, 
more  than  one  welding  pass  is  required,  necessitating  the  use 
of  the  modifications  deseril>ed  at  an  earlier  section  (that  is,  the 
(Kisitioning  of  the  |>oint  source  at  the  center  of  each  pass). 
Moreover,  one  should  note  that  as  in  the  line  heat  source  case 
the  temperature  predictions  close  to  the  weld  centerline  will 
not  he  verv  accurate  due  to  the  singularity  induced  by  the 
presence  of  the  |xiint  source  To  remedy  the  situation  some¬ 
what  one  should  limit  the  maximum  lem[K*rature  calculated 
to  a  fixed  value,  such  as  the  liquidus  temperature;  otherwise, 
artificially  high  values  will  Ire  obtained. 

To  demonstrate  the  aforementioned  points  the  first  pass  of 
a  multipass  welding  experiment  on  a  25.4-mm-thiek  (1  in.) 
111-130  plate  welded  using  the  GMAW  process  was  analyzed. 
Details  of  the  experiments  can  Ire  found  in  |2,4,49|.  Figure  8 
presents  some  of  the  results  obtained  Shown  are  the  cases  of 
constant  material  pro|rerties  w  ith  the  point  heat  source  located 
on  the  top  surface  of  the  plate  1  conventional  solution)  and  at  a 
distance  1 1  1  mm  it)  4375  in  I  from  the  top  (simulating  the  first 
pass  of  tin  actual  welding  experiment)  and  the  case  of  variable 
properties  (iterative  solution) 

The  results  show  <  learlv  the  overestimation  of  temperatures 
if  the  conventional  point  heat  source  is  applied  even  if  the 
temperature  variation  ol  material  properties  is  considered 
Ibis  overestimatiou  is  even  more  pronounced  in  the  high- 
temperaturc  region  close  to  the  weld  centerline  Finally, 
comparison  of  the  modified  solutions  predictions  with  experi- 
mentallv  obtained  data  revealed  good  correlation  { 2 j 


sw 


Fig.  7  Comparison  of  finite  heat  source  solution  with  experimental 
data  for  a  3.2-mm-thick  welded  plate 

Finite-element  solution — To  test  the  predictive  capabilities 
of  the  more  accurate  nonlinear  finite  element  heat-transfer 
program  ADINAT,  a  series  of  experiments  involving  the  mul¬ 
tipass  GMA  welding  of  25.4-mm-thick  (1  in.)  HY-130  plates  was 
analyzed. 

Figure  9  shows  the  test  plate  arrangement  for  these  experi¬ 
ments.  The  weld  joint  consisted  of  a  double-V  groove  with  a 


Fig.  8  Investigation  of  point  heat  source  solutions  ( 12.7  mm  away  from 
weld  centerline) 
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Fig.  9  Test  plate  arrangement 


60-deg  included  angle  in  accordance  with  U  S.  Navy  specifi¬ 
cations.  The  test  plate  support  arrangement  consisted  of 
knife-edge  supports  located  through  the  whole  specimen  length 
at  127  and  203  mm  (5  and  8  in.)  from  the  weld  centerline  on 
either  side;  this  way  the  experiments  were  performed  with  the 
specimens  completely  unrestrained.  Twelve  passes  were 
needed  to  complete  the  welding  with  an  arc  voltage  of  25  V, 
a  travel  speed  of  5.1  mm/sec  (0.2  in. /sec)  and  a  heat  input  of 
14.6  kj/cm  (37  kj/ in. ).  Temperatures  and  transient  strains 
were  measured  throughout  the  welding  operation  using 
Chromel/Alumel  adhesive  bonded  thermocouples  and  90-deg 
Rosette  electric  resistance  strain  gages,  respectively,  located  at 
preselected  distances  from  the  weld  centerline  (the  instru¬ 
mentation  area  is  shown  in  Fig.  9). 

A  cross  section  of  the  plate  at  its  midlength  was  analyzed. 
Figure  10  shows  the  finite-element  mesh  used  for  the  analysis. 
Four-  to  six-node  isoparametric  elements  were  used,  special  care 
was  taken  not  to  include  any  triangular  elements  in  the  mesh. 
A  total  of  77  nodes  and  47  elements  was  utilized 

Table  3  summarizes  the  four  analyses  performed  Two 
parameters  were  varied,  the  arc  efficiency,  rj„,  and  the  tem¬ 
perature  variatioii  of  the  convection  coefficient,  h,  given  the 
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Fig.  10  Finite-element  mesh  used  in  the  analysis  of  the  welding 
problem 


fact  that  some  uncertainty  exists  regarding  their  true  values. 
All  other  variables  were  kept  constant  in  the  four  cases. 

Cases  A2,  A3,  and  A4  are  compared  with  experimentally 
obtained  results  in  Fig.  1 1  for  a  point  on  the  plate’s  top  surface 
12.7  mm  (0.5  in.)  away  from  the  weld  centerline  and  for  the 
first  welding  pass.  Considering  the  same  heat  input  (rja  = 
0.60),  little  difference  is  found  between  Cases  A2  and  A3.  This 
is  due  to  the  fact  that  although  higher  values  for  the  convection 
coefficient  were  chosen  in  Analysis  A3,  these  values  were  not 
high  enough  to  significantly  alter  the  heat  losses  from  the  top 
and  bottom  surfaces  of  the  plate  and  consequently  the  tem¬ 
perature  distribution.  A  substantial  increase  in  the  convection 
coefficient  was  therefore  chosen  for  Case  A4;  at  the  same  time, 
however,  an  increase  in  the  arc  efficiency  was  made  to  partially 
compensate  for  the  higher  h  and  thus  to  obtain  a  good  estimate 
of  the  maximum  temperature  reached.  As  seen  in  Fig.  1 1,  the 
combination  of  values  used  in  this  latter  analysis  succeeded  in 
bringing  the  cooling  rate  much  closer  to  the  experimental  one. 
At  the  same  time  a  4  percent  overprediction  of  the  maximum 
temperature  is  observed.  The  difference,  however,  was  very 
small  so  that  any  further  analysis  was  not  felt  necessary. 

Similar  results  were  found  at  other  points  of  the  plate  and  for 
the  consecutive  welding  passes.  Case  A4  always  came  closest 
at  matching  the  experimentally  obtained  results. 

Transient  strain  and  residual  stress  distributions 

One-dimensional  analysis — The  same  experiments  as  the 
ones  mentioned  in  the  line  heat  source  solution  were  analyzed 
to  test  the  predictive  capabilities  of  the  one-dimensional  com¬ 
puter  program.  Figure  12  shows  the  longitudinal  strain  history 
at  a  point  25.4  mm  (1  in.)  away  from  the  weld  centerline. 
Presented  are  experimental  data  and  predictions  based  on  the 
old  |21 1  and  new  |22)  versions  of  the  program  (the  most  im¬ 
portant  modification  is  a  change  in  the  numerical  integration 
scheme).  It  can  be  seen  that  the  analytical  results  have  a  good 
correlation  with  the  experimentally  obtained  ones.  Similar 
conclusions  can  also  lx*  drawn  with  regard  to  the  residual  stress 
distribution  [22|. 

To  investigate  the  limits  of  the  predictive  capabilities  of  the 
one-dimensional  model,  experiments  on  25.4-mm-thick  (1  in.) 
plates  welded  using  the  electron  beam  (EB)  and  the  CM  AW 
processes  were  analyzed.  The  EB  specimens  had  the  same 
geometric  configuration  as  the  CM  A  welded  ones  (described 
previously  w  hen  the  finite-element  heat-transfer  analysis  was 
discussed),  the  only  difference  lx*ing  the  joint  shape  which  was 
square  butt  in  the  cast*  of  EB  welding. 

Figure  13  shows  the  comparison  between  experiment  and 
analysis  of  the  longitudinal  transient  strain  history,  e,.  in  the 
EB  specimen  at  two  [joints  17  and  25.4  mm  (0.67  and  1  in.) 
away  from  the  weld  centerline.  Despite  the  thickness  of  the 
specimens  1 25.4  mm  ( 1  in.)|  the  results  show  a  remarkably  good 
correlation  of  analysis  and  experiment.  This  can  lx*  attributed 
to  the  physics  of  the  EB  welding  process  which  result  in  a  rel¬ 
atively  uniform  through-t hickness  temperature  distribution  and 
a  narrow  heat-affected  zone.  As  a  consequence,  stresses  ol(, 
<7j„,  and  o.  are  small  compared  with  <t,.  a  fact  substantiated 


Table  3  Summary  o(  analyses  performed 
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"  rj„  is  the  arc  efficiency  utilized. 
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Fig.  1 1  Comparison  of  finite-element  analysis  results  with  experimental  data  (weld  pass 

No.  1) 


by  the  experimental  data.  The  assumption  of  the  one-di¬ 
mensional  theory  that  <t,  is  the  only  stress  present  is,  therefore, 
approximately  valid.  It  should  be  emphasized  once  more, 
however,  that  these  results  apply  for  the  midlength  of  a  rela¬ 
tively  long  plate  | longer  than  457.2  mm  (18  in.)),  where  the 
maximum  possible  stresses  have  been  developed  and  where  the 
end  effects  can  be  neglected. 

The  comparison  between  experimental  data  and  analytical 
results  for  the  longitudinal  transient  strain  history  in  the  GMA 
welded  specimens  at  two  points  25.4  and  57.2  mm  (1.0  and  2.25 
in  )  away  from  the  weld  centerline  is  shown  in  Fig.  14.  Here 
the  story  is  completely  different.  The  correlation  is  not  good, 
especially  for  the  closest  point;  at  around  45  sec  from  the 


commencement  of  welding,  for  example,  the  experimental  data 
show  a  maximum  positive  longitudinal  strain  of  about  0.9, 
whereas  the  analysis  predicts  a  negative  longitudinal  strain  0.4 
in  magnitude.  It  is  thus  apparent  that  the  assumptions  involved 
in  the  one-dimensional  analysis  break  down  in  the  case  of 
multipass  GMA  welds  of  thick  plates.  Experimental  data  ac¬ 
tually  confirm  this  by  showing  that  the  transverse  strains  are 
of  the  same  order  of  magnitude  as  the  longitudinal  ones. 

In  conclusion,  then,  the  one-dimensional  program  seems  to 
be  appropriate  for  analyzing  thin  plates  welded  by  conventional 
recesses  or  medium-thickness  plates  welded  by  one  of  the 
igh-energy-density  processes. 

Finite -element  analysis — The  same  experiments  as  the  ones 


/\ 


12  Longitudinal  strain  history — analysis  versus  experiment  (bead-on-plate  weld; 
one-dimensional  program) 
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Fig.  13  Longitudinal  transient  strain  histories  for  EB  welded  specimen 
(experiment  and  analysis) 

analyzed  for  the  finite-element  heat-transfer  problem  were  also 
used  in  the  stress  analysis.  The  finite-element  mesh  used  was 
also  the  same  (see  Fig.  10);  now  the  appropriate  nodal  restraints 
were  also  imposed. 

Comparison  of  the  experimentally  measured  transverse 
transient  strain  history  at  a  point  located  at  the  top  surface  of 
the  plate  is  made  with  the  numerically  obtained  results  in  Fig. 
15.  The  correlation  is  relatively  good  if  one  takes  into  account 
the  various  assumptions  involved  in  modeling  the  complex 
welding  problem.  A  delay  in  the  transition  from  tensile  to 
compressive  and  again  from  compressive  to  tensile  strains  is 
observed  in  the  analysis.  The  same  delay  as  far  as  the  occur¬ 
rence  of  the  maximum  strain  is  also  exhibited.  It  is  believed 
that  this  phenomenon  is  primarily  due  to  the  relative  coarseness 
of  the  finite-element  mesh,  and  the  complex  loading  history 
present  in  the  welding  problem 

Figure  16  shows  calculated  transient  longitudinal  stress  dis¬ 
tributions.  Compressive  stresses  exist  in  the  weld  metal  prior 
to  melting  When  the  metal  is  in  its  molten  stage,  negligible 
compressive  or  tensile  stresses  were  calculated.  As  cooling 
commences,  tensile  stresses  start  appearing  in  the  weld  metal. 
These  stresses  then  build  up  to  the  residual  stress  pattern  when 
ambient  temperature  is  reached.  For  self-equilibrating  pur¬ 
poses,  compressive  stresses  exist  in  areas  removed  from  the  weld 
centerline.  Note  in  Fig.  16  the  effect  of  phase  transformation, 
causing  a  sudden  decrease  in  the  tensile  stresses. 

Summarizing,  it  can  be  said  that  the  finite-element  model 
developed  based  on  AD1NA  modified  for  phase  transformation 
effects  captures  most  of  the  important  aspects  of  the  welding 


Fig.  14  Longitudinal  transient  strain  histories  for  GMA  welded  spec¬ 
imens  (experiment  and  analysis) 

stress  analysis.  Although  some  discrepancies  have  been  ob¬ 
served  between  the  obtained  results  and  experimental  data, 
considerations  of  cost  and  the  assumptions  made  do  not  allow 
much  more  sophisticated  analyses  at  present.  It  is  expected, 
however,  that  the  dramatically  decreasing  computer  costs,  to¬ 
gether  with  developments  in  the  areas  of  coupled  thermo¬ 
plasticity,  will  enable  investigators  to  perform  more  accurate 
analyses  in  the  future. 


Part  2 

Analysis  of  consequences  of  residual  stresses 

This  part  of  the  paper  addresses  the  problem  of  analyzing  the 
consequences  of  residual  stresses  on  the  service  behavior  of 
welded  structures.  A  brief  introduction  on  the  subject  is  first 
iven,  mostly  qualitative  in  nature,  followed  by  a  discussion  on 
ow  the  various  problems  in  the  area  can  be  analyzed.  Because 
of  the  vastness  of  the  field,  however,  and  due  to  length  limita¬ 
tions,  one  particular  case  only  has  been  singled  out  and  analyzed 
in  some  detail.  It  consists  of  an  investigation  on  the  fracture 
characteristics  of  weldments  made  of  high-strength 
quenched-and-tempered  steels. 

Finally,  and  in  accordance  with  the  general  philosophy  of 
always  testing  the  predictive  capabilities  of  any  analysis,  the 
analytical  results  are  compared  with  data  obtained  during  a 
series  of  experiments. 


Fig.  15  Transverse  transient  strains  on  plate's  top  surface  (experiment  and  analysis;  first  welding 

pass) 
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Effects  of  residual  stresses 

Besides  the  basic  issues  of  structural  and  material  strength, 
several  additional  factors  have  to  be  taken  into  account  when 
one  wishes  to  analyze  the  service  behavior  of  a  welded  structure. 
One  of  the  most  important  of  these  factors  is  the  presence  of 
residual  stresses  due  to  welding. 

Unfortunately,  there  are  no  general  rules  that  can  guide  the 
designer  of  a  welded  structure  in  the  determination  of  the  ef¬ 
fects  these  residual  stresses  will  have  on  the  reliability  of  the 
structure.  This  is  mainly  due  to  the  fact  that  these  effects  can 
be  very  complex  in  nature.  Furthermore,  and  depending  on 
the  nature  of  the  structural  or  material  characteristic  one  is 
interested  in,  controversial  opinions  might  be  found  in  the 
available  literature  regarding  the  true  effect  the  residual  stresses 
have. 

In  this  section  an  effort  will  be  made  to  briefly  discuss  some 
aspects  of  the  consequences  of  residual  stresses  on  the  service 
behavior  of  welded  structures  in  a  more  or  less  qualitative 
manner.  This  discussion  will  thus  lead  to  the  more  specific 
analysis  of  the  fracture  characterization  of  weldments. 

Brittle  fracture 

It  is  widely  known  that  to  avoid  brittle  fracture  in  a  welded 
structure,  the  material  must  have  adequate  notch  toughness. 
Differently  stated,  and  using  concepts  from  the  fracture  me¬ 
chanics  theory  |1|,  unstable  fracture  occurs  when  stresses  are 
applied  to  a  structure  containing  a  crack  longer  than  a  given 
value  (critical  crack  length).  For  low-carbon  steel  this  critical 
crack  is  several  inches  long  at  the  yield  stress.  Brittle  fractures 
in  welded  structures  have  been  observed  to  originate,  however, 
from  relatively  short  cracks  and  under  an  overall  stress  equal 
to  about  one-third  of  the  material  yield  strength. 

These  catastrophic  failures  of  welded  structures  from  sub- 
critical  cracks,  called  low  applied-stress  fractures,  are  due  to 
presence  of  high  tensile  residual  stresses  in  areas  where  the 
cracks  are  located.  The  crack  can  thus  grow  uncontrollably 
even  though  the  level  of  applied  stress  is  low  ( 1 ). 

Fatigue  fracture 

Only  limited  studies  have  been  made  on  the  analysis  of  the 
effects  of  residual  stresses  on  the  fatigue  strength  of  welded 
structures.  This  stems  from  the  major  difficulties  that  come 
from  the  change  in  stress  distributions  around  a  crack  as  it  grows 
due  to  the  repeated  loading. 

Because  of  these  and  several  other  reasons  1 1  ],  our  knowledge 
of  the  relationship  between  residual  stress  and  fatigue  strength 
is  confused  An  observation  that  has  been  reported  by  several 
investigators,  although  not  conclusive,  is  that  the  fatigue 
strength  of  welded  specimens  seems  to  increase  when  com¬ 
pressive  residual  stresses  exist  in  regions  near  the  surface  of  a 
plate  It  should  lie  pointed  out,  however,  that  the  fatigue 
strength  depends  greatly  on  the  condition  of  the  surface;  the 
effect  of  residual  stress  is  thus  secondary  and  is  overshadowed 
by  such  major  factors  as  weld  geometry  and  surface  irregular¬ 
ities  1 1 1. 

Stress-corrosion  cracking  and  hydrogen 
embrittlement 

As  previously  mentioned,  residual  stresses  significantly  affect 
those  phenomena  that  occur  under  a  low  applied  stress.  Since 
stress-corrosion  cracking  and  hydrogen-induced  cracking  of 
a  weldment  has  Ireen  observed  to  occur  without  any  external 
loading,  it  can  be  concluded  that  residual  stresses,  especially  the 
high  tensile  residual  stresses,  do  significantly  affect  these  phe¬ 
nomena.  leading  to  catastrophic  failures  1 1  ]. 


Fig.  16  Longitudinal  stress  distribution  at  several  time  instances  (first 
welding  pass) 


Buckling  strength 

It  is  generally  known  that  residual  compressive  stresses  de¬ 
crease  the  buckling  strength  of  a  welded  structure,  as  measured 
for  example  by  the  externally  applied  critical  stress.  This  is  true 
for  columns,  plates,  stiffened  plate  structures,  spherical  and 
cylindrical  shells  1 1,20,38]. 

Methodologies  for  analyzing  residual  stress 
effects 

In  Part  1  of  this  paper,  and  in  particular  when  discussing 
methods  for  predicting  residual  stresses,  it  was  pointed  out  that, 
due  to  the  complex  nature  of  the  problem  and  before  the  wide 
use  of  the  electronic  digital  computer,  most  investigations  in 
the  area  were  done  experimentally.  This  statement  is  even 
more  true  in  the  case  of  analyzing  the  effects  of  residual  stresses 
on  the  service  behavior  of  welded  structures. 

Very  few  analytical  studies  have  actually  been  performed 
in  the  area.  Moreover,  these  studies  appear  to  be  rather  isolated 
attempts  at  addressing  the  various  problems.  A  more  sys¬ 
tematic  effort  has  been  followed  by  M.I.T.  investigators  aimed 
at  providing  designers  of  ships  and  other  welded  structures  with 
tools  capable  of  aiding  in  the  reliability  assessment  of  these 
structures. 

The  methods  used  by  M.I.T.  investigators  are  both  analytical 
and  numerical.  The  former  ones  can  be  used  whenever  a 
rougher  estimate  of  the  residual  stress  effects  is  desirable, 
whereas  numerical  techniques,  such  as  the  finite-element 
method,  should  be  used  when  more  detailed  and  accurate  es¬ 
timates  are  required  or  when  the  problem  is  too  complex  to 
analyze  otherwise. 

Examples  of  cases  that  have  already  Ireen  studied  using  an¬ 
alytical  techniques  include  the  corrugation  damage  of  ship 
bottom  plating  |20],  the  effect  of  residual  stresses  on  the  buek- 
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metal  and  welding  electrodes  through  tests  on  various  restrained 
welded  specimens.  The  results  of  these  tests  are  used  on  a  “go” 
or  “no  go”  basis,  which  means  that  when  the  restrained  joint 
can  be  welded  successfully  without  cracking,  then  the  material 
and  the  welding  procedures  used  in  the  test  can  also  be  used 
successfully  during  actual  construction.  Very  little  effort  has 
been  done  to  obtain  quantitative  information  on  the  fracture 
potential  of  such  weldments. 

A  study  was  thus  undertaken  at  M.I.T.  to  investigate  the  ef¬ 
fects  of  residual  stresses  on  some  of  the  fracture  characteristics 
of  welded  structures,  especially  those  made  of  high-strength 
quenched-and-tempered  steels  (3j.  The  methodology  used  and 
some  results  obtained  in  this  study  are  reported  next,  serving 
as  a  case  study  of  how  one  can  analyze  the  effects  of  residual 
stresses  on  the  service  behavior  of  a  welded  structure. 
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ling  strength  of  welded  plates  |38|,  the  systems  analysis  ap¬ 
proach  to  the  fatigue  crack  growth  of  surface  and  embedded 
transverse  cracks  in  a  butt  weld  in  the  pressure  hull  of  a  deep 
submersible  [  1 ),  and  others. 

On  the  other  hand,  the  more  sophisticated  numerical  tech¬ 
niques  have  been  used,  for  example,  in  the  fracture  analysis  on 
weldments.  This  analysis  will  be  the  topic  of  the  next  sec¬ 
tion. 


Fracture  analysis  of  weldments:  a  case  study 

The  pursuit  of  improved  technology  for  designing  high- 
performance  marine  vessels  has  led  to  the  use  of  new  materials 
having  high  strength-to-weight  ratios,  as  is  the  case  with  the 
H  Y-80  steel  used  in  present  day  U  S.  Navy  applications  or  the 
even  more  advanced  HY-130  steel. 

Unfortunately,  however,  these  high-strength  steels  can  be 
susceptible  to  catastrophic  failures  in  the  presence  of  small 
flaws,  making  the  initiation  and  propagation  of  cracks  in  the 
presence  of  residual  stress  fields  of  the  utmost  concern.  Pres¬ 
ently  the  U  S.  Navy  evaluates  the  cracking  sensitivity  of  the  base 


Tabic  4  Parameters  of  specimens  for  fracture  experiments 


Specimen 

Number 

Type  oi 
Steel 

Thickness 
mm  (in.) 

R 

mm  (in.) 

K, 

(MPa/mm) 

1 

SAE  1020 

12.7(1/2) 

305(12) 

32.5 

2 

SAK  1020 

2.r>.4  ( 1 ) 

457 (18) 

343, 

SAE  1020 

25.4  (D 

305(12) 

65.1 

4 

SAE  1020 

50.8(21 

457 (18) 

68.0 

5 

SAE  1020 

50.8(2) 

:’05(12l 

130.2 

6 

HY-130 

22.2 (7/8) 

457 (18) 

29.8 

7 

HY-130 

22.2(7/8) 

305(12) 

56.9 

8 

HY-130 

50.8(2) 

305(12) 

130.2 

Fracture  experiments 

As  previously  mentioned,  one’s  confidence  in  the  predictive 
capabilities  of  any  analysis  can  grow  only  if  the  analytical  results 
correlate  with  experimentally  obtained  data  in  a  satisfactory 
manner.  For  this  reason  a  series  of  experiments  was  performed 
and  later  analyzed  as  discussed  in  the  next  subsection. 

Figure  17  shows  the  typical  configuration  of  the  simple  re¬ 
strained  specimens  used  in  the  experimental  program.  A  total 
of  eight  similar  specimens  was  used  as  detailed  in  Table  4. 
Welding  was  performed  using  the  GMAW  process  according 
to  U.S.  Navy  specifications. 

Both  temperature  changes  and  transient  strains  were  mea¬ 
sured  during  the  welding  operation  using  thermocouples  and 
electric  resistance  strain  gages,  respectively,  at  preselected  lo¬ 
cations.  The  joint  transverse  shrinkage  was  also  measured  after 
each  welding  pass  using  a  clip  gage.  Results  of  these  mea¬ 
surements,  which  are  not  of  interest  in  the  present  discussion, 
can  be  found  in  (3,4,49). 

After  the  completion  of  welding,  the  crack-opening  dis¬ 
placement  (COD)  was  measured  in  all  specimens.  It  should 
be  noted  that  the  COD  concept  can  be  used  instead  of  the 
stress-intensity  factor  since  both  can  provide  similar  information 
regarding  the  fracture  characteristics  of  a  weldment. 

The  COD  was  recorded  while  a  notch,  inserted  in  a  direction 
parallel  to  the  weld  (x-direction)  and  located  on  the  weld  cen¬ 
terline,  was  continuously  propagating  in  the  specimen’s 
thickness  direction.  No  external  stresses  were  applied  during 
this  process;  the  only  stresses  present  were  the  residual  stresses 
due  to  welding.  Under  such  conditions  the  crack  would 
propagate  through  zones  of  relatively  constant  fracture 
toughness  (that  of  the  weld  metal)  under  the  influence  of  high 
values  of  transverse  residual  stresses,  rry,  which  can  be  consid¬ 
ered  to  be  independent  of  x  in  the  central  portion  of  the  spec¬ 
imen. 

Notches  were  made  with  0.8-mm-thick  ('/is  in  ),  15.25- 
cm-diameter  (6  in.)  circular  blade  saws.  The  saw  tips  were 
ground  to  a  90-deg  angle  so  that  a  sharp  notch  tip  could  be 
obtained.  The  COD  was  measured  using  a  specially  designed 
clip  gage  with  four  Type  SP- 133-20-35  semiconductor  strain 
gages  bonded  on  its  spring  steel  arms.  To  improve  accuracy 
and  sensitivity  of  measurements,  all  four  strain  gages  were  ac¬ 
tive  in  the  Wheatstone  bridge  arrangement 

Results— Curves  of  COD  versus  notch  depth  were  obtained 
for  all  specimens  tested.  The  relevant  curves  for  two  HY-130 
specimens  (Nos.  6  and  7  in  Table  4)  are  reproducer!  here  in  Fig. 
18. 

It  is  interesting  to  mention  that  the  relatively  high  values  of 
COD  shown  were  obtained  with  no  applied  external  stresses 
Furthermore,  the  nonlinear  character  of  the  curves  signifies  the 
presence  of  a  complex  through-thickness  residual  stress  distri¬ 
bution 
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Numerical  analysis 

Joint  degree  of  constraint — As  a  first  step  towards  analyzing 
the  fracture  experiments  on  the  simple  restrained  specimens, 
the  joint  degree  of  constraint  has  to  be  calculated.  This  is  de¬ 
fined  by 

K,  =  <70 16 

where  Oq  is  a  tensile  stress  uniformly  distributed  along  the  joint 
and  5  is  the  joint  transverse  dispacement  (for  complex  joints  an 
average  value  is  used).  K,  is  a  kind  of  spring  constant  of  the 
structure.  Its  value  can  thus  be  assigned  to  a  number  of  linear 
springs  which  model  the  restraint  imposed  on  the  central  part 
of  the  structure  by  the  surrounding  material. 

The  degrees  of  constraint  for  all  the  specimens  were  calcu¬ 
lated  using  a  simple  linear  elastic  finite-element  analysis  [3]. 
The  results  are  given  in  Table  4. 

Estimation  of  residual  stress  field — The  fracture  experi¬ 
ments  were  performed  in  the  absence  of  any  externally  applied 
stresses;  the  only  stresses  present  were  the  residual  stresses  due 
to  welding.  To  successfully  analyze  these  experiments  using 
linear  elastic  fracture  mechanics  (LEFM),  therefore,  it  is  nec¬ 
essary  to  know  the  residual  stress  field,  and  in  particular  the 
residual  stresses  perpendicular  to  the  crack  surfaces,  <Ty. 

One  way  of  calculating  these  residual  stresses  is  through  the 
thermal-elastic-plastic  finite-element  analysis  outlined  in  Part 
1  of  this  paper.  Such  an  analysis,  however,  of  all  eight  speci¬ 
mens  would  be  prohibitively  expensive.  Instead,  a  simplified 
approximate  approach  was  used  for  determining  the  surface 
and  through-thickness  distribution  of  transverse  residual 
stresses. 

This  approximate  method  is  based  on  the  fact  that  the  COD 
measured  during  the  fracture  experiment  was  due  to  the  elastic 
release  of  the  residual  stresses  present  around  the  notch.  Thus 
by  applying  the  superposition  principle  the  distribution  of  the 
transverse  residual  stresses  can  be  obtained  by  finding  the  stress 
which  if  applied  on  the  notch  surface  would  result  in  the  same 
COD  values  as  the  experimental  ones  [3|. 

A  finite-element  code,  utilizing  several  subroutines  of  the 
finite-element  program  FEABL  [50|.  was  developed  to  im¬ 
plement  the  method. 

Figure  19  shows  the  finite-element  mesh  used,  where  only 
triangular  and  quadrilateral  assumed  displacement  elements 
were  used.  From  symmetry  considerations  only  one-half  of 
each  specimen  was  considered;  the  restraint  was  simulated  by 
equivalent  elastic  springs. 

The  notch  depth  increase  was  simulated  by  16  incremental 
steps,  each  equal  to  ‘/ih  of  the  plate's  thickness.  Steps  1  and  10 
are  shown  in  f  ig  19.  In  each  step  the  nodal  force  distribution 
that  would  give  the  same  COD  as  the  experimentally  obtained 
one  was  calculated  iteratively  using  the  method  of  false  posi¬ 
tion 

Figure  20  shows  the  calculated  through-thickness  residual 
stress  distribution.  rr„,  for  two  HY-130  specimens  (Nos.  band 

7). 

Calculation  of  stress-intensity  factor— To  calculate  the 
stress- intensity  factor  for  the  specimens  used  in  the  fracture 
experiments  a  linear  elastic  finite-element  analysis  was  per¬ 
formed  Instead  of  using  standard  finite  elements,  however, 
something  that  would  require  up  to  1500  degrees  of  freedom 
to  obtain  reasonable  accuracy,  a  hybrid  crack  element  was 
developer!  enabling  one  to  simulate  the  crack  growth  in  a  much 
more  inexpensive  way  The  element  developed  is  an  extension 
of  the  one  proposed  by  Tong  et  al  [51 1  It  introduces  stresses 
on  the*  crack  surface  to  simulate  the  residual  stress  field,  whereas 
the*  original  element  assumed  stress-free  crack  surfaces  Details 
on  the  development  of  this  element  are  not  included  here  The 
interested  reader  is  referred  to  |3,4j. 
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The  hybrid  crack  element  was  used  in  conjunction  with 
FEABL  to  perform  a  linear  finite-element  fracture  analysis  of 
the  experiments  described  earlier. 

Figure  21  shows  the  finite-element  mesh  used  in  the  analysis. 
Only  the  central  part  of  a  specimen  is  represented  in  the  mesh; 
a  series  of  springs  applied  on  both  sides  of  the  mesh  simulates 
the  restraint  the  rest  of  the  specimen  exhibits  on  this  central 
portion  of  the  plate. 

The  crack  growth  was  simulated  by  two  ways:  first,  by  in¬ 
creasing  the  size  of  the  crack  in  the  hybrid  crack  element 
without  changing  the  rest  of  the  mesh;  and  second,  by  moving 
the  hybrid  crack  element  deep  into  the  mesh,  something  that 
required  changes  in  the  position  of  some  of  the  other  elements, 
increases  in  the  number  of  nodes,  and  a  renumbering  of  ap¬ 
proximately  half  of  the  nodes  in  the  mesh.  The  crack  growth 
was  done  automatically  by  the  program  in  increments  of  Vib 
of  the  plate  thickness,  I,  from  VirI  to  14/i«t.  Figure  22  presents 
schematically  several  of  these  crack  growth  steps. 

Equivalent  nodal  forces  and  distributed  stresses  on  the  crack 
surface  were  determined  and  prescribed  automatically  by  the 
program  at  each  step. 

Numerical  results  of  the  Mode  I  stress-intensity  factor,  K|, 
obtained  by  the  finite-element  analysis,  are  presented  in  Fig. 
23  for  Specimens  6  and  7,  where  K i  is  plotted  versus  the  crack 
depth,  a. 

It  is  interesting  to  note  the  nonlinear  character  of  the  curves, 
which  usuallv  show  high  values  of  K|  for  cracks  with  depths 
equal  to  about  one  quarter  of  the  plate’s  thickness;  some  of  the 
curves  even  show  negative  values  for  K |. 
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Fig.  19  Finite-element  mesh  used  for  the  approximate  calculation 
of  the  through-thickness  residual  stress  distribution 


RESIDUAL  STRESS  (MPo) 


The  calculated  values  of  the  crack-opening  displacement  at 
a  distance  6.35  mm  (0.25  in.)  from  the  crack  centerline  are 
compared  with  the  experimentally  obtained  ones  in  Fig.  18  for 
the  same  specimens.  Close  agreement  between  the  two  is  ob¬ 
served  for  the  case  of  small  cracks.  The  divergence  between 
the  numerical  and  experimental  curves  observed  for  longer 
cracks  is  probably  due  to  three-dimensional  effects  as  the  crack 
grows  deeper. 


Fig.  20  Through-thickness  distribution  of  transverse  residual  stress, 
tr y  (Specimens  6  and  7) 

Summary,  conclusions  and  future  possibilities 

This  paper  presents  computer-aided  methodologies  for 
solving  the  problems  of  predicting  residual  stresses  due  to 
welding  and  of  analyzing  their  consequences  on  the  service 
behavior  of  welded  structures.  These  methodologies  represent 
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Fig.  22  Schematic  representation  of  crack  growth  process 


the  efforts  of  M.I.T.  investigators  over  a  period  of  more  than 
ten  years  in  this  area  and  can  be  thought  of  as  the  state  of  the 
art  on  a  subject  that  is  very  complex  but  at  the  same  time  very 
important  for  the  naval  architectural  community. 

Part  1  of  the  paper  discusses  the  subject  of  residual  stress 
prediction.  General  information  is  first  presented,  including 
how  residual  stresses  can  be  classified,  how  they  are  formed 
during  the  welding  operation,  and  which  are  the  physical 
mechanisms  that  should  be  accounted  for  in  any  accurate 
analysis  capable  of  predicting  them.  This  is  followed  by  a 
historical  overview  of  the  different  methodologies  developed 
for  their  analysis.  It  is  concluded  that  the  most  accurate 
analysis  is  the  one  based  on  the  simulation  of  the  whole  welding 
operation  using  a  computer. 

The  welding  simulation  consists  of  first  calculating  the 
temporal  and  spatial  temperature  distributions  resulting  from 
the  application  of  the  intense  welding  heat  source,  and  subse¬ 
quently  of  using  this  information  as  one  of  the  inputs  in  a 
thermal-elastic-plastic  stress  analysis.  Both  analytical  and 
numerical  methods  for  performing  each  of  these  two  steps  are 
presented  These  are  then  investigated  with  respect  to  their 
predictive  capabilities  and  their  limitations  by  comparing  re¬ 
sults  from  sample  analyses  with  experimentally  obtained 
data. 

It  is  concluded  that  closed-form  analytical  solutions  and 
simple  numerical  ones  (that  is,  the  one-dimensional  computer 
program)  can  be  used  with  reasonable  degree  of  confidence  in 
the  analysis  of  thin  plates  welded  by  conventional  processes  (for 
example,  GMAW,  GTAW)  and  in  the  analysis  of  medium¬ 
thickness  plates  welded  by  high-energy-density  processes  (for 
example,  electron  beam  and  laser).  On  the  other  hand,  if  high 
accuracy  is  required  (that  is,  when  fabricating  critical  structures 
such  as  deep-sea  submersibles)  or  when  analysis  of  thick 
structures  is  desired,  the  only  currently  available  technique  that 
can  be  used  is  a  numerical  one  based  on  the  finite-element 
method. 

Following  completed  discussion  of  the  problem  of  residual 
stress  prediction,  Part  2  of  the  paper  concentrates  on  the  con¬ 
sequences  residual  stresses  have  on  the  service  behavior  of 
welded  structures.  This  is  a  very  important  problem  since  it 
is  well  known  that  residual  stresses  adversely  affect,  among 
others,  the  brittle  fracture,  fatigue  life,  buckling,  stress  corrosion 
cracking,  and  hydrogen  embrittlement  characteristics  of 
weldments. 

The  vastness  of  this  subject  area  does  not  allow  for  a  detailed 


description  of  all  the  methods  that  can  be  used  for  analyzing 
the  various  problems  encountered.  For  this  reason,  after  a  brief 
qualitative  general  discussion,  the  emphasis  of  Part  2  is  placed 
on  a  particular  subject,  that  of  the  quantitative  characterization 
of  the  fracture  of  welded  high-strength  quenched-and-tem- 
pered  steels.  A  linear  elastic  fracture  mechanics  analysis  using 
the  finite-element  method  is  first  developed.  Its  predictions 
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Fig.  23  Stress  intensity  factor  versus  crack  depth 
(Specimens  6  and  7) 
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are  then  compared  with  experimentally  obtained  data. 

A  major  conclusion  can  be  drawn  from  the  discussion  pre¬ 
sented  in  this  paper.  Among  the  various  disciplines  in  naval 
architecture,  including  hydrodynamics,  structural  mechanics, 
and  dynamics,  the  technology  of  welding  fabrication  has  long 
been  primarily  empirical.  One  of  the  major  reasons  for  this 
may  be  the  complex,  transient  phenomena  involved  during 
welding.  In  the  authors’  opinion,  however,  computer-aided 
analyses  like  the  ones  presented  in  this  paper  prove  that  it  is 
possible  for  the  designers  and  fabricators  to  start  incorporating 
manufacturing-related  parameters  in  their  analyses.  It  is  hoped 
that  such  practice  can  start  soon  because  it  makes  a  lot  of  sense 
economically. 

Future  research  efforts 

The  question  that  one  might  ask  is,  Where  do  we  go  from 
here?  Although  there  is  no  single  answer,  it  is  the  authors’ 
belief  that  very  exciting  opportunities  lie  ahead. 

For  one,  the  fact  that  an  analysis  of  such  a  complex  phe¬ 
nomenon  as  welding  has  been  proven  to  be  feasible  leads  to  the 
possibility  of  analyzing  several  other  welding-related  or  more 
general  fabrication  processes.  The  flame  heating  and  flame 
bending  operations,  the  thermal  stress  relieving  process,  and 
the  problem  of  weld  metal  overmatching  or  undermatching 
are  just  a  few  that  come  to  mind.  In  addition,  it  will  be  well 
worth  analyzing  several  of  the  quality  acceptance  tests  currently 
used  and  which  are  mostly  empirical  in  nature;  some  of  them 
do  not  even  provide  any  quantitative  information  but  are  rather 
performed  on  a  "go”  or  "no  go”  basis,  such  as  the  explosion 
bulge  test  or  the  restrained  cracking  window  test.  There  is  a 
high  probability  that  these  analyses  will  aid  in  the  rationaliza¬ 
tion  of  the  tests  or  even  in  the  development  of  new,  more  sci¬ 
entifically  meaningful  ones.  We  believe  that  the  ever-in¬ 
creasing  power  and  speed  of  the  large  digital  computers  that 
is  spurred  by  such  developments  as  parallel  processing  will 
enable  analyses  of  this  type  to  become  more  economical  and 
hence  more  widespread 

On  the  other  hand,  the  tremendous  increase  in  the  power  of 
microprocessors  and  microcomputers,  referred  to  by  many  as 
the  second  industrial  revolution,  will  S(x>n  enable  a  single,  in¬ 
expensive  chip  to  perform  calculations  based  on  some  of  the 
simple  models  descrilred  in  this  paper  If  one  couples  this  with 
appropriate  sensors,  it  could  lx-  possible  to  develop  welding 
machines  with  fully  adaptive  automatic  control  that  can  always 
perform  perfect  welds  in  the  presence  of  any  external  distur¬ 
bances  Such  machines,  the  dream  of  all  fabricators,  will  lx* 
[tosxible  some  time  in  the  future  if  enough  effort  is  applied  in 
this  direction. 
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Appendix  1 

Closed-form  solution  of  three-dimensional  fir 
heat  source  model 

The  solution  to  the  governing  partial  differential  equa 
of  heat  transfer  during  welding  under  the  assumptions  discusseu 
in  the  main  text  was  found  using  the  method  of  separation  of 
variables,  the  appropriate  boundary  conditions,  and  by  dividing 
the  plate  into  two  regions,  inside  and  outside  the  heat  input 
circle  of  radius  r/,  [2],  The  obtained  solution  can  be  written 
as 

0  =  do  +  ~  •  (\/l  +  2>  •  e~Xvi  ■  d>(r,z)  —  1) 
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hi  and  /»2  the  heat-transfer  coefficients  at  the  top  and  bottom 
surfaces  of  the  plate,  respectively;  H  the  plate  thickness;  J o(x) 
and  Jj(x)  the  Bessel  functions  of  first  kind  and  of  zero  and  first 
order,  respectively;  I o(x)  and  / i(i)  the  modified  Bessel  functions 
of  first  kind  and  of  zero  and  first  order,  respectively;  and  Ko(x) 
and  Ki(x)  the  modified  Bessel  functions  of  second  kind  and  of 
zero  and  first  order,  respectively.  All  other  variables  have  been 
previously  defined. 

A  computer  program  has  been  written  in  the  FORTRAN  IV 
language  to  perform  the  necessary  calculations. 
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